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The backseattering of fast (0.3, 1, 3, 6, 9, and 14 Mev) neutrons from a semi-infinite 
water medium has been calculated by the Monte Carlo method. The information obtained 
includes the joint angular and spectral distribution of the reflected neutrons, the dependence 
of the number albedo and energy albedo on the source energy and obliquity, and the con- 
tributions to the albedo of successive orders of seattering. The spectra were calculated 
down to epithermal energies (~0.5 ev). The results for each case are based on the analysis 
of 3,000 neutron histories, generated with the use of an IB M-—704 computer. In the random 
sampling, elastic scattering from hydrogen and oxygen, inelastic scattering from oxygen, 
and absorption due to n-a@ and n-p processes were taken into account. The cross sections 
for some of these interactions are not well known. Parallel calculations with different 
assumptions about the cross sections were made in order to estimate how sensitively the 
albedo depends on the cross sections. The paper inciudes a self-contained description of 
the Monte Carlo method, its application to the calculation of radiation diffusion and in 
particular to the neutron albedo problem. Emphasis is placed on the technique of cor- 
related sampling which makes possible an accurate estimate of albedo differences resulting 
from different assumptions about the cross sections. The random sampling computations 
were supplemented by analytical calculations of the single-scattering albedo. This was 
useful for an understanding of the Monte Carlo results because a considerable fraction of 
the reflected neutrons return after only one collision. 


1. Introduction 
1.1. Statement of the Problem 


Neutrons that are incident on an extended medium make one or more collisions with 
nuclei in each of which they are deflected and lose a fraction of their energy. Thus they 
perform a “random walk” until they are either absorbed or re-emerge from the medium with 
reduced energy. The probability that an individual neutron will eventually re-emerge from a 
semi-infinite medium (i.e., be reflected) is called the albedo, or the differential albedo if the 
energy and/or direction of the reflected neutron are specified. 

The purpose of this paper is threefold: (1) To calculate the differential albedo and related 
quantities for neutrons incident on a semi-infinite water medium; (2) to give a self-contained 
exposition of the methods of calculation used (Monte Carlo, single-scattering analysis); (3) 
to analyze how the albedo depends on various factors, particularly the cross sections for elastic 
and inelastic seattering and absorption. 


1.2. Background 


Exact calculations of the neutron albedo have so far been limited to situations in which 
the slowing down of neutrons could be disregarded. With this limitation, and with further 
restrictive assumptions, Halpern, Lueneburg, and Clark [1]? have given an analytic solution 
based on transport theory. A more general solution of the one-velocity problem is also avail- 
able from Chandrasekhar’s [2] theory of radiative transfer. Fermi [3] and more recently 


Grosjean [4] have given approximate solutions which are in good agreement with the more 


rigorous calculations. 
When the complicated energy-dependence and directional dependence of the collision 
processes must be taken into account, an exact analytical treatment becomes all but impossible. 
1 Supported jointly by Armed Forces Special Weapons Project, OCDM, and the Department of the Navy (Bureau of Yards and Docks), 
2 Figures in brackets indicate the literature references at the end of this paper. 








Spinney [5] has given an approximate treatment based on age theory which is applicable 
mainly to energies lower than those considered in the present work.’ Another possible approach 
would be to incorporate the exact one-velocity treatment of reference [2] into an energy-group 
scheme similar to the multigroup diffusion schemes used in reactor calculations. This has not 
been attempted so far. 

The Monte Carlo method recommends itself as the simplest and most promising approach. 
It has the advantage that the complexity of the calculations is not essentially increased by 
adding to the number of variables or by introducing complicated cross sections. Increased 
detail merely creates a data-handling problem which requires a computer with sufficient ca- 
pacity. A disadvantage of the Monte Carlo method is its purely numerical character which 
makes necessary a very careful scrutiny of the results before one can understand their signifi- 
We have found that such an analysis is greatly facilitated by parallel analytical 
In spite of its promise, the Monte Carlo method has not 
Some results for the backseattering of 


cance. 
calculations of single scattering. 
vet been used extensively for albedo calculations. 
0.89-Mev neutrons from water have been given by Foderaro and Obenshain [6], but their 
work was mainly concerned with transmission through slabs. 

One reason for the searcity of Monte Carlo applications to the albedo problem appears to 
be the lack of adequate information on neutron collision cross sections, which may have dis- 
For few if any elements (other than hydrogen) are the 


couraged extensive numerical work.* 
Only the knowledge of 


cross sections known in the detail and with the accuracy desirable. 
the total collision cross section is in definitive form for most elements. Information about 
the contributions of different processes to the total cross section is reasonably adequate. Much 
less is known about differential cross sections. The angular distribution of elastic scattering 
which is highly anisotropic at high energies—is not well known for many elements. 


deflections 
There is uncertainty as to the nuclear levels 


Even less is known about inelastic scattering. 
which are excited, the relative excitation probabilities, and the angular distribution of the 
inelastic scattering deflections. 

We have chosen water as the medium in our calculations because the cross sections for 
hydrogen and oxygen are known better than those of most other elements. In addition to 
being of interest in their own right, the results for water will also be typical of those expected 
with other hydrogeneous substances. Thus the present work should provide guidance for 
future calculations for media such as concrete which are held up now by lack of cross section 
data. 


2. Variables of the Problem 


We consider the reflection of neutrons from a semi-infinite medium occupying the region 
z>0. Figure 1 shows a typical trajectory of a neutron which is incident on the medium with 
obliquity 6 and which is reflected with obliquity 6. These obliquities are angles between the 
direction of motion of the neutron and the positive z-axis. Thus 4 lies between 0 and 7/2, 
and @ lies between 2/2 and z. 

We denote by A(E,6;/),%) sinédédE the probability that a neutron incident with energy 
FE, and obliquity 4 will be reflected from the semi-infinite medium, after one or more scatterings, 
with an energy between F and /+-dE and an obliquity between 6 and 6+-dé. The function 
A(E,0;£,0)) is called the differential albedo. Alternatively, instead of regarding the albedo as 
a probability, we may consider it to be the ratio of the average number of neutrons reflected 
to the number of neutrons incident. Thus if a beam of Np neutrons is incident with energy 
Fi) and obliquity 6, and if N(4,@) sinédédE neutrons are reflected on the average, then 


4 


N(E,6) (2.1) 


A(E,0;E,,0) =~ > . 
040) No(£20,00) 


3 See, however, sec. 7.8. 
4 For a survey of the neutron cross section situation in relation to transport calculations, see the recent book by Goldstein [7]. 
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By integrating the differential albedo one obtains the number albedo Ay (ratio of the number 
of all reflected neutrons to the number of incident neutrons): 


4 F Eo ,’ aj . , , ‘ 
Ary(E',9o) | dE q sindd0A(E,0;E,0). (2.2) 
JD) 


J 1/2 


Similarly, the energy albedo Ag (ratio of the average reflected to the incident energy) is defined as 


; ee ee orn 
Az(Eo,o)- | ik | sinéd@ = A(E6; Eo ,00). (2.3) 


0 1/2 Ky 


We shall have occasion to consider the probabilities of reflection after exactly 1, 2,.. . 
collisions. The corresponding albedos will be indicated as A (E, 6;/5,0)), Ay”? (45,0) and 
Ap? (1,9), n=1,2, .... The case n= 1 (single-scattering albedo) is particulary important. 


9 a9 


Figure 1. Neutron trajectory. 








3. Neutron Cross Sections of Hydrogen and Oxygen 
3.1. Integral Cross Sections 


The most important parameter characterizing the neutron transport process is the mean 
free path between collisions, \(/2). It is inversely proportional to the number -‘% of target 
nuclei per unit volume and the total collision cross section o(F), 


] 
\(E)=— Tex (3.1) 
) ANo(E) ns) 
The reciprocal of the mean-free path (the attenuation coefficient) is additive for a homogeneous 


mixture of different nuclei. For water, 


l l l 


tes hy ia lou ( 12) { ‘ (ox90x(E), (3.2) 


where . [76.69 X10”? em~* and AG, =3.34 X10” em~?, 
We have taken the values of o4(/2) and oox(/) from the standard Brookhaven compilation 
[8]. The attenuation coefficients for hydrogen and oxygen thus obtained are shown separately in 
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figure 2 as functions of the energy. The plots are made in histogram form corresponding to the 
approximate manner in which the data were used in the machine computations. It can be seen 
that the energy-dependence for hydrogen is smooth, but that the histogram for oxygen has 
several sharp peaks and irregularities at high energies. 

Next we consider the relative contributions of various processes to the total collision cross 
section in the energy region of interest (0.5 ev to 14.0 Mev). For hydrogen, only elastic seat- 
tering is important in this region. The (n,y) cross section is smaller than 0.1 percent of the total 


cross section at 0.5 ev, and is even less significant at higher energies. 


60 F 
J 


4 20 
55 + | > | 
1 | 

2 


+ 
+ HYDROGEN t 


——+—_——} 


ao Oo nN@wOo 


> 
+ 


wot 1g 


25 + 
HYDROGEN 


T ed ; OXYGEN 


.@) 5 10 5 2.0 25 3.0 35 40 
\ 
09 E(Mev) 


FiGuRE 2. Attenuation coe fhicients for hudrogen and orygen, 
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For oxygen, the processes that matter are elastic scattering, inelastic scattering (above 
~6.4 Mev), and absorption due to (7,@) reactions (above ~4 Mev) and (n,p) reactions (above = 
10 Mev). Cross sections for these processes, derived from the best experimental evidence and 
theoretical considerations, have been compiled by Lustig, Kalos, and Goldstein [9].° In table 
3.1 we list the values of the various cross sections at selected energies to indicate their relative 
importance. For use in machine calculations, we have fitted cross-section ratios with poly- 
nomials in the energy variable. The representations are a smoothed-out version of the data of 
reference [9 which is sufficient considering the uncertainty in 





accurate to 5 percent or better, 


the data. 


TABLE 3.1. Neutron cross sections (barns) for oxygen (from reference [9]) 
———- | 
Energy Cee "ae 76s C+ Ox 
Met Total Na n,p Inelastic 
14 1. 58 0. 31 0.09 0.59 
y l ts os 1 
6 1.48 19 
3 1.30 
Ratio of inelastic to total collision cross section: : 


On". OX___1_05173-+.0.24196E—0.013666 /2-+-0.00024009 13(6.4 Mev</<18 Mev) 


Fox 
Ratio of absorption to total collision cross section: 


Fa, Ox TF 


POX — — (). 404334-0.132544—0.0085614/2?+-0.000174133(4 Mev< H<18 Mev). 


J 





Tox 


The probabilities that a given collision will be with a hydrogen or an oxygen nucleus are 
Ano(/2) Ay (Ee) and Ay o(2)/Aox(/2), respectively. A collision can have four possible outcomes: 


1. With probability 
Muto 


. (3.4 
Pr he ? ) 


the neutron is scattered elastically by a hydrogen nucleus. 
2. With probability 


Ano 0oOx On! Ox Ta,Oox Tn. Ox 
= - - Pg 


) ff) 
Pe (9.0) 
Nox Fox 
the neutron is scattered elastically by an oxvgen nucleus. 
3. With probability 
Att,0 On! Ox 9a 
V3 — , (3.6) 
Nox Fox 
the neutron is scattered inelastically by an oxvgen nucleus. 
4. With probability 
Ant,o OF a,Ox Tt Tp, Ox "7 
Ds = 5) (3.7) 


Nox Fox 


the neutron is absorbed. Note that Pit Pot Pst Pr=l. 


5 Due to theoretical and experimental uncertainties, the cross sections given in this publication are subject to revision, We understand that 


Dr. Goldstein and his associates are now engaged in making an improved evaluation, 
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3.2. Differential Cross Sections 


The outcome of a collision in which a neutron of energy / is scattered can be specified in 
terms of the probability distribution /;(0;/) of the scattering deflections 0, in the center-of- 
mass coordinate system. In conformity with the notation in the preceding subsection, we 
shall use indices 7=1, 2, and 3 to indicate the distributions for elastic scattering from hydrogen, 
elastic scattering from oxygen, and inelastic scattering from oxygen, respectively. 

For elastic scattering from hydrogen, the distribution is assumed to be isotropic,® L.e., 


fi(O;F)=1. (3.8) 


The angular distribution for elastic scattering from oxvgen is also isotropic at energies 
P P Pee ) SD 
below about 300 kev, but markedly anisotropie at higher energies, particularly in resonance 

. tn) DP . 
regions. Above 300 kev we have used the compilation of reference [9] which gives the coef- 
ficients f, (2) of the Legendre expansion 


: —, 2/+1 ' 
f.(0;E)=D5 —5— Jo (E)P (cos @,). (3.9) 
i] “ 


(Note that the normalization is such that /:,0(/2)=1.) Typical angular distributions at various 
energies are shown in figure 3. 

The angular distribution for inelastic scattering from oxygen depends not only on the 
energy of the neutron before the collision but also on the excitation energy which the struck 
nucleus receives. In O'®, the four lowest levels are 6.06, 6.14, 6.91, and 7.12 Mev above the 
ground state, the next is at 8.6 Mev and there are numerous levels between 9 and 14 Mev [11]. 
It appears from the discussion in reference [9] that neither the relative probabilities for the 
excitation of these levels nor the corresponding angular distributions are known. In order 
to explore the effect of inelastic scattering on the albedo, we have done three sets of calculations 
on the basis of the following assumptions: 

(1) The nucleus is excited to a level 6.1 Mev above the ground state. The angular dis- 
tribution is isotropic; L.e., 
f3(0.;£)=1. (3.10) 

(2) The neutron loses all its energy in an inelastic scattering, so that it is, in effect, 
absorbed. 

(3) An inelastic scattering is treated as if it were an clastic scattering; Le., the excitation 
energy of the nucleus is set equal to zero. 

The assumptions of case (1) seem to us to be the most realistic, and may be expected to 
yield reasonable answers when the neutron source energy is not too much above 9 Mev. The 
assumptions of case (2) lead to an underestimate of the albedo, and the assumptions of case 
(3) lead to an overestimate of the energy of the reflected neutrons. 


4. Single-Scattering Albedo 


4.1. Formulation 
The probability that a neutron, incident on a semi-infinite medium with energy E) and 
obliquity 4, will be reflected with energy / and obliquity @ after exactly one collision is 


’ 7 er, ls 3 , ’ , e f “ 
A (E60 ;E,90): -| ens = 7 pT, (E,0;E,0)e-8' if r/2<0<n, 
0 No j=l (4 1) 


0 if 0<0<m/2. 
See Goldstein, Zweifel and Foster [10]. 


6 This assumption is open to doubt at energies above, say, 5 Mev. 
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Figure 3. Angular distribution of elastic scattering from oxygen (according to ref. [9)). 


. . ds “re . . . 
In this expression, the factor te =e represents the probability that the incident neutron will 
0 


travel a pathlength s and then undergo its first collision between s and s+ds, \=A(/) being 
the mean-free path. The collision can be of type 1 (elastic scattering from hydrogen), type 2 
(elastic scattering from oxygen), or type 3 (inelastic scattering from oxygen). The fourth 
type (absorption) need not be considered. The factor p; (defined by (3.4, 3.5, and 3.6)) 
represents the probability that the collision will be of type 7. The function H;(4,0;£5,9) is 
called the scattering function and denotes the probability that in a collision of type 7 the energy 
of the neutron will change from /, to Z and the obliquity from 6) to 6. Finally, exp(—s’/A) 
is the probability that the neutron will travel from the point of collision to the boundary 
without making any more collisions (A is the mean free path at the energy ). The pathlength 
to the boundary is 

P _ 2 098 Bo. (4.2) 

cos 6 
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Using this relation, and carrying out the integration with respect to s in (4.1), we obtain the 
equation 
wh Pj (hy) (E26; 15,00) cos 8 


— 


A‘? (EO; Fo,00) 
1 COS O—(Ao/A) COS Ao 


if r/2<0<7n, 


0 if 0<0<72. 


4.2. Scattering Function 


The function /7,(/2,4;45,9)) to be inserted into (4.3) can be derived from the corresponding 
distribution /,(0,.;/)) for the deflections in the center-of-mass coordinate system. To establish 
the connection, we need the relations between 0,, the corresponding deflection angle O in the 
laboratory coordinate svstem, and the ratio of energies before and after the collision, 4/2. 
These relations are a consequence of the conservation of momentum and energy in the collision 
of a neutron with a nucleus. We state them for the case of inelastic scattering, characterized 
by the transfer of an amount ©; of energy to the target nucleus as excitation.’ The same 


relations also apply to elastic scattering if we set “;—0. 
The energy of the neutron after the collision is given by 


1+-2M,K,cos 0.4 A?K:? 


. y (4.4) 
E=E 1+2' +A 
and the deflection in the laboratory svstem Is 
fk, 1+-M,K,;cos 0. M,-+1 keH M?Kk?-—1 Ki 7 
cos O - ae le = LS) 
VE 1+, 2 WE, 2(1,+1) VE 
where M, is the mass of the target nucleus and where 
; Q,1+-M 
k = | ——: - . (4.6) 
tg ee 


Let F( hk dk denote the probability that the energy of the neutron after the collision 
lies between Land k-+dhk. Using (4.4) we have 


Pees , 10, (M,+1 

F (Ey) =f(03E) sin 0, a= 545 ps S(O5F)- (4.7) 

lai “4e IM KE," " 
According to (4.5), the deflection 0 is determined by the value of / /,. The joint distribution 


of and 0 can, therefore, be expressed formally with the use of a Dirac delta funetion, as 


P(E) ( cos 0— 


$.S) 


Mtl /F ,MjKj-1 /K ) 
2 WS, ' 24,4) WE) 


Let the direction of the neutron before the collision be represented by a unit vector 
(sIn@, Cos go. SIN4, SING), cos#), and its direction after the collision I ao unit vector (sind 


COSY, sing sing, COSG). Then the cosine of the deflection angle is the inner product of the 
two vectors,” 

COSO = COS, COSA SINB, SINP COS( GQ) — @). (4.9) 
(The angles g, and g are azimuths with respect to the w-: plane.) If we insert (4.9) into 


(4.8), introduce the abbreviation 


‘OS @ : .— og (. ()) 
” 2 VE 2M4DVE -_ 


’ For a derivation, see, e.g , Bethe [12] 


§ See also figure 5 
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, 


and carry out an azimuthal average with respect to g’=g—¢, we obtain the seattering 
function 


; ; oe oar , ‘ 
HT (0; £0,600) =F, E; Eo) | de’ (cos I COS 8+-sin 4 sin 6 Cos gy’ —Ccos a;) 
us 


J 0 


: PF (EE) Paes (4.41) 


7. COS(a@;— 69) — COS@ | >! ¢osd COS(@;-+ 4p) 


The square roots in (4.11) must not become imaginary, which requires that 
COS(@;-+- 4) <cosd <cos(a;—A%). (4.12) 


In the case of perpendicular incidence (@)=0), the region (4.12) shrinks to a single value 


cos a), and the scatterme function becomes 
IT ,(£2,6:4,,0) I’ (1;£))6(cos 6—cos a@;). (4.13) 
In the limiting case of grazing incidence (@=90°), the scattering funetion becomes 


F (EE) 


mT sin’ @—cos’ a; 


HE, 0; Eo, 7/2) (4.14) 


4.3. Spectra and Angular Distributions 
When using the single-seattering results as a guide for interpreting the Monte Carlo 
results, it is useful to have the differential albedo with one of the variables integrated out; 
Le., the energy spectrum 


A” (E0; Fo,40) | sinddo.A' (12,0; 5,40), (4.15) 


and the angular distribution 


A (9: E69) ea" (1,0: y,0). (4.16) 

In the evaluation of these expressions the range of integration must include® all values 

of @ for given /¢ (or of / for given 6) for which the condition (4.12) is satisfied. Moreover, 

cosé must be negative. The resultant formulas are somewhat complex. We give them in 
appendix A, confining ourselves here to some qualitative remarks. 

The chief characteristic of the energy spectrum is that it is confined to limits which are 
wide for the case of hydrogen but very narrow for somewhat heavier nuclei such as oxygen. 
The upper limit of the spectrum follows from the condition that the neutron must be deflected 
sufficiently to aequire a velocity component in the negative z-direction. This implies, as a 


consequence of (4.12), that /’ cannot exceed 


kK... K, \ APA? cos? 09+ sin A . (4.17) 
' M41 


The lower limit of the spectrum is determined by the conservation of momentum and energy 


which implies, according to (4.4), that 7 cannot be smaller than 


| Bop? | | ). (4.18) 
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The allowed energy regions between /,,, and /yyin, are indicated in figure 4 as functions of the 
obliquity of incidence 4%. 

The angular distribution cannot be described by any one simple law. For perpendicular 
incidence, the distribution follows approximately a cosine law; 1.e., it peaks in directions corre- 
sponding to perpendicular emergence from the medium. As the direction of incidence becomes 
more oblique, this peak tends to disappear. In the limiting case of grazing incidence the angular 
distribution becomes flat or may even peak in directions corresponding to emergence parallel 
to the boundary of the semi-infinite medium. The angular distribution under these circum- 
stances also becomes a sensitive function of the distribution /(@,;/) of scattering deflections. 


For details, see appendix A. 
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Figure 4. Possible energies which a neutron may have if it is reflected after a single scattering. 


The allowed energy regions for elastic scattering hold for all Ho, whereas the indicated region for inelastie scattering applies only to Ho=14 


Mev and 6.1-Mevy excitation. 


5. Monte Carlo Method 
5.1. Generalities 


Although the general principles of the Monte Carlo method are well known,? we give a 
brief review to make this paper self-contained. The usual assumption is made that a neutron 
undergoing multiple scattering travels ina straight line until it makes a collision with a nucleus, 
is then suddenly deflected, travels again in a straight line till the next collision, and so on. 
The trajectory of a neutron is, therefore, completely described by a specification of the initial 
conditions and of the energies, obliquities, and positions of the neutron immediately after 
successive collisions. In the physical process, these characteristics are chance variables whose 
probability distributions are determined by the neutron cross sections. The Monte Carlo 
method consists of imitating nature by plaving a game of chance. With the use of random 
numbers, the parameters of model neutron trajectories are calculated according to the appro- 
priate probability distributions. Such a procedure is called random sampling, and a neutron 
trajectory thus generated is usually called a neutron history. 


9 See. e 
see, € 


g., reference [13] for background and bibliography. 
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By examining a sufficiently large number of neutron histories, one can obtain statistical 
estimates of various transport phenomena. For example, if No histories are generated in each 
of which the neutron starts at the boundary of a semi-infinite medium, and if in N of these 
histories the neutron eventually emerges from the medium, then the ratio N/N>o is an estimate of 
the number albedo Ay. 


5.2. Random Sampling of Neutron Histories 


If we let /,, 6,, and 2, represent the energy, obliquity, and distance from the boundary 
immediately after the w’th collision, the neutron trajectory or history can be described by the 
array 


ey Pays % «5 Pee 
A, 6, ° . °y 6, ° . . (5.1) 
<0) ets . ° “9 eno e e 


in which /25, 6, and z)=0 indicate the initial conditions. The history ends when the neutron 
is absorbed or reflected from the semi-infinite medium. We have adopted an additional termi- 
nation rule in the present work, ending a history when the energy of the neutron drops below 
0.5 ev. In order to generate a history, one must proceed from collision to collision, computing 
(1 41,Ons1.2n41) from (£,,0,,2,) recursively by random sampling. 

By “random sampling from a distribution /(7)’’ we mean, in principle, the following 
procedure. A set of numbers, w;, #2, . . ., is prepared whose frequency distribution is made 
to approximate f(r) as closely as possible. One of these numbers is then selected at random 
and is called a random variate from the distribution f(#). It is, in fact, sufficient to prepare 
only a single set of random variates, from a uniform distribution f(7)=1(0<a2<1). These 
numbers, distributed uniformly between zero and one, are called random numbers. When 
suitably manipulated they can be made to yield random variates from any distribution, as 
shown below. 

Ina machine calculation requiring many random numbers it would be awkward to prepare 
and store a large set of random numbers. — Instead, it is customary to use so-called pseudorandom 
numbers which are generated systematically but have, nevertheless, sufficient randomness 
for Monte Carlo calculations. We have employed the method of congruential multiplication ® 
in order to produce a sequence of pseudorandom numbers 


Str 
Ser 
Str 
— 
ey | 
bo 
—_ 


These numbers are fractions (between zero and one) which are obtained by normalizing a set 
of integers 1, 


£, —2-5y,. (5.3) 


The integers, in turn, are obtained from the following recursive scheme: 


Io CG ~ 
; eK» | 235 (5.4) 
Mn41—9 4, MOQ 2", 


where ¢ and & are odd integers. The last equation states that ,,, is congruent to 5‘r, modulo 
2%- an other words, that w,.; is the remainder obtained upon division of 5*r, by 2%. In the 
course of our Monte Carlo caleulations we have used three pseudorandom number sequences 


of the type (5.2), with A=11, 13, and 15. This will be discussed more fully in section 5.4. 
See, e.g., the article by J. Todd in reference [13c]. 
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It can be proved (see app. B) that the sequence (5.2) contains 2° different numbers which 
are distributed with uniform density between zero and one. The successive random numbers 
used in the Monte Carlo calculation have been obtained by picking success’re numbers in the 
pseudorandom number sequence (5.2). This is equivalent to a random choice, because the 
pseudorandom number sequence closely resembles a true random number sequence in the sense 
that there is no noticeable correlation between successive numbers. This lack of correlation 
has not been proved, but is an empirical fact that can be demonstrated by statistical tests. 
In appendix B, section 9, we present the results of such tests for 237,000 pseudorandom numbers 
actually used in our calculations. 

One procedure for selecting a random variate «from a distribution f(r) by the manipulation 
of random numbers depends on the possibility of finding a random variate y(2) which is dis- 
tributed uniformly between 0 and 1. We let 


y(r= fi x’ )dxr’. (5.5) 
le., y(w) Is the cumulative probability distribution associated with f(r). Note that 
probability {y<y’ probability {2 >.2(y") is Wer =" (5.6) 


so that ¥ is equally likely to have any value between 0 and 1. We therefore identify y with a 
random number &, and determine a random variate « by inverting the equation 


~] 


Fa’ dr. (; 


It 


A second method, called the “rejection technique” applies to distributions with a finite 


range, say O<r<a. Let L be the maximum of f(z). The sampling procedure consists of the 
following steps. First, two random numbers, say & and &, are selected. If Lé& <f(aé,) we 
chose the number ag, as a random variate from f(r). If Lé, >f(aé,), we reject the random num- 


bers & and & and repeat the testing process with a new pair of random numbers. The 


justification for this selection procedure is supplied by the observation that 
probability (acceptance of «—aé, probability — Lé, > flag) = a dé, I 


The average number of random number pairs that must be tested before finding an acceptable 
random variate Is given by 

I al, (5.9) 
Ls Jo 


flag jdé, 


We shall now describe the actual sampling procedure for generating a neutron history: 


i.e., the rule for going from (12,,6,.2,) to 16, ).8, «1 
a. Location of the Next Collision 


The probability distribution of 2,.; is 


l | 9 
, Sn, /ME, ; Cos 9, ALE,), (5.10) 


where 


Pee ian ead. (5.11) 


is the pathlength traveled by the neutron between the nth and n+ Ist collision. The distribu- 
tion (5.10) is derived from the assumption that the spatial distribution of target nuclei is random 
and that the probability of a collision per unit pathlength is 1/A(/2,).. Applying the cumulative- 
probability method, we select a random number &, and determine s,,; from the equation 


C 


Str 


re , ds ree ; 
« A( En X( E n) 


Sn+l 
The pathlength determines the position of the next collision, which is given by 
in = Ent COSOnSn 41 = Zn — €O88,A(EL,) log Eq. (5.13) 
b. Type of Collision 


The probabilities p)(y==1, 2, 3, 4) for the occurrence of each of the four possible types of 
collisions are given by (3.4), (8.5), (3.6), and (3.7). To make a random selection, we divide 
the region between O and 1 into four intervals proportional to the probability for each process. 
A random number &3 is then picked and the type of collision is determined by the interval into 


which &3 falls: 


O<iscp : clastic scattering from hydrogen 
Pic Sa Spi} pro : elastic scattering from oxygen 
Pit poscés < pit pot ps : inelastic scattering from oxygen 
Pit Pot ps Sa S pit pot pst psr= | : absorption. 


c. Energy Loss 


If the w+ Ist collision results in a scattering, the energy of the neutron is changed from £, 
to /,,, according to the single-seattering law (4.7); i.e., the probability distribution of /,,; 1s 


PF (E23; 2,): pt = . TAB: waa; B,). (5.14) 

where 
(M4+-1)°F,4,-U4+M2K3)E, 
2M,K,E, , 


COS O¢ p41 (5.15) 


The choice of the energy loss can be made indirectly through the choice of the center-of-mass 
deflection 0,.,.;. When the distribution of deflections is isotropic (for collisions of type 1 and 
» 


od), We Sel 
COS 9. n41= 2&—1, (5.16) 


where & is a random number. For clastic scattering from oxygen, the angular distribution is 

anisotropic and is given by the Legendre expansion (3.9). In this case it is convenient to use 

the rejection technique. Provisionally, a deflection angle is chosen according to (5.16). An 
ae ERAN: a : P of 
additional random number, €& 7, is then selected and a test is made to determine if 

; o¢ rd rye) — 

Joi cos” *(2&,—1);#,,} > L(A, )E, , (5.17) 

where L(/2,) is the largest possible value of f,(0,;/2,). If the inequality (5.17) is satisfied, the 

provisional random variate given by (5.16) is accepted; if not, the selection procedure is re- 


peated with two new random numbers.” 


1 Using formula (4,9), noting that the range of the angular distribution (in cos) is a= 2, and using the differential cross sections as shown in 
figure 3, we find that the average number of random number pairs that must be tested to make a selection has the following values: 6.8 at 14 
Mev, 4.3at 9 Mev, 2.2 at 6 Mev, 2.6at 3 Mev, 2.8at 1 Mev, and 1.03at 0.3 Mev. The sampling efficiency at higher energies could be improved by 


more elaborate schemes, but we have not considered this worthwhile. 
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d, Change of Obliquity 


The prescription for selecting the obliquity after the ”-+-Ist collision can be obtained in 
two different ways: by a direct geometrical argument, or by an application of the scattering 
function derived in section 4. Starting with the latter, we note that, according to (4.11) the 
probability distribution of 6,,; (for fixed 6,,,.; and 6,) is given by 


l as 
senitpictads ; 3 = —— (5.18) 
1, COS(a;,—86,,) —COS6,, 4, } 7. COSA, +1;—COS(a;, +8, ) }? 
where 
MA fe, A111] & ~ 16 
COS @» _— pens (5.19) 
——— BE,  2(M,+) VE, 
and 


COS(Q;, +-0,) [€COS6,,4.1 <COS(aj,—9G,,). (5.20) 


Applving the cumulative-probability method, we pick a random number, &, and set 


1 (‘int sinddé 
1—é; | } n 
’ TO On+1 COS(Q;,, 6,,)—cosé >?) COSA— COS(a,, + Din) 3? 
] 1 J COSO,41—COSa;, COSA, " 
=] COs : : . (o:21) 
1 sina,;, sing, 
The inversion of this equation yields the rule for ealculating @,4;: 
COSO, 4.1 =COS@;,COS 8, -+SiNa;, Sind, COS(1é;5). (§.22) 


Let 0,41 be the polar deflection and Xn+1 the 


The geometrical derivation goes as follows. 
Ist collision. 


associated azimuthal deflection in the laboratory system, resulting from the n-4 
ynii iS assumed to be measured with respect to a plane that contains the z-axis and the direction 
The various deflection angles and 


of motion of the neutron before the n-+-Ist collision.) 
The application of the law 


obliquities form a spherical triangle that is sketched in figure 5. 
of cosines to this triangle vields the relation 


to 
w 


COSO, = COSO, +, COSA, +-SiINO, .; SING, COSX, «1. (5. 


NEUTRON DIRECTION 
Zz - DIRECTION AFTER n'th COLLISION 


NEUTRON .DIRECTION 
AFTER n+i'st COLLISION 








POINT OF nel'st 
COLLISION 


hicure 5, Spherical triangle illustrating the kinematics of a neutron scattering. 
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Now, according to (4.5) and (4.10), 0,4; 1s identical with @;,. Furthermore, it can be assumed 
that the azimuthal deflection angle x,;; is distributed randomly between 0 and 27, so that it is 
legitimate to replace cosx,,; In (5.23) by cos, where £6 is a random number. Thus (5.23) 
is in fact identical with (5.22). 


5.3. Albedo Estimates 


The estimate of the number albedo Ay involves merely the determination of the fraction 
NUN, of the number of histories that end in reflection.’ To obtain an expression for the 
statistical error of such an estimate we shall be somewhat formal and assign to each sampled 
history a “score” 22, which has the value 1 in case of reflection and is 0 otherwise. Thus, 2? is 
a random variable which assumes the value 1 with probability Ay and the value 0 with proba- 
bility 1— Ay. 

The mean value of 2? is 


R Ay: Tt (1 —Avy)-0 =Ay (5.24) 


and its variance (defined as 222— 2 ) is 
var R=Ay(1—Ay). (5.25) 


Suppose now that we sample and examine Ny histories with scores P,,/2o, .. . Ryo. An esti- 
mate of the number albedo is given by the average score 


| N 
Az , (R, +- PP, | ha fe ear iy )=y> (5.26 
v=, : vo) =N, ) 
The variance of this estimate 3 is 
var AX : var 2=- _ 1y(1—Ay) (5.27) 
< A“ N r < y, «4 N 4 N . Seka 
N 0 No 


An estimate of the energy albedo can be obtained by multiplying the number albedo 


estimate NON, by the quantity /7/), where / is the average energy of the reflected neutrons. 

. : . m ; : sae. e 
More formally, we assign to each history a score T which is equal to /7//) in case‘of reflection 
and which vanishes otherwise. The random variable 7 has the mean value 


> E 
'=Ay KE? (5.28) 
<0) 
and a variance 
——_ on? ] — —2 ; ; 
var T=17—1 pe vl ARE ). (5.29) 
“0 
An estimate of the energy albedo is provided by the average score 
oe = NE 
Ag=— (7,4+-7,+-... +7y)=> >=: (5.30) 
7 N, 7 fe No Ky 
The variance of this estimate is 
l — —2 
var Af=——— (AvE*— AE ) 
Nok 
l —_ ‘ ‘ _ 
KY? var R-+-Ay, var EF). (5.31) 


=— — ( 
Nok 


12 Crude estimates of the differential albedo can be obtained by counting the neutrons that are reflected with energies and obliquities in speci- 
fied energy-angle intervals. The statistical accuracy of such estimates is given by formulas analogous to those for the number albedo. 

3 Actually one does not know the exact value of .1y, but an estimate of var (14) is obtained through replacing Ay in (5.27) by An, and No 
by No—1 in order to make the estimate unbiased. 
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The two terms contributing to this variance are both non-negative. The first arises from 
fluctuations of the number of reflected neutrons and the second from variations of the energies 
of the reflected neutrons. 

The significance of the variances var Ay and var .1; must be interpreted with reference 
to the probability distributions of Ay and AZ. Both of these quantities are the averages of 
N, identically distributed independent random variables. According to the central limit 
theorem of statistics [14] their distributions approach a Gaussian distribution in the limit of 
large Not One can therefore make the following statements about the deviations of the 


albedo estimates from the true albedo values: 


) »/9 
Probability fl AX—Ay| > hy var 1%} Probability Afs—Ag| >kyvar At \ - | e *''"dz. 
7. 


5.4. Correlated Random Sampling 


To establish the dependence of the albedo on the assumed conditions (source energy and 
obliquity, collision cross sections, ete.), it is better to use a difference method than to make 
two or more separate calculations. The difference method can be designed to eliminate 
irrelevant statistical fluetuations, by keeping constant all random elements in the calculation 
which are not affected by the change of assumed conditions. Thereby greater accuracy can 
be achieved with a given amount of computation. 

We recall from section 5.2. that the following random numbers specify the nature and 
outcome of a collision: € determines the quantity s,., A(/¢,) which is the distance (in mean 
free paths) traveled by the neutron between the previous collision and the collision under 
consideration; &€ determines the type of collision; €, and, for anisotropic elastic scattering 

1) 


from oxygen, additional numbers &,”, £7, ... determines the center-of-mass deflection 0,413 


Ey , by; 
finally, € determines the azimuthal scattering deflection x,.;. The above characteristics, 
together with the initial conditions, specify a history completely. 

Suppose now that we want to make a direct determination of the difference in the albedo 
brought about by a change of the obliquity of the incident radiation. We ean do this by 
generating simultaneously two correlated sets of histories, one set with one initial obliquity, 
the second set with another. The sets should be correlated in the sense that for each history 
in one set there exists a counterpart in the second set in which exactly the same random numbers 
are used for determining the outcomes of corresponding collisions. Thus two matched 
histories are alike, except that one of them is “rotated with respect to the other” and may 
be shorter because of the different initial obliquity. We call this case one of marimal cor- 
relation.” 

In order to find out how a change in the assumed cross sections affects the albedo, we 
may again generate two (or more) correlated sets of histories, using identical random numbers 
for corresponding choices. Then two matched histories will at least be identical in regard to 
the distances (in mean free paths) between successive collisions, and in regard to successive 
azimuthal deflections. However, energy losses in corresponding collisions may be different. 
Consequently all the parameters of a history which depend on the energy may also be different. 
We call this case one of partial correlation. Even in this ease, the use of identical random 
numbers tends to minimize accidental differences. 

It is most economical, but not absolutels necessary, to generate several correlated sets 
of histories simultaneously. In the present work we wanted to compare the albedo under 
many different conditions, and used the neutron histories also for the solution of other prob- 


lems, so that simultaneous computation of all correlated histories was not practicable. We 
‘In fact, the distribution of Ay is the No-fold convolution of a binomial distribution with itself, and is very closely approximated by a Gaussiar 
even for No as low as 10. The numerical results for the energy albedo indicate that the distribution of Age is a reasonable approximation to a 
Gaussian for Ne= 100 
18 This case also arises in the comparison of the reflection of neutrons after 1, 2,... 7 collisio 
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therefore had to devise a technique for doing separate but correlated calculations. This was 
facilitated by the use of pseudorandom number sequences, which according to (5.3) and (5.4), 
are completely determined by a “starting value” ¢ and a “multiplier” 5*, and can be conven- 
iently regenerated provided we remember these two parameters. 

The number of different random numbers required to generate two histories is usually 
not the same, even when they are correlated. In the first place, even in the case of maximal 
correlation, the two histories will tend to differ in length. Secondly, in the case of partial 
correlation there is a statistical variability in the number of random numbers &”, &, 
that may be required to select a scattering deflection by the rejection technique, according 
to (5.17). Thus it could happen that in spite of the use of the same pseudorandom number 
sequence (with the same e and k), different portions of this sequence would be used for 
supposedly correlated histories. This difficulty was circumvented by the use of many 
sequences of pseudorandom numbers instead of just one. A primary sequence was introduced 
with an arbitrary odd starting value ¢ and a multiplier 5". Suecessive members of the primary 
sequence were then used as the starting values of secondary sequences generated with a 
multiplier 5". Successive secondary sequences were used for the calculation of successive 
histories. When the need arose for random sampling by the rejection technique, the requisite 
random numbers &?, &?,... were obtained from a tertiary sequence. The tertiary sequence 
used a multiplier 5“, and a starting value equal to the last member of the secondary sequence 
currently being used.!® 

We now discuss the effect of correlated sampling on Monte Carlo estimates of Ain tee and 
Ap ae, where sly and vty, are albedos under condition 1, and As and Ag the albedos under 
condition 2. Proceeding as in section 5.3, we assign to histories generated under condition 1 
the scores Rand 7, and to the histories generated under condition 2 the scores Rand T. Albedo 
differences can then be estimated as follows: 


. me 1 ~ Pe a ie = 
A*— At “yr (R,— Ry) + (R.— hz) +... +CRy,- xy) }; (5.33) 

. ae _ m Ff mn ci 
A} pe YS \ (T; 7) (T,- i>) (Ty, Ty) S (o 34) 

0 
The variances of these difference estimates are 
rey ( 1* {* I Feces R 9 (R PR — R r¢ PR lL we PR (5 QF) 
Val Lin~ 4 \) \ , Val i—a.p i, wy Val i Vali i Val i 9 0.00) 
oe l m9 er et = a 2 aiais, 
var (AS—A5)= ye var T—2p,(T, T)V¥ var T var T+ var T}- (5.36) 
+ oO 


The quantits 


o( RR) =~ RR—-RAK (5.37) 


Vvar Rvar R 


is the correlation coefficient of R and R. Similarly, 


~ 
ryy 


_ e707 ’ 
ps( 1 1) = (5.38) 


~ 


yvar 7 var T 


Q 


~— 


6 We have set kyl, ko=15, and k 13. 
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is the correlation coefficient of T and 7. For uncorrelated pairs of histories, the correlation 


coefficients vanish. The objective of correlated sampling is to reduce var (A*—A¥*) and var 
(A‘- Ax) by making the correlation coefficients positive and as large as possible. 

The magnitude of the actual correlation coefficients must, in general, be estimated from 
the sampled data, and some typical values will be given in section 7. It is possible, however, 
to put a theoretical upper limit on p(?, I). Using (5.25) and (5.37) we find that 


~ 


~ fe. AA. o 
o(R,R)=— es ; (5.39) 


~ 


¥ Ay(1— Ay) Ay(1— Ay) 


~ 


Without loss of generality we may assume that sly<-ly. 


Therefore 


probability R=1 and R l } <probability R li=Aw, (5.40) 
so that 
RR < Ay. (5.41) 


Accordingly the correlation coefficient must obey the inequalits 


A (1—Ay) 
x = <1. (5.42) 
Ay(1 Ay) 


eR) <¥ 


The greatest possible value of p(/?, 2) corresponds to the case of maximal correlation discussed 


5 ; * + » 
above. It leads to a reduction of var (Aly—«Ay) by a factor 


2 
|— -* — 
Ay) Io (5.43) 
‘iy 1—Ay 


6. Computations and Results 


6.1. Machine Computations 


The Monte Carlo calculations were programed for the IBM-—-704 computer at the National 
Bureau of Standards.'* The neutron histories, generated in order to determine the albedo, 
were concurrently used for the calculation of the neutron flux in infinite media and in slabs,'* 
so that it is difficult to state how much computer memory and computing time were required 
for the albedo problem alone. Altogether, a little more than 5,000 words of magnetic core 
memory were required of which not less than one-third were allocated to the albedo problem. 
For each case treated, 3,000 neutron histories were generated and analyzed. This amount of 
computation gave the desired accuracy and required approximately 45 min of computer time. 
Not more than half of this time was required by the albedo problem. The use of magnetic tapes 
was not necessary, but was an optional device for storing the results of one machine run and 
combining them with the results of a later run. Intermediate print-outs of the results were 
made after groups of 1,000 histories had been processed, in order to follow the course of the 
computations and to have a check on the performance of the computer. 

In line with the exploratory nature of the investigation, the machine program was kept 
flexible. This made possible the exploration of new aspects of the problem as they emerged in 


17 The program was based in part on an earlier code [15] developed at the Westinghouse Atomic Power Division, The authors are indebted to 
Dr. B. H. Mount and Mrs. Heidi Kuehn who very graciously made their code available to us, 
1@ To be reported in a later publication. 
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the course of the work. In the simplest version of the program only elastic scattering from 
hydrogen and oxygen was taken into account, and the print-out of results was limited to tables 
of the differential albedo. In the latest version all relevant collision cross sections were in- 
cluded, and much additional information was printed out that had been found necessary for an 
understanding of the physical and statistical significance of the Monte Carlo results. This 
additional information included the mean and mean square energy of the reflected neutrons, 
the contributions of different orders of scattering to the number albedo and energy albedo, 
listings of the particular histories in which reflection occurred and of the energies of the reflected 
neutrons, and a statistical analysis of the pseudorandom numbers used. Intermediate versions 
of the program vielded only some of this information. 

The analytical formulas for the single-scattering albedo, developed in section 4, were used 
mainly for the interpolation of results, and single-scattering computations were done with the 
Monte Carlo code. This was simpler to do once the machine program was available but slightly 
less accurate. 

Altogether 38 different cases were treated which are listed in table 6.1. This table contains 
the source energy and obliquity as well as a brief indication of the assumptions about the cross 
sections. The computations fall into two groups: 30 cases were done with one set of pseudo- 
random numbers (analyzed in app. B), and the remaining 8 cases with another set. A high 
degree of correlation was introduced by the repeated use of the same set of pseudorandom num- 
bers in order to increase the accuracy of intercomparisons. We believed this to be a sound 
procedure, because the set was large enough (consisting of approximately 237,000 numbers) 
to be statistically representative. 

TaBLeE 6.1. Summary of calculations 


The entry in the left column serves as index to the differential-albedo tables in appendix C. Each result for the number albedo Aw and the 
energy albedo Ag is based on the analysis of 3,000 neutron histories, with 64» and 6.1¢ indicating the statistical standard deviations. All runs 
were done with the same set of pseudorandom numbers (described in app. B), except the runs indicated by an asterisk (*) which were all done with 
another set of pseudorandom numbers. 





Table Eo( Mev) 6. at 6A, A, 6A, Assumptions 
1 0.3 0 0. 296 0. OOS 0. 048 . ang) 7 
None 1 0 . 370 009 (. 135) (*). 
2a 1 0 . 361 009 (. 115) (*); elastic scattering from oxygen assumed isotropic. 
2b 1 0 381 009 . 148 0. 005 
2c 1 0 213 . 006 . 091 Single scattering. 
3a 1 15 . 450 009 159 : 
3b ] 15 114 006 . O85 _.| Single scattering. 
43 1 60 518 009 . 183 
ib ! 60 . 142 . 006 . 099 Single scattering. 
5a 1 75 619 . 009 . 241 
5b l 75 213 O07 . 144 . x Single scattering. 
ba l 90 767 . 00S .415 ; ‘ ' 
ib 1 Hn) 493 009 . 347 = Single scattering. 
7a 3 0 061 004 . 023 ‘ 2-cm slab; see section 7.9; (*). 
Tb 3 0 .131 O06 . 038 5-em slab; see section 7.9; (*). 
7c 3 0 . 163 007 . 044 10-cm slab; see section 7.9; (*). 
None 3 ° . 168 . 007 (. 044) pha ig 2 
Sa 3 0 178 . 007 (. 049) (*); elastic scattering from oxygen assumed isotropic. 
Sb 3 0 . 167 . 007 . 041 . 003 
SC 3 0 28 . 003 A | ..-| Single scattering. 
Ya 3 45° | . 2a2 OOS . 052 ers | 
9b 3 45 035 . 003 . 023 ae Single scattering. 
10a 3 60 solo 008 (. 080) }-.--- 
10b 2 60 066 005 (. 039) | -- Single scattering. 
lla 3 75 472 . 009 . 144 beomaccaee 
Ilb 3 75 . 143 . 006 . O84 _| Single scattering. 
l2a 3 v0 . 728 . 008 (, 339) < 
12b 3 0) . 480 009 (. 290) _.| Single scattering. 
1: 3 (t) 371 009 (.110)|_- 
l4a 6 0° . 157 007 . 052 . 003 
l4b 6 0 171 007 . 059 .004 | Absorption disregarded. 
l5a 9 0 107 . 006 . 027 .002.| Excitation of 6.1-Mev level in oxygen by inelastic scatter- 
ing. 
15b 9 0 . 138 006 . 041 .003 | Absorption disregarded. 
lie 9 0 .131 | . 006 . O47 __. Absorption disregarded. Inelastic scattering from oxygen 
treated as elastic. 
Ifa 14 0 . 124 . 006 . 040 .003 | Excitation of 6.1-Mev level in oxygen by inelastic scatter- 
ing 
16b 14 0° . 163 007 . 038 .002 Elastic scattering from oxygen assumed isotropic; excita- 
tion of 6.1-Mev level in oxygen by inelastic scattering. 
l6e 14 0 . 165 . 007 .079 .004 | Absorption disregarded. Inelastic scattering from oxygen 
treated as elastic. 
16d 14 0 O64 004 . 023 .002 Inelastic scattering treated as absorption. 


(t) Isotropie source, 
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6.2. Differential Albedo 


Tables 1 to 16d of appendix C contain the Monte Carlo results for the differential albedo 
in the form in which they were printed out by the computer. The number of reflected neutrons 
(for N,=3,000 neutrons) is given in a two-way classification according to the spectral energy, 
FE, and the cosine of the obliquity of emergence, cos 6. By normalizing the entries in the tables, 
i.e., by dividing through by No, one obtains an estimate of the differential albedo averaged over 
a small energy-angle interval. Unnormalized numbers are given in the tables in order to pre- 
serve their statistical significance, the standard deviation of each entry being proportional to its 
square root. The energy classification is divided into 10 intervals, labeled 7=1, 2,...10, such 
that the 7th interval pertains to neutrons with energies between (11—7)/104) and (10—7)/104). 
The last interval extends from /)/10 to 0.5 ev rather than to zero energy. The obliquity 
classification is also divided into 10 intervals, labeled #==1, 2,..., 10, such that the Ath 
interval pertains to neutrons with obliquity cosines between (1—/)/10 and —//10. Two sets 
of marginal totals are also given which represent the angular distribution integrated over all 
energies, and the spectral distribution integrated overall directions, respectively. 


6.3. Number Albedo and Energy Albedo 


Table 6.1 contains all the results obtained for the number albedo Ay and the energy albedo 
Ay.’ The number albedo was obtained by adding the marginal totals in the corresponding 
differential-albedo tables and dividing the result by Nj=3,000. The energy albedo was obtained 
as the product of the number albedo and the average energy /¢ of the reflected neutrons. In 
most of the calculations «4, was computed directly in the course of the machine calculations. 
In some of the earlier cases, Z had to be estimated from the differential-albedo tables.2?. Nu- 
merical experimentation indicated that such estimates could be made with an accuracy of 
approximately 1 percent. Values of the energy albedo based on such estimates of are indicated 
in table 6.1. by being placed in parentheses. 

The standard deviations 61, yvar Ay and 6:1, qvar ely, are also shown in table 6.1- 
They were computed according to (5.27) and (5.31), respectively. The determination of 6.1, 
required knowledge of / as well as of 4%. Accordingly, 614, is given only for those cases in 
which £? was obtained as direct output of the machine program. The relative standard devia- 
tion é1y uy, obtained with 3,000 histories, is of the order of 0.05. The corresponding relative 
standard deviation é:4,/A,, in the cases examined and presumably also in the other cases, is 
roughly 135 times greater than 6:1, «Ly. 


7. Analysis of the Results 


7.1. Single-Scattering Contribution to the Albedo 

ry. -_ ° . . 1) 1 : ° 
fable 7.1 contains asummary of theratios 1; /Ayand a1, /A, for various cases. In general, 
single scattering contributes about 20 percent of the number albedo and more than 50 percent 
of the energy albedo. The contributions are particularly large at 1 Mev, because at this energy 
the ratio of the scattering cross section of oxygen to that of hydrogen is large, and the differential 


oxygen scattering cross section has a peak in the backward direction, 


vIn sec. 5.3. the albedo estimates were indicated by A*y and A*g in order to distinguish them from the true values 4y and w4¢. In the 
following, we shall ignore this distinction and simply use Ay and 21g to Jabel the Monte Carlo results. 


oF ne 
20 The most accurate estimation procedure was found to be one based on the “cumulative albe "| ar'| sindd6.A(E’ 0; 0,85), which ean 


be obtained from the differential-albedo tables for F.=i/10Fo,i=1, 2. . 10, and by interpolation for other values of &. It was also necessary to 


take into account the characteristics of the spectrum of the single-scattering albedo, in order to obtain good accuracy. 
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The estimation of the albedo ratios is a problem with maximal correlation as detined in 
section 5.4. We have estimated the standard deviations of these ratios in a few typical cases 
and find, for example, that AY? /Ay=0.3040.03 and A,?/A,=0.61--0.02 for Ay=1 Mev 
and @=0°. 

The ratio Aj’ /Apgisarelatively slowly varving function of the source energy, whereas Aj 
and sl, separately vary rapidly. This suggests that the energy albedo for arbitrary energies 
could be obtained with fair accuracy by combining a single-seattering calculation with an 
interpolated value of Aj /Ap. 


Taste 7.1. Patio of single-scattering albedo to total albedo. (The results for 6, 9, and 14. Mev correspond to cases 
I4a, 15a, and 16a, respectively, as listed in table 6.1). 


ly Ag 
Ly 1K 
i al 
if) 15 60 75 90) 0 15) 60 to 90 
FE BE 
Mey (Mev 
0.3 0.11 0.3 0. 54 
] 30) 0.25 0.27 0.34 0. 64 1 Hl 0.54 0. 54 0. 60 0.84 
3 17 «15 1 30 ‘6 4 55 45 4s _ 58 gh 
i 22 t a 
y 17 y $2 
14 aT 14 AO 





7.2. Contribution of Higher Orders of Scattering 


Figure 6 shows the contributions of successive orders of scattering to the albedo for a few 
typical cases. As many as twenty orders are shown to contribute to the number albedo,” 
Whereas only about five orders contribute significantly to the energy albedo. The relative 
importance of the first few orders of scattering increases with the source obliquity and energy. 
In general, an analytical orders-of-scattering approach does not seem to provide a practical 
alternative to a Monte Carlo caleulation. The contribution of the second order of scattering 
is not sufficient fora precise albedo determination, and calculations for the higher orders would 


be prohibitively complex. 
7.3. Energy Spectra of Reflected Neutrons 


Figure 7 contains histograms of the energy distribution of reflected neutrons for incident 
The contributions of single scattering and multiple 


I-Mev beams with various obliquities. 
The spectra are characterized by two rather sharp 


scattering to the spectra are indicated. 
peaks: one, primarily due to single scattering, occurs at the upper end of the spectra; the other, 
due to multiple scattering, at the low-energy end. The single-scattering contribution can in 
turn be resolved into two components, due to scattering from oxygen and hydrogen. For 
perpendicular incidence, only single-seattering reflection from oxygen occurs. As the source 
obliquity is increased, hydrogen scattering becomes possible, and the oxygen contribution is 
At 45° and 60° the oxygen and hydrogen components are 


shifted toward higher energies. 
In the limiting case of grazing incidence 


clearly separated, whereas at 75° and 90° they merge. 
there isa hydrogen component distributed uniformly over all energies, and an oxygen component 
Further details about 


The 


energy spectra for other cases are similar to those shown in figure 7 and can be obtained from 


Which gives rise to a peak at an energy slightly below the source energy. 
the oxygen single-scattering contribution can be obtained from appendix B (fig. 15). 


the tables in appendix C, 


Some contributions from orders beyond the twentieth would have occurred if our calculations had not been cut off at an energy of 0.5 ev. 
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Figure 6. Contributions of successive orders of scattering to the aibedo, 
7.4. Angular Distribution of Reflected Energy 


Figure 8 shows plots, in histogram form, of the quantity 


- p 

A,(0; Eo, 90) - | "ACE, 6; Hy, 6) 524, (7.1) 

J0 4) 

for £,=3 Mey and various values of 6. The angular distribution is dominated by the single- 
scattering contribution. In particular, it should be noted that for large source obliquities the 
distribution does not follow the cosine-law which holds for perpendicular incidence and has 
often been assumed to be universally valid. Rather the emerging radiation tends to peak for 
cosé=0; 1.e., in directions parallel to the boundary of the semi-infinite medium. This is due to 





__# The histograms are a somewhat crude estimate of the angular distribution, because each angular element was obtained by multiplying the 
—— ial element of the corresponding table in app. C by the midpoint of the appropriate energy interval and then summing over-all energy 
ntervals, ; 
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FiGuRE 7. Energy spectra of reflected neutrons. 
The shaded areas indicate the contribution of single-scattering. Source energy 1 Mev. 
For 4,=3 Mev and @=90°, we find from 


“0 ’ 


the behavior of the single-scattering component. 
eq (A.21), (A.22), and (A.24) of appendix A, and from the cross sections given in figure 2, that 


AY ~0.049(1.00+-1.05 sin 2 6+-7.14 sin * 6), (7.2) 


which is in agreement with the Monte Carlo calculations. 


7.5. Dependence of Albedo on Source Energy and Obliquity 
The dependence on the source energy is irregular. It is correlated with the complicated 
energy-dependence of the attenuation coefficient for water, and is further influenced by the 
relative probabilities for scattering by hydrogen and oxygen, and by the energy-dependent 
anisotropy of elastic scattering from oxvgen. Just like the attenuation coefficient, the number 
albedo and energy albedo have a peak at 1 Mev where they are 2 to 3 times larger than at the 
other energies considered. The correlation between attenuation coefficients and albedo is 
complicated, however, by the interaction of the various factors, so that interpolation with 
respect to source energy is not easy, and is perhaps done best by the single-seattering analysis 
indicated in section 7.1.7 
"Ee tliee alsin of an albedo experiment with a broad source spectrum (such as the Po-Be source used by Strickler, Gilbert, and Auxier 
[16]), it might be preferable to make a Monte Carlo calculation with a realistic source spectrum instead of combining the results of calculations for 


monoenergetie sources, 
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Figure 8 Angular distribution of reflected energy. 
Source energy 3 Mev. A single-scattering curve for 40=90° is shown which corresponds to eq (7.2). 


The dependence on the source obliquits is smooth and interpolation should be easy. By 
way of example, figure 9 contains curves of Ay and 1, as functions of 4, for source energies of | 
and 3 Mey. All four curves shown have similar shapes, starting out with a rather flat slope 
near 4=0° and rising verv steeply near @=—90°. The closer @ is to grazing incidence, the 
smaller are the differences of the albedos for different source energies. 


7.6. Effect of Absorption and Inelastic Scattering 


If absorption were disregarded, the number albedo at 6 Mey would be increased by 9 per- 
cent and the energy albedo by 13 percent. At 9 Mev, the corresponding increases would be 29 
percent and 52 percent, respectively. The effect of absorption is greater on the energy albedo 
than on the number albedo because the relative probability of absorption Is an increasing fune- 
tion of the energy (up to 14 Mev). The effect of absorption on the energy spectrum is shown in 
figure 10. Although the effect is large, it introduces no large uncertainty into the ealculations 
because the cross section for 1—a@ and n—p processes are reasonably well-known, 

The effect of inelastic scattering is more difficult to determine because so little is known 
about the cross section for this process. We have made computations under three different 
assumptions: (1) Absorption taken into account, inelastic scattering assumed to excite oxvgen 
to a level 6.1 Mev above the ground state; (2) absorption taken into account, inelastic scattering 
also treated as absorption; (3) absorption disregarded, inelastic scattering treated as if it were 
elastic, 

As has been discussed in section 3.2, assumption (1) is presumably the most. realistic. 
Assumption (3) is the least realistic and has been included mainly because it has often been made 
in previous neutron transport calculations. Assumptions (1) and (2) are limiting assumptions 
in the sense that they allow for the minimum and maximum possible energy loss that a neutron 
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Picture 10, Reduction of energy spectrum due to absorption, indicated by shaded areas. 


Perpendicular incidence, 
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could suffer in an inelastic collision. The albedo calculated on assumption (3) is almost twice 
as great as that calculated on the realistic assumption (1), whereas the albedo calculated on 
assumption (2) is only about half as large. These extreme variations point up the need for 
accurate knowledge of the inelastic scattering cross section. 

The energy spectra obtained with a 14-Mev source under the three assumptions are shown 
in figure 11. The spectrum in case (1) shows three peaks, due to elastic single scattering, in- 
elastic single scattering, and multiple scattering. The inelastic scattering peak disappears in 
case (2) and merges with the peak due to elastic single scattering in case (3). More accurate 
treatment of inelastic scattering would probably give a spectrum similar to that of case (1), but 
with the middle peak reduced in size and “smeared out” over the lower part of the spectrum. 
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Ficure 11. Effect of inelastic scattering on the energy spectrum, 


Source energy 14 Mey. Histograms (i), (ii), and (iii) correspond to the different assumptions in sec. 7.6, and were obtained from tables 16a 
16d, and 16c of app. C, respectively. Perpendicular incidence, 


7.7. Effect of the Anisotropy of Elastic Scattering 


The differential elastic scattering cross section is not well-known for most elements. — It is 
therefore of interest to know how much the albedo would be changed by a change of the assumed 
distribution of elastic scattering deflections. To explore this situation, we have compared the 
albedo calculated with the best available cross sections for oxygen (reference [9], see also fig. 3) 
with the albedo calculated on the assumption that the elastic scattering deflections from oxygen 
are isotropic in the center-of-mass svstem. The difference between the two assumptions is 
greater than any expected realistic modification of the oxygen cross sections. Comparisons 
were made for source energies at 1,3, and 14 Mev. At 1 Mev, Ay is decreased by 2.5 percent 
and Ag by 11 percent. Finally, at 14 Mev, Ay is increased by 31 percent, whereas Ay is de- 
creased by 5 percent. In figure 12 the spectra obtained under the two different assumptions 
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FictrRE 12. Effect of anisotropy of elastic scattering from oxygen on the energy spectrum. 


Perpendicular incidence, 


are compared, There is relatively little change, the main effect being a shift of the single- 
scattering peak toward higher energies in the case of isotropic scattering. 

The albedo differences quoted above were obtained by correlated sampling (case of partial 
correlation). A listing of all significant statistical parameters for the comparison at 14 Mev is 
given in table 7.2 from which it appears that by means of correlated sampling an accuracy was 
achieved that would have required a sample size 3 to 4 times as large with ordinary sampling. 


TapueE 7.2. Estimate of albedo difference by correlated sampling 


Source energy 14 Mev, perpendicular incidence. Elastic scattering from oxygen assumed to be anisotropic in case 16a, isotropic in case 16b. 


Number albedo: Energy albedo: 
Aw, case 16a, Ar, case lfa, 
Aw, case 16b, Ar, case 16b, 


Ar—A g=—0.0012+0.0016 (with correlated sampling), 
+(0.0033 (without correlated sampling), 


Correlated sampling increases effective sample size by factor 3.9, 


Ay—Ay=0.039+-0.005 (with correlated sampling), 

+0.009 (without correlated sampling), 
Correlated sampling increases effective sample size by factor 2.9, 
Correlation coefficient p=0.664, Correlation coefficient pz=0.771. 


Maximal correlation coefficient pmax =0.853, 
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7.8. Low-Energy End of the Spectrum 


Details of the low-energy part of the spectrum were obtained as a byproduct of the detailed 
information obtained for the statistical analysis of the correlated calculations at 14 Mev. 
Figure 13 shows a plot of the spectrum on a logarithmic scale which emphasizes the low energies. 
The spectral shape at low energies is quite similar to that obtained by Spinney [5] who calculated 
by age theory the albedo for 2-Mev neutrons incident on concrete. This suggests that com- 
putations of the spectrum down to very low energies could advantageously be done by a 
combination of Monte Carlo and age theory, with Monte Carlo covering the energy region 
not very far below the source energy in which age theory does not apply. 
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Figure 13. Energy spectrum of reflected neutrons on semilog scale. 
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Source energy 14 Mev, perpendicular incidence, 


7.9. Slab Albedo 


The deeper a neutron penetrates into the medium, the smaller is its chance of returning. 
Therefore, the collisions leading to reflection take place rather close to the surface of the 
medium, and thin slabs have nearly the same reflecting power as a semi-infinite medium, 
This is shown by the results in figure 14 pertaining to the reflection of 3-Mev neutrons incident 
on slabs with thicknesses of 2, 5, 10, and 20 em. (These thicknesses are to be compared with 
the transport mean-free path in water at 3 Mev which is 5.1 em.) The dependence of the 
number albedo on the slab thickness is similar to that found by Foderaro and Obenshain [6] 
at 0.89 Mev. The energy albedo reaches a limiting value at an even smaller slab thickness 
than the number albedo, because the energy degradation in successive collisions makes it 
necessary for a neutron to be reflected very early in its history if it is to carry back an appre- 
ciable amount of energy. 
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Figure 14. Dependence of slab albedo on slab thickness. 


Source energy 3 Mev, perpendicular incidence. a, Number and energy albedo and b, Energy spectrum. 
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8. Appendix A. Integrals Over the Differential Single-Scattering Albedo 
8.1. Energy Spectrum 


The evaluation of (4.15) leads to the following expression for the energy spectrum of neu- 
trons reflected from a semi-infinite medium after exactly one collision: 


3 
A (B3Eo,0) = >> pj (Lo) P(E; 4) { U(a;,b;) —kV(a;,6;)}, (A.1) 
j=l 
where 
aj;=COs (a,+- 4%) 
b;=cos (a,—) (A.2) 


k= (Xo/A) COS Oo, 


and where U(a;,b;) and V(a,,b;) are defined as follows: 


(1) If B k M? K?—cos? @)+sin _ <E<E, 


M,+1 
U'(a;,b;) = V(a;,b,;) =0. (A.3) 
ay [VIB R3—cos? %—sin 6 Fo yo pp PV MGK3—cos? +sin 4% | 
(2) ez k At |seseol i} ) 
M,+1 M,+1 
; . 1 1. _, by+a;, , 
UL" (a5,b5) 57 Sin ho; (A.4) 








. y \ I ; | 
ifk<b,, J (nb) pda (b,;—a,)k—log [(b;+-a,)k—2a,b, a 
—2y a,b ;(b;—k) (a,—h)]}; 





77 , he [=4,, (A.5) 
if k=b,, Vayb)=—5 ay 5 4 
sie me I a (6 ,+-a,;)k—2a,b; 
fk>b,, Va,,b)=—= 57—— Su he 1+ -! ‘|}. 
am, >; Is 1 Pea 2 7 (b,—a,)k J 
. -y PAGK AVP co ope pp [ VVGK3—cos? —sin af 
(3) If Ky “M,4 1 ]sestl = M+1 ——— 
(A.6) 
U'(a;,b;)=1 
V(a,b,) (A.7) 
al (k—b,) (k—a,) ro 
nae | eee Hl 
(4) If K<k ra | 
U(a;,b;)=V(a;,b;) =0. (A.8) 


In certain limiting cases, the formula for the energy spectrum becomes much simpler than 
the above expressions. For perpendicular incidence, the energy must lie in the range 


M,k ,—1 7 M}K3—1 
%| =z | SE <b AF A.9 
a ie es |s#se qe Sa 
and the spectrum is 
—— z a l 
AM (FE; Eo,0) =>) p;(Ey) F; (E; E)) oo (A.10) 
j=1 (0 
bcd : 
A cos a; 
For grazing incidence (6)—>7/2), the energy range is 
MK ,—1f ,[M,K,417 
401 aa C*| Sw E< Eo | = A. 
Ds M41 Sk<si M41 ] (A.11) 
and the spectrum is 
ae 2 oe eee ; 
A® (E; By, =) Si s( En) 5B, (Es Ba) (A.12) 


Figure 15 shows energy spectra of neutrons reflected after a single scattering from a semi- 
infinite medium of pure oxygen (4,=1 Mev, @%=0°, 45°, and 90°). These spectra were 
obtained by numerical evaluation of the factor 


Fy ( E; Fy) € "(a, b2) —kV(a2,b2) } 


in (A.1). They are characterized by a peak at the low-energy limit which is rather independent 
of the source obliquity @. In addition there is another peak at the high-energy limit in the 
case of grazing incidence. 

For perpendicular and grazing incidence the spectrum for reflection from water can be 
derived very simply from the corresponding spectrum for oxygen. For @)=0°, reflection from 
hydrogen is impossible after one scattering, and it is sufficient to renormalize the oxgyen 
spectrum in figure 15, multiplying it by p.(/))=0.485. For 6)=90°, one must not only normal- 
ize the oxygen contribution but also add a hydrogen contribution which is energy-independent 
and has the value p,(/))4F,(E;4,) =0.515 (1/2) =0.2575. 
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Figure 15. Energy spectrum of neutrons reflected after one collision. 


Source energy 1 Mev, oxygen medium. 


8.2. Angular Distribution 


The angular distribution of neutrons reflected after one collision is even more complicated 
than their energy spectrum. We shall confine ourselves to the discussion of some limiting 
cases. The characteristics of the angular distribution are expressible in terms of the integrals 


H,(E.0;E,0) C086 | 
hs(0;Eo,8o) = =|: es aves} dE (A.13) 


which are needed for the evaluation of (4.16). The angular distribution of reflected energy, 
Aj, (0; ),0)), depends on the integrals 


me, es HT, (E,6, Fy; 9) cosé ‘ 
hes 0;E,0x)—= |“ cos 6 — ( o/d) oe _ (A.14) 


The integrations in (A.13) and (A.14) are facilitated by changing the variable of integration 
from £2 to cos a;, noting that 





E-E cos aj+- cos ay+ a,+M? 3kG—1 | (A.15) 
gees: M,+1 
and that 
dE 2K, ate | ./Nf?2? 2. cos* 
jen, OEP {2 cos a;+ M3 K3— sin? a va ro =} (A.16) 


We now take up limiting cases, and drop the index 7 for simplicity. 


Case (1): Heavy nucleus, > 1 
Perpendicular incidence, 6)=0° 
Elastic scattering, A= 
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Note that 
ExEy, X¥=Xo, . 


dk 2Ko 


dcosa M’ ‘a (A.17) 





ie M., ' 
te ~a, F(E-:E)) ~OF. AOE I) 


Substituting expression (4.13) for /7(/2,0; /5,0), we find that 


cos Of(a; fs) 


5) (A.18) 
cos 6— 1 


h(6;£y,0 = 


and that 
hy(0:1),0) = Eyh(0;4,0) (A.19) 


Cast (2): Heavy nucleus, M>1 
Grazing incidence, 6)=90° 
Elastic scattering, A=1. 


. : re ‘ - ’ < : . 
The approximations (A.17) may again be used. For iI E.9;E.,= ) we substitute expression 
(4.14), the result being 


- ®\_1 ("9 fla; Ey) deos | 
h( 6:4, 5 ) =: Fla; Ko) d cosa, (A.20) 


= T —sin Oy sin’ @ — COS” a 


Assuming for f(a;/)) a Legendre expansion of the type (3.9), with expansion coefficients 
fl (ky), we can integrate (A.20) and find that 


l (2/ )! 


: wHi+1. 4. 
h( 8; =D g Sal Bo) (— I)" gee p ryt 2 (00s 4)» (A.21) 
and that 
h ( 0; Fy, ) = Kh 0;Eo, = ) (A.22) 


Case (3): Hydrogen nucleus, \/= 1 
Grazing incidence, 9,— 90°. 


We have /(0@,;/)) =1/2 and 0<a<a2. Using (A.14), (A.15), and (A.16), we find that 
) *sin a . ) 
, 7 Z cosad cosa 2 5 
h( 6;Eo,5 )= — ==" sin 8, (A.23) 
a TJ 0 ysin*é—cos-a 7 
and that 
ey *sin @ ; : ’ 
Tua 2Eo COS ad COs a 1h, A 
he( 0;E0,5 ) : | —=—— gin? 6. (A.24) 
- 7 JO \ sine 8— COs" a oT 


9. Appendix B. Pseudorandom Numbers Generated by Congruential 
Multiplication 


9.1. Proof of Their Uniform Distribution 


We consider the pseudorandom number sequence {&,—2>%r,| (n—0.1, 233__] ) 
| iS j 


derived from a sequence of integers (r,} which is generated according to (5.4) or by the 
equivalent rule 


kn 


7. =¢5"" mon 2”. (B.1 
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where ¢ and & are odd integers. The sequence (é,{ consists of 2° different numbers which are 
distributed with uniform density over the interval (0,1). This can be proved with the use of 
standard number-theoretical results [17]. However, Dr. Morris Newman*! has given a simple 
proof proceeding from first principles which we reproduce here in a somewhat amplified form. 

The uniformity of the distribution of the pseudorandom numbers obtained by congruential 
multiplication is a consequence of the following theorem. 


Tureorem: The se (uence of inte gers 
{ 33 | y \ 
Wns Ly Ba, - « sy wy ly (B.2) 
generated according to (B.1), is a permutation of the sequence 


Lek... 3 (B.3) 
or of the sequence 


a, 7, ERs. 2, BT} (B.4) 


depending on whether the starting number #,=e has remainder 1 or 3 when divided by 4. 
The proof of this theorem is based on the following lemma: 


ae anal ; ie a be ; — 7 P ‘ 
The integer 5°" —1 is divisible by 2™°* but by no higher power of 2, if \k is an odd integer. 


To see this we consider the factorization 


Hk2M | — (Hh2M—1 4) (5 M—1__ 1) — (5k2M-1 4-1) (5M? 4-1) . 2. (5*+1)(5*—1). (B.5) 


The factor 


5*§—1=(1+4)'—1 weer (“)¥ (B.6) 


is divisible by 4 but by no higher power of 2... The other factors in (B.5) have the form 


-.om kom 9] a ho” 
5" 4-1 —(14+4)""+1—24 > E. 


)4,0<m<M-1 (B.7) 
pd Ne 


and are each divisible by 2 but by no higher power of 2... By combining all the factors we obtain 
the lemma. 
We now consider the sequence (e5'"}, n=O0,1, .. . 2%—1. The difference between any 


ek(n-+n/ 


two terms in this sequence, say cd and ¢5*", can be expressed as 


A(n,n’) =05**(5*"”— 1). (B.8) 
We may set 
NH ‘ 2°). (B.9) 


Where a is an integer <2%) and where 6 is an odd integer. By application of the lemma it 
follows that 5’ —1, and therefore A(i.n’), is divisible only by powers of 2 smaller than 2%. 
This implies that the integers w, in the sequence (B.2) are all different from each other. 

It can be seen from (B.8) that the differences between the .,’s are all multiples of 4. Thus 
two different sequences (r,} can be formed, one of which is a permutation of (B.3) and the other 
of (B.4), depending on the starting number ¢. Each of the two sequences contains 2% num- 
bers, and there is no overlap between them, so that together they exhaust the 2* odd integers 

2%) This completes the proof of the theorem. Finally, we note that by computing the 
r,s modulo 2° rather than modulo 2*°, and by using a starting number ¢, that is not odd but 
contains a factor 27, one would obtain a sequence of 2%~*~* different numbers. 


‘Informal memorandum, Applied Mathematics Laboratory, NBS, Oct. 27, 1951. 
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9.2. Statistical Indications of Their Randomness 


In order to be acceptable as a source of random numbers for a Monte Carlo calculation, 
the pseudorandom number sequence {g,; must fulfill the following requirements: 

(1) In any subsequence of reasonable length, such as would be used in generating neutron 
histories, the &,’s must be distributed with uniform density in the interval (0,1). 

(2) Suecessive &,’s must be uncorrelated. 

These two conditions cannot be satisfied rigorously, but empirical statistical tests indicate 
that they are satisfied in very good approximation. For example, elaborate tests on 16,384 
pseudorandom numbers have been reported in the memorandum mentioned in footnote 24. 

We have made some rough-and-ready statistical tests on 237,000 pseudorandom numbers 
that were used to generate a set of 3,000 neutron histories. These numbers were divided 
into three sets, corresponding to the three types of pseudorandom number sequences (primary, 
secondary, and tertiary) described in section 5.4. Set P includes the numbers in the primary 
sequence, set S the numbers in all the secondary sequences, and set 7’ the numbers in all the 
tertiary sequences. Separate tests were applied to each set. In order to gauge the uni- 
formity of the distribution of the numbers, we computed the averages £, and &?. In order to 
check on the lack of correlation between successive pseudorandom numbers we computed the 
mean square difference of successive numbers, (£,.;—€,)°. Asa further test of possible correla- 
tion, the algebraic sign of the differences between successive pseudorandom numbers was 


S 


examined, and an index ¢ was computed which was defined as the excess of positive over nega- 
tive signs. Table B.1 compares the actual values of all these quantities with the values that 
one would expect if the pseudorandom numbers were truly random. We infer from the table 
that the pseudorandom numbers are sufficiently random for our purpose. 


TABLE B.1.) Analysis of the randomness of the pseudorandom numbers 
In order to indicate the significance of the deviations of the empirical from the expected values of the various statistics, the standard devia- 
tions of the empirical values are given. The pseudorandom numbers were obtained by starting a primary sequence with the initial value 
¢=16,743,785,845. The number of pseudorandom numbers in each set is indicated within brackets. 


Empirical value 


Statistic Expected 
value 
Set P [3,000] Set S [220,776] Set T [13,384] 
ee ly 0. 49402449. 00527 0. 49955+0. 00061 0. 50233+0. 00250 
3 1 320374 . 00044 333594 . 00064 335804 . OOZSS 
(En+ En 1 147434 00360 165794 00042 . 1H4614 O17] 
t 0 as +55 H3S +470 34 +116 


10. Appendix C. Tables of Differential Albedo 
The following tables are indexed by table 6.1 and described in section 6.2 of the text. 


TABLE 1. 








Cos @ 
1 2 3 t 5 " 7 S 9 10 sum 
i 
] 0 
2 3 7 4 s 12 14 14 6 OS 
3 1 2 1 3 l ] 3 15 18 0) 
4 2 ] 4 1 2 3 3 14 
E 3 l 3 10 
f 2 5 ' 2 5 7 25 
7 1 } 5 2 3 3 IS 
s l 2 I ; 2 2 2 13 
9 ] l 2 } t 6 1 in 35 
oe cia ia 16 9 4s : 91 101 124 142 656 
sum 19 24 4 72 112 125 127 157 Is4 SAU 
25 This set of numbers was used repeatedly for most of our caleulations, 
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TABLE 2a. 
COS @ 
6 i) 10 sum 
| I 1 
) 5 10 16 20) 32 37 18 30 3 J 201 
$ 7 } 7 1 14 40) 4] 117 
i | ! 2 2 2 4 3 3 3 5 26 
5 2 1 l 5 3 5 4 2 5 28 
6 ! 3 2 4 5 2 11 5 33 
7 2 l 5 +) S 3 9 6 43 
s ] 2 l ] 5 5 (! 3 10 35 
y ] ] 3 5 4 6 15 1] 10 56 
10 3 10 17 13 52 is 52 94 102 123 544 
Sum 10 26 11 SI 108 121 136 172 184 205 1, O84 
TABLE 2b. 
cos @ 
6 7) 10 Sum 
1 0 
2 ! 5 10 21 32 1) ay) 1 172 
3 ] 2 2 5 6 5 12 34 SH 113 266 
} 2 5 3 3 6 8 12 39 
5 l 3 i] 5 11 3 é 35 
6 1 2 1 i 4 8 7 7 34 
7 1 3 7 5 1 4 7 4 32 
s 3 S 5 10 2 5 sS 11 52 
v 7 5 1 4 5 11 5 9y 10 55 
10 4 7 22 37 44 5l 51 4 86 8S 459 
Sum s 17 34 69 101 120 135 193 215 252 1, 144 
TABLE 2c. 
cos @ 
5 6 9 10 Sum 
l . 0 
= l 4 9 19 31 42 56 5 162 
3 21 68 88 177 
4 = 0 
5 0 
t) =" 0 
7 ae 0 
% 0 
qy a 0 
10 0 
sum l 0 4 y 19 31 42 77 68 SS 339 
TABLE 3a, 
cosdé 
5 6 9 10 Sum 
! 1 2 1 3 2 9 
2 4 4 13 +] 16 1 17 21 18 37 150 
| 1 11 14 18 32 21 33 49 47 45 271 
1 1 1 2 7 7 6 7 10 6 8 55 
5 1 4 1 a 4 10 +!) 2 4 t 42 
H l 2 2 2 7 S Ff 10 7 46) 
7 l 10 2 i 5 2 5 3 Y 4 44 
s 1 4 9 6 12 10 7 6 14 y 78 
9 2 5 9 7 9 10 12 10 19 86 
10 S 11 24 42 56 62 5Y S4 100 123 569 
Sum 21 2 77 100 139 138 155 194 218 256 1, 350 
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TABLE 3b. 


Cosy 
| k 1 2 5 4 n 6 7 S 9 10 sum 
i 1 1 2 1 2 6 
- 1 3 10 6 12 10 If IS 16 a7 132 
3 7 6 5 22 14 22 31 35 27 73 
} 0 
BK 5 0 





7 1 
s 1 5 tb 
y ] 1 5 l 11 
10 2 if ] i) 
sum . 7 17 2b 24 1] 25 3s 14 5 64 342 





TABLE 4a. 

















cost 
| 2 3 } 7 ih 7 S 9 10 sum 
i 
1 6 10 th i) S r 2 l ou 
2 5 s “ Is 1 1 It 13 22 44 15s 
3 j 14 14 27 19 ta) it) 4] 3h 1s 261 
4 2 t Y } ty t) q th 3 4 Ee 
iE 5 1 10 s b 11 10) 5 s 4 7) a6 
6 1 3 12 ‘ 7 7 s 10 f ] jaY 
7 2 2 4 12 ri t j 4 5 } ay) 
s l } s 12 4 13 11 a) y 13 St) 
9 1 y 12 3 14 17 13 12 15 ht) 120 
10 7 15 29 55 44 7Y s4 Sd 119 121 637 
Sum 30) 79 114 151 137 19s 19s 140) 21S 234 1,549 
TABLE 4b. 
COSH 
} l 2 3 } 5 Hi 7 Ss 4 10 sun 
1 th , ) 7 t $8 
2 4 i) 2 14 st) ] 12 12 21 4:3 136 
3 3 s 13 ] y 2 27 28 21 $ 15S 
4 1 l 
I 4 l 5 ’ 
‘ 1 2 7 l 1 
“6 ps 3 7 $ 15 
s 1 2 2 2 4 l 
y 5 l 1 7 1] 2 27 
10 2 2 l 2 } } ) 17 
Sum ‘ 17 41 3S M4 39 62 15 12 12 Th i261 
TABLE 5a 
COS 
l 2 3 1 5 t s 4 WW) sun 
i 
l 6 24 2H 21 s 12 in| 7 1 ] 127 
2 12 20 24 1 11 l 17 24) $5 24 1] 
3 10 22 37 41 1 ta 34 31 10 10 IT 4 
t 3 } It th y 6 a re ’ - =; 
E 5, 4 s st) 13 20 10) th) 4 i a SS 
6 4 y 6 20 13 10) 12 4 2 SI 
7” 1 5 s i) 11 1s 16 7 8% | 3 SH 
& 5 5 17 10 14 14 10 21 s 15 119 
y 2 4 Is 12 12 7 27 22 25 12 141 
10 3 Is 27 17 5b Sf) ] 90 13 116 HOT 
Su $f) 114 190 195 212 219 227 228 234 192 1,857 
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TABLE 5b. 























cose 
1 2 3 6 7 8 9 10 Sum 
l 6 21 23 20 It +] Y 5 109 
2 12 17 13 14 10 i) 11 25 35 21 167 
3 S Ib 2 22 25 27 18 17 160 
4 3 13 1] 27 
5 2 2 t 6 12 26) 
6 l 3 10 6 20 
rf 3 ] 3 6 Y 5 27 
Ss 2 2 6 2 6 4 6 5 33 
Y 2 3 6 3 2 2 15 10 4 47 
10 l l } 1 1 3 9 4 24 
sum 35 65 95 SS SS 67 67 71 43 21 640 
TABLE 6a. 
Cos? 
1 3 ) 6 9 10 sum 
l 111 M4 61 i] 52 42 30 19 15 3 458 
2 15 33 4] Ath 20 30 35 30 24 Ss 322 
B} 17 53 57 $3 57 38 y s 4 3 319 
1 11 13 Is 14 20 24 5 $ 7 4 120 
5 11 H IS 15 15 13 22 2 1 |} 1 109 
t) Ss 13 10 s 11 18 20 2 1 96 
7 i 11 15 1 y 17 1S 26 S 1 120 
s iT) 6 i) It 16 15 11 21 14 4 121 
v 7 6 11 14 17 10 13 21 30 17 146 
10 7 15 28 4t 42 7 67 75 78 85 490 
Sum 260 237 71 At) 256 257 228 226 183 127 2, 301 
TABLE 6b. 
cos? 
) 6 9 10 sum 
1 107 XU 5S 39 18 39 29 1S 14 3 435 
2 1 30 39 52 17 24 34 28 21 a 293 
3 14 it) 17 37 18 29 l 252 
} i] 11 12 s IS 21 2 Sl 
5 i) s 13 S s 10 20 76 
6 2 i Ss 4 5 6 12 15 56 
7 3 7 y 5 6 10 s 20 3 7l 
s 7 3 6 i) 7 7 5 ll 11 66 
q 6 5 7 6 s 4 s 16 19 11 90 
10 4 3 Z } 5 t S 5 12 12 5Y 
sum 232 197 201 ‘2 170 154 127 113 SO 33 1,479 
TABLE 7a. 
eos? 
) 6 9 10 sum 
i 
1 0 
2 l 2 2 if) 4 Ss Ss 7 38 
3 2 1 3 13 | 14 33 
} 2 2 4 
5 1 2 2 5 
j 1 1 3 1 6 
7 1 1 l l i 
s 1 1 1 3 
4 2 1 1 3 7 
10 5 1 4 7 7 13 19 13 14 8&3 
Sum ! st] 3 12 16 21 26 34 31 30 183 








TARLE 7b. 





























cost 
3 4 5 if 7 q 10 sum 
‘ 
] 0 
2 2 2 s 11 14 v n2 
3 ] 3 l t ih) 24 54 
4 3 l 4 ! l 10 
1 2 3 l 1 2 10 
ti 1 2 1 3 2 1 ! 1] 
i l l 2 ] 2 1 ‘S 
S a) l 1 ; { 4 13 
4 l 3 2 $ 5 4 3 21 
10 s if 10 13 Is 2g 42 3S 1s 215 
Sum 12 4 2» a4 133 rl 7) Go X] 3g9 
TABLE 7( 
Co 
3 { ‘ q te) Srl 
1 0 
2 2 2 s 7 12 14 11 7 
3 1 1 » t 20 oT “Ho 
| 4 ] 5 l 5 1 2 13 
5 1 2 3 2 ] 3 12 
iD ! 2 3 2 3 2 15 
‘ 1 ] 3 l 2 ] u 
S 1 2 2 l 3 5 3 17 
9 1 l 3 2 1 6 5 6 28 
10 4 4 12 15 24 4] 51 Al 4 ra | 
Sum 13 12 27 $2 ; eS 87 ST 107 {ss 
TABLE Sa 
COsSH 
4 1 ty 7 ” 10 sum 
l () 
- 2 5 5 7 13 Is 1 0) 75 
3 1 2 l 3 ] a Is 22 a) 
4 2 2 ] 2 2 2 1] 
5 l 3 ] 3 3 3 14 
6 l 2 2 5 1 5 2 Is 
4 1 4 4 3 3 i 1Y 
s 2 5] 2 l 1 5 3 17 
4 l 2 2 ; s 6 7 5 34 
10 10 1] 17 23 27 1 Wy) is oY 4 
Sum 17 20 32 1} 63 S38 SI] a] 100 35 
TABLE Sb 
cost 
3 4 5 6 7 4 10) Sum 
l () 
, » bs 
2 2 if) 4 11 10 It) st] l ov 
3 1 1 1 1 2 6 11 22 415) 
4 2 1 2 é 4 2 l 3 14 
) v4 2 3 } 1 ] Ls 
6 l 1 1 1 l 4 4 ] 17 
‘ l ] 2 ] } 4 12 
s l | 3 4 l l 3 t of) 
4 1 l 2 3 2 } 3 f) if) 30 
10 5 10) 16 22 3 1 1] 53 6 293 
Sum s 1) 27 45 2 80) 71 &2 112 2 





TABLE 8&e. 























COs? 
ny be 9 10 sum 
| 0 
2 2 5 } 11 10 16 S 56 
3 3 Ss 17 28 
i 0 
5 0 
if) 0 
7 0 
S 0 
0] 0 
10 0 
Sum 0 2 5 } Il 10 16 11 S 17 S4 
TABLE 9a. 
cos? 
6 s 9 10 sum 
1 1 ] 2 | 2 7 
2 2 5 4 s } 6 7 S ll 55 
3 1 l 2 6 } 2 15 | 10 ) 50 
} ] 2 5 3 l 3 15 
4 ] 3 l 2 1 $ 2 } 2 20 
2 2 3 l 4 } 3 $ 2 25 
7 2 l 3 1 4 3 6 2 22 
S ] j by 3 7 5 3 4 9 2 47 
se) ] z 5 7 a] s ri 10 6 11 66 
ta) { s 17 t} 36 3S 45 56 &2 78 390 
sum 11 23 13 ) 78 68 73 101 130 120 697 
TABLE Yb. 
COSA 
2 6 S 9 10 sum 
1 ] l 2 | 2 7 
2 3 2 6 2 5 6 S 11 43 
3 | 2 3 l 10 7 6 30 
4 0 
5 0 
t ] l 
7 1 1 
s 1 ] a) Ss 
y 3 2 | 6 
th 1 » 6 8] 
Sum 3 1 14 Ss 17 5 6 16 15 V7 105 
TABLE 10a. 
Cos? 
6 Ss 4 10 Sum 
1 3 6 5 4 6 3 l 33 
2 3 l i } 5 6 4 S 10 15 62 
3 2 ] 5 5 } 6 11 12 7 5 5S 
4 2 3 } ] 3 2 1 4 3 23 
5 S 6 3 5 5 3 } ] 2 37 
i 2 Zz 14 3 % 5 5 5 3 2 1} 
7 2 3 a) +] 7 l i 5 4 3 43 
s 1 5 \ s 12 13 10 y i) 5 sO) 
iF) 3 2 ) 10 10 +] 10 17 16 S3 
10 2 4 IN 27 34 ‘4 70 77 92 92 480) 
sum 17 it) 73 i3 S7 108 120 131 147 143 945 
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ABLE 1O0b. 























ecosd 
7 st) 10 sum 
1 3 5 5 3 4 h 2 1 29 
2 Z l 2 3 3 3 § in 10 15 1s 
3 2 l 3 6 9 7 2 l 31 
H Z y 4 
5 4 3 7 
6 1 2 S 3 14 
rf l 2 2 7 2 14 
S l 3 5 3 5 i 19 
4 2 2 5 7 7 
10 1 4 3 ; z lt) 
sum 12 0) 24 26 4 29 20 16 12 16 107 
TaBLE Ila 
4 10 | 
l 5 st] 12 14 10 10 6 *) 7H 
2 10 13 10 10 z i] 12 Is 13 107 
3 $ 12 20) 4 1 11 t) 4 0 1 a2 
4 2 2 2 ia, 6 1 2 2 | 3 jb 
5 l 1 i 10 22 if) 7 ; 3 2 2 
H 11 7 7 22 13 0) ‘ ' 4 at) 
7 ) 4 sS 7 1Y 10 ] 4 5 77 
s ] 5 17 1 14 1s 17 4) 11 10 128 
8] 3 4 4 14 7 14 1 29 14 1Y 142 
10 1 11 28 thy 42 Hb 74 Y2 Bat) 115 54 
Sum 27 74 128 a0) 153 170 17 IST Ist 175 1,417 
TABLE IIb. 
7 ot) 1 sum 
1 5 4 1] 13 4 7 t) } 64 
2 10) 10 3 % 6 } 3 12 I 12 7s 
3 3 12 Is 5 s ry $ iD 
} 2 l 10 Is 31 
5 1 F 4 $ 5 13 20 
6 1 2 16 s 32 
7 $ l 3 7 1] i) 24 
S ] ] ty 3 7 } 10 37 
4 ] 2 } 4 ] 4 13 7 ] 37 
10 3 | 2 2 3 1 12 7 34 
Sum 23 44 ail 61 “7 17 1 + » 12 131) 
TABLE 12a 
co ‘ 
7 1] tp Sul 
1 s| 62 50 aay dQ ss 10 15 Ww s 416) 
2 38 30) 38 57 19 13 If} 16 14 11 252 
4 1Y 23 22 30 4} Is $ 3 1 3 171 
} 17 20 14 If 2 3S y 2 145 
5 16 20 19 22 17 20) 29 2 l l 147 
f 7 12 15 12 10) 14 33 20 5 Z 135 
7 7 17 1 12 15 22 17 36 7 2 14 
x 1] 9 13 19 21 IS 19 28 20 2 160) 
ua 13 10 16 11 HA) 1] 3 Pa! 48 Is ISS 
It i 14 23 38 44 v4 73 71 43 GY Slt 
sul 213 22 234 245 247 $] 240) 224 194 146 2, 184 





TABLE 12b. 























cos9 
k 1 2 3 1 5 6 7 8 9 10 sum 
! 
1 77 60 50 26 27 25 10 15 10 8 308 
2 34 27 35 53 1S s 13 14 13 10 225 
3 18 20 20 29 45 (ie ee = 4 150 
4 $2) } 17 15 15 25 37 5 —_ 129 
/ 5 14 17 16 14 11 15 29 p e 116 
6 6 y 10 6 7 9 25 21 93 
5 10 15 9 10 14 s 31 ee 107 
S 4q 3% s s 10 10 9 19 17 93 
y 11 8 10 6 i) 6 14 20 28 ll 123 
10 2 6 4 s 6 7 12 of) 20 21 95 
Sum. 191 177 183 174 168 148 126 129 93 50 1, 439 
TABLE 13. 
COSA 
| 1 2 3 { { 5 6 7 8 9 10 sum 
? 
] ) 15 b Ss ) b ) 2 1 3 56 
2 t) 5 7 13 5 10 10 6 12 21 95 
3 5 10 9 ll 7) 5 9 6 5 10 79 
j 1 3 S 3 ) 6 2 3 1 33 
| 2 3 ) 9 4 4 2 2 4 42 
ti » 6 7 ti 5 14 12 s 6 4 70 
7 1 3 s 7 11 10 10 9 6 7 72 
S 2 . 3 10 9 15 20 10 7 4 87 
’ 3 4 5 6 13 19 7 1S 10 13 YS 
lt h 10 4 } 3S 59 62 57 OS 106 481 
Sum 33 70 S) 90 110 148 141 121 147 173 1,113 
TABLE l4a. 
cosd 
1 2 % $ 5 6 ij S 9 10 Sum 
] 0 
2 3 $ 11 5 15 16 15 11 3 82 
] ] l 3 2 2 S 13 23 54 
} ] 1 2 4 3 2 3 16 
| l 2 l 2 2 1 2 11 
t 1 l I 5 2 1 7 1 4 23 
7 2 3 I 3 l 3 5 18 
s 2 2 2 3 | 5 3 Is 
4 1 2 2 j 6 | a) 7 30 
10 a, 1 6 10 14 22 30 31 50 50 219 
Sum (} 4a 99 25 16) 53 O4 66 83 | 97 471 
‘TABLE 14b. 
cCos4# 
- 
k l 2 3 4 5 6 7 S +) 10 sum 
| 
1 5 0 
» { 3 12 6 if IS 16 17 3 ; 5 
3 ] ] 1 3 4 2 10 | 16 28 66 
{ ] | 2 4 3 2 3 16 
/ ¢ ] 2 1 2 2 1 3 12 
( 1 1 1 5 2 3 7 2 4 26 
7 ] 2 3 1 3 3 3 5 21 
s 2 2 3 3 1 1 5 3 20 
) 1 3 3 4 6 3 5 8 33 
lk 2 i 7 10 14 22 $1 33 l 51 225 
sum 7 10 24 27 40 57 | 6S 72 SS 105 514 
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TABLE 15a. 


cosé 























» 5 7 i) 10 Sum 
1 0 
2 2 2 3 4 2 ‘. 3 *4 23 
3 2 3 1 h 21 33 
$ 1 1 ] P4 2 2 i 14 
5 2 3 l 2 s 
‘ l 3 2 3 2 i | 16 
7 ] ! 3 1 3 1 14 
& 1 } ] 7 3 3 3 22 
4y 1 2 3 4 3 u 10 34 
10 2 ) 6 14 17 2b 17 23 28 21 159 
Sum 2 10 15 25 30) +4 35 | 5S oo) 320) 
TABLE L5db. 
Cos 
, ; 9 10 Sum 
1 0 
2 ] 2 3 5 s 2 5 6 3 1 3H 
3 l 2 t | l 2 1] 27 14 
4 l | l 3 i i 2 } 21) 
5 4 4 1 3 3 2 4 19 
6 1 3 1 4 3 3 4 1 21) 
7 2 2 t 4 | 3 2 Is 
S 1 l 4 1 7 3 4 3 2 26 
y 3 3 3 6 3 10 12 4() 
10 2 6 7 10 1S 24 2) 24 30 27 182 
sum 3 12 1Y 38 42 53 $5 50) iS SO 10 
TABLE 15e. 
cose 
2 7 9 10 Sum 
i 
l {) 
2 3 4 4 b s 2 7 4 ] 44 
3 ] 4 5 ! 2 1 4 14 33 63 
4 1 1 l 3 $ 5 2 3 20 
5 ] 3 t 1 3 3 3 3 21 
i) l 3 2 } oe 3 4 20 
7 2 2 } 4 2 3 v4 1Y 
s ! ] 3 1 4 3 2 15 
Y 4 3 2 4 2 5 7 25 
10 <4 7 6 14 15 4) 1s 1 28 24 165 
sum 5 15 1¥ 32 34 A 40) 51 66 75 302 
TABLE I6a 
2 7 4 IK Sum 
] i) 
2 3 { } t $ i 4 l 2G 
3 4 1 ] 7 Ps | 33 
$ 3 | » (} 
5 l l l ! l 1 6 
t 2 1 3 s 10 10) 4 lth 4 7 sy 
7 ] l 4 ( S s 11 44 
S l ] 2 2 l 2 4 a) 23 
4 | | 3 2 F 5 4 6) $2 
10) 1 3 2 11 14 15 1 13 25 22 12 
Sul 4 s 12 $2 3u 42 17 1 4 75 372 
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TABLE 16b. 

















cosé 
} l 2 5 } ) 6 7 Ss a 10 Sum 
1 
1 0 
2 1 6 10 7 16 10 13 13 2 78 
3 1 2 7 13 13 3Y 
} ] ] 2 1 2 1 8 
l 4 1 2 8 
t 2 | 3 s i) 7 s 17 13 | 10 78 
7 2 i 2 Q 6 10 5 | 14 52 
s 1 6 3 2 6 7 4 6 35 
y 1 ) 6 6 6 1 4 y 38 
10) 4 4 3 12 11 a3 20 23 27 28 153 
Sum_. 7 12 1Y 413 50 60 65 Ss) 72 SI] 480 
TABLE 16ce. 
COSA 
k ] 2 3 4 Fi 6 7 8 9 | 10 Sum 
i | 
1 i 0 
2 5 S 11 12 14 if) 13 13 4 2 SS 
3 ] 3 +] 4 3 5 14 34 56 129 
$ = 1 3 1 8 | 5 2 4 26 
5 1 3 > 3 4 3 45 4 24 
fi 2 4 2 4 3 4! 10 1 28 
7 2 3 4 4 2 4 4 2 
S 2 1 3 4 3 ] 5 4 23 
4 l 2 2 3 5 2 4 S y 36 
10 2 1 2 13 10 15 22 15 | 21 | 16 117 
Sum 10) 14 20 44 44 $5 64 | 61 92 | 100 494 
TABLE 16d. 
coseé 
l 2 3 { 5 6 7 S i) 10 Sum 
l 0 
2 $ 4 4 ih) $ 4 2 l 27 
3 4 3 l 1 7 21 33 
4 3 | 2 6 
1 1 1 I 1 1 6 
6 l l 1 3 l 2 3 1 13 
7 ! 2 1 | 2 2 8 
s l Y 4 l l 1 6 
st) 1 1 3 1] | 4 } 3 14 
10 1 l | S +) 1] 14 S 15 | 11 79 
sum 2 5 5 17 21 25 25 17 33 42 192 
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A parallel band at 


been observed under high resolution and analysed. The 
revealed the presence of l-type doubling in the doubly degenerate excited state. 
analysis of the parallel band it is found that Fo 
active 


been observed near 2,600 em-!. Both 
Ivsed, 
vious determinations. 


1. Introduction 


A study of the spectra of deuterated methanes 
was undertaken as part of a program for studying 
the spectra of deuterated molecules. In spite of the 
apparent simplicity of the CH3,D and CD;H mole- 
cules, both beng svmmetric tops, very little work 
has been done on their spectra. Early work on thé 
spectra of these molecules was done at rather low 
resolution {1,2].2 More recently one band of each 
molecule has been studied under somewhat higher 
resolution by Thompson and his coworkers [3,4]. 
Two bands of CDH have been studied in the photo- 
graphic infrared [5| and the overtone of the C—H 
stretch, x, has been studied with high resolution [6]. 
A comple te study of the spectra of the two molecules 
has been made using prism instruments by Wilms- 
hurst and Bernstein [7]. 

The resolution available with the instruments of 
the Radiometry Section made it feasible to study 
several bands of these molecules between 3 and 6 yw in 
greater detail than had been possible before. No 
perpendicular band of either molecule had previously 
been studied with high resolution. In this work a 


perpendicular band of CH;D was observed and 
analysed. Likewise a hybrid band of CD,H was 
observed and analysed, one component of which 


obeys the selection rules of a perpendicular type 


band. 


2. Experimental Method 


The spectra were recorded on a 2.35-m_ grating 


spectrometer. For the wavelengths shorter’ than 
4 uw the spectrometer was equipped with a 10,000 
linesin. grating which was double-passed, and: 


cooled PbS detector. Higher orders of the grating 
were eliminated by using a cut-off filter that started 
transmitting at 2.7 uw. The resolution in this region 
is 0.02 to 0.03 em-'. For the CH,D band at 2,200 
’ l . stn , "a0 A 1 rd itl , PbTe 
em! the spectrometer was equipped with a Pb T% 
detector and the grating was not double-passed. 


The resolution in this band was ~0.05 em~'. The 
The work reported herein was supported by the U.S. Atomic Energy Com- 
Mission 


? Figures in brackets indicate the literature references at the end of this paper. 


components, A 
The ground state /’y value found from this analysis is in good agreement with pre- 


1959) 


2,200 em-! and a perpendicular band at 2,780 em~! of CH3D have 


analysis of the perpendic ular band 
From the 

A hybrid band of CD3;H has 
and FE are observed and ana- 


3.880 em-!, 


spectra were measured by using the fringes of a Fabry- 
Perot interferometer as prev iously described [8]. 

The CH,D and CDH were obtained from Merck 
and Co., Ltd., and had a stated purity of 98 percent. 
The gases were used without further purification and 
no troublesome impurities were found. The gas was 
placed in a cell with a pathlength of 6 m and the gas 
pressures used ranged from a few millimeters to a 
few centimeters of Hg. 

Two regions of absorption of CH;D were studied 
in detail. The parallel band at 2,200 em~! which is, 
in Herzberg’s notation [9], ™, and a_ perpendicular 
combination band near 2,780 cm~! which has the 
excited state v3+v;. Recorder traces of the observed 
spectra of CH3;D are shown in figures 1 and 2. 

One region of CDH was studied in detail. This 
region is of particular interest since the absorption 
arises from a hybrid band. Both infrared active 
components of this band are resolved. The observed 
absorption is shown in figures 3 and 4. 


3. CH;D 


Parallel band. The absorption in the 2,200 em 
region is clearly that arising from a typical parallel 
band. Unfortunately the short wavelength portion 
of the absorption is masked by the absorption of 
atmospheric CQ,. In spite of this handicap it is 
possible to identify transitions as high as (10), 
although some of the A- components are not observed 
because of the atmospheric absorption. The ob- 
served wave numbers and quantum number assign 
ments are given in table 1. Some of the assignments 
are indicated on figure 1. The intensity alternation 
of the A-components arising from the nuclear spin 
statistics of the three identical H-atoms is clearly 
seen in figure 1. The statistical weights of the A 
and / rotational levels are in the ratio 2:1, and the 
observed intensity alternation simplifies the process 
of assigning quantum numbers. The rotational- 
vibrational energy levels of a symmetric rotor in a 
nondegenerate vibrational state are given by 


F(J,K)=BI(J+1)+(A—B) kK? 
+ DyJ?(I+1)?+DigI (J +1) KPI+Dgh*. (A) 
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Frevre 1. The parallel band, v2, of CH,D at 2,200 emt. 


This energy expression together with the selection 
rules 


r 


AK=0 Kk +0 
AK=0 K=0 


AJ=0,+1 
AJ=+1 


enables the caleulation of the difference relations 


involving only the ground state energies. 


1,A)—P(J+1,AK) 
(Ze T 1) 


-4D5( J+ 


w he _ sv 
AF’ Said -(2B 

] < ” ro 

=) 4D, Kk’, (2) 


3D)5) 


and a similar one involving only excited 


energies, 


, RiSI,K)—POS,K) 


A _(9 2’ __ J)’ 
Al 2 J 1) 28 3D35) 


4D;(J+5)?—2DixK?, (3) 


which were used to determine the B’s and centrifugal 
distortion constants. Since the transitions in the P- 
and R-branches with A=0O, 1, 2 were not resolved 
these transitions were not used in the calculations. 
Thus 35 equations of type (2) and 29 equations of 
type (3) with A>3 were solved by the method of 
least squares for the best estimates of the constants. 
The combination sum, 


R(J—1, K)+P(J, K)=2%+2[(A’—B’) 
(A’’— B’’)|K?+2(De—Dr) K'+2 [(B’— B"’) 
tL (Dix—D jx) K |S? +2(DS—D5) J?(J+1)?, (4) 


was used to determine the band center. The results 
of these determinations are given in table 2. These 
molecular constants compare favorably with pre- 
vious results [4] although a different method of analy- 
sis was used here. The present results are believed 
to be more precise and are more complete than those 
reported previously for this band. Earlier workers 
neglected the effeet of Dyx whieh accounts for the 
apparently different values of (A’— B’)— CA’ — Bb’) 
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Figure 4. The ®Rx branches of the -component of the 2v; band of CD3H. 


TaBie 1. Observed wave numbers of 2,200 em band of CH;D in the P- and R-branches, a result they noticed but 
eetemes did not explain. Starting at P(10) there is a pertur- 

Jr R JK R P bation which affects the A=9 component. This 

: : arises from an interaction of the excited state with 











Oo. -- 2, 207. 70 . ear ar | a aa os another higher vibrational level of CH;D. Some of 
” 2, 192. 29 the transitions to this other level can be seen in 
2, 215. 32 9 rp 2, 127. 29 . Y . . : 
% Yosh figure 1. Unfortunately not enough of this overlap- 
2 oo) wees 2s 6(ping band can be identified to enable a complete 
2 2, So. Se U4 a, led. de e » ° ° 
5 2is2s understanding of this perturbation. However such a 
3 2, 230.23 | 2,176.53 || 9 2, 128. 72 : ; 2 i . 
eerste ee Saees ee 34299, perturbation cannot have any effect on the ground 
3 2, 200. 54 9 2123 state constants as determined by the method out- 
’ 2,237.56 2,168.50 100.1 2279.48 21883 lined here. The perturbed levels were not used in 
4 2, 237.66 | 2,168.63 102 2,279.58 | 2,118.95 | eS } oe ls , ; : 
{ 2 937.85 2,168.84. 10 297980, Zig determining the excited state constants. 
4, 2, 238. 13 104 2, 280. 04 2,119. 44 
105 2, 119. 80 TABLE 2. Constants of CH;D from v2 at 2,200 cm-! 
2, 244. 78 2, 160. 44 10; 2, 280). 75 2, 120. 26 NS Se ad oe ee 
2, 244. 89 2, 160. 55 10; » IRL YI 2 120.78 == 
2,245.09 2,160.77 10 2, 121.38 Excited | Ground 
2, 245. 34 2, 161.04 109 2, 121. 9 
2, 245. 68 
es se Ho, 1 a 4 B __| 3.8372 em-! 3. 8809 cm 
‘ ) 4 ») by yA ‘ ) . 2 ) £ | 
6: ean 2, 152.39 TE 2,110. 60 ( J --. 0.000055 0. 00005 
6 2, De. 2d 2, 152. 59 lly 2,110. 90 Dyk OOO12 _ OOO12 
6 2252.48 21ST 1s 211. 26 
th: 2, 252. 80 2, 153.23 11 2,111. 70 
tie 2, 253. 21 11; 2, 112. 23 (A’— B’) —(A”— B”’) =0.03858 vo= 2200.03. 
Mk 22.4 
: face | ce oe ae Perpendicular band. The absorption near 2,780 
z 2,259.31 2,144.35) cm—' arises from the perpendicular band y;+ y;, and 
7 2, 144. 63 12; 2, 101. 62 * ° . ‘ rl ns . : 
; 2,250.84 2144.98 12, 2110177 1S shown in figure 2. The excited state of this band 
c oe | ee |e 3100-9 is doubly degenerate and the rotational contribution 
‘ -- ) ‘ amt ~s, meme . e* : . . rm 
12 2:02.65 Is not given by the simple expression (1). The levels 
s 2, 260. 90 2,135. 68 12 2 103. 10 5 ; ‘ : ‘ 
8 2266.02 | 2135.81 |, 12. Size are perturbed by the interaction of rotation and 
8 2, 266.20 | 2,136.03 || 125 210425 vibration which leads to an additional term in the 
s 2, 266.4 136. 31 12 2 105. 48 
. 2,136.37 | 12 2105.75 energy, thus, 
igs. Ky) EBLE J T 1) T (A B)k°F2ACh 
centrifugal terms. (5) 
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For a perpendicular transition the selection rules on 


K are somewhat relaxed, 


a 


AJ=0, +1 AK= +1. 
The overall band is thus made up by a series of sub- 
bands. [9, pp. 425 and 431] The prominent feature 
of the overall band is the series of Q-branches of the 
subbands whose spacing approximately  2[A’ 
(1—¢)—B’]. 

The series of (@-branches characteristic of a perpen- 
dicular band are clearly evident in the 2,780 ¢m7! 
region. Again the nuclear spin statistics provide the 
clue which enables the assignment of K-values to 
the Q-branches. Those @-branches with ground 
state A’s 2 multiple of three are stronger due to the 
nuclear spin statistics and this intensity pattern is 
clearly evident in the figure. This is the first per- 
pendicular bend to be observed in which all the 
individual (V-branch transitions are resolved. Once 
the quantum numbers have been assigned in the 
Q)-branches, each subband can be considered sepa- 
rately and effective constants can be calculated from 
the Q-branch using formulas applicable to linear 
molecules, 1.e., 


IS 


Vv Veub (B’ BY’) SIS T 1) T (Di — D7) J? (J T 

These effective constants together with the ground 
state constants determined from the parallel band 
can then be used to calculate the ?- and R-branches 
of the subband associated with the particular Q- 
branch in question by using 


TABLE 3. Subbands of CH3D 
kK=0 K=1 K=2 
J 
AK=1 AK=!1 AK 1 AK=1 
AJ=0 
1 2, 780, 42 2, 774. 85 
2 2, 780. 20 2, 785. 96 2, 774. 70 
3 2, 779. 82 2, 785. 74 2, 774. 48 2, 791. 43 
H 2, 779. 35 2. 785. 44 2, 774. 19 2,791. 11 
5 2, 778. 78 2, 785. 11 2, 773. 85 2, 790. 75 
t 2, 778. 13 2 784. 78 2,773. 45 2, 790. 34 
2, 777. 40 2, 784. 43 2, 773. 00 2, 789.95 
s 2, 776. 61 2, 772. 52 2, 789. 57 
+] 2. 140. 82 yA FG Pa 2, 789. 23 
10 2, 775. 04 2 til, ae 
AJ=+1 
if) 2 4 ee 
1 2, 795. 99 2, 801. 48 
» 2,803. 64 2, 809. OD 2, 797. 76 2,814. 70 
3 2,811. 28 2,816. 47 2, 805. 26 2, 822. 15 
H 2, 818. 93 2, 823. 89 2, 812. 65 2, 829. 51 
- 2 S96. 61 2,831.35 2, $20.00 2, 836. 85 
ti 2, 834. 32 2, 838. 71 2, 827. 29 2, 844. 20 
\J ] 
l 
2 2, 764. 99 2, 759. 33 
3 2, (ot. 15 2, 751. 44 
4 2, 749. 35 2, 754. 70 2, 743. 47 
5 2, 741. 47 2, 746. 67 2, 735. 43 2, 752. 36 
6 2, 733. 63 2, 738. 61 2, 727. 34 2, 744. 27 
7 2, 725. 85 2, 730. 54 2, 719. 20 2, 736. 12 
s 2, 718. 12 2,711.05 2, 727.97 
u 


pe 


New 


' Wwe 


NNN 


NN te te bo 


NNN 


Your + (B’+ B’’)m-+ (B’ — B’’)m?+2(D, 


yp 

+ D7)m3+ (D5—D7)m*, (7) 
with m=—J for the P-branch and m=J+1 for the 
R-branch. 


It should be remembered that in a given subband 
the branch in which AJ=AK is the stronger and 
becomes very much stronger than the branch 
with AJ#AK as K increases. In this manner as- 
signments could be made for all the subbands except 
that for which A=0. The observed wave numbers 
and quantum number assignments made in this 
manner are given in table 3. Unfortunately the short 
wavelength end of the band is overlapped by a 
much stronger band of CH;D so the weaker perpen- 
dicular band could not be analyzed in this region. 
The /?- and /f-branch transitions of the A=0 
subband were located by a trial and error procedure. 
Since both branches are fairly strong in this mole- 
cule they could be picked out with a minimum of 
difficulty. 

Ground state constants and band centers were 
determined using equations (2) and (4) with A=0. 
The ») found from (4) agreed exactly with that 
found from (6) but the values of B’—B’’ were 
different by nearly a factor of 3, that found using 
(6) being larger. This behavior is reminiscent of 
the behavior of a perpendicular transition of a linear 
molecule where the levels of the degenerate state 
have the so-called l-type doubling, one component 
of the doublet being accessible from the ground 
state with AJ=0, and the other component being 
accessible from the ground state when AJ=+1, 





rpendicular band at 2,780 cm7! 
kK=3 K=4 K=5 K=6 
=—] AK=1 AK=-1 AK=-1 AK=-1 AK=-1 
769. 15 , 
768. 91 2, 763. 40 
768. 64 2, 796. 86 2, 763. 14 2, 757. 66 e 
768. 33 2, 706. 44 2, 762. 82 2, 757. 38 2, 751. 91 2, 746. 15 
767. 98 2, 762. 47 ys 2, 751. 59 2, 745.85 
767. 60 2, 795. 54 2, 762. 09 z 2, 751. 27 2, 745. 50 
767.18 2, 795. 11 2, 761. 64 2, 2, 750. 92 2, 745. 19 
766. 73 2, 794. 70 2, 761.14 2, 2, 750. 61 2, 744. 90 
766, 23 2, 760. 58 aul 2, 750. 30 F ; 
792. 19 
799. 6S 2, 827. 87 2, 794. 16 
S807. 12 2, 835. 20 2, 801. 62 
S14. 50 2, 842. 51 2. 809. OO 
$21. 83 2,816. 29 
753. 76 
745. 85 2, 740. 35 
737. 91 2, 732. 39 2, 726. 91 
729. 89 2, (58. 13 2, 724. 40 2, 718. 91 2, 713. 46 
721. 83 2, 749. 96 2, 716. 34 2, 710. 89 2, 705. 43 2, 700. 03 
713. 76 2, 741. 78 2, 708. 27 2, 702. 85 2, 697. 39 2, 691. 97 
705. 64 2, 733. 63 2, 700. 14 2, 694, 77 2, 689. 32 2, 683. 91 
2, 681. 27 
» 


673.17 
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thus leading to a different effective Bb’ value for the 
Q-branch than for the R- and /-branches. 

The possibility of this type of behavior in sym- 
metric top molecules was first suggested by Wilson 
(9, 10]. The interaction to a second order approxi- 
mation has been worked out by Nielsen [11] who has 
found that for molecules with C3, symmetry the 
interaction matrix elements are of the form 


(v,,l;,,K\Aelr,l;42,K42= [ei Fl) (£14 2)] 
<[(J(J+1)—-A(K+1)|[(J(J+1)-A(K+1 )(K+2)]}? 
(7) 


As can be seen levels are mixed pair-wise by this 
interaction but only the levels A=1, /=1; A=—1, 
/=—1 are actually split. All other levels are merely 
shifted in position by this interaction to this approxi- 
mation. For K=1, /=1, (7) reduces to 


glv+ 1) J(J+1) (8) 


which is the same as the expression for a liear 
molecule. The actual splitting is twice the quantity 
given in (8). Since the interaction in the A—0O 
subband varies as J(J+-1) it has the effect of pro- 
ducing two effective B’ values in the subband. Thus 
from the two effective B’s found for the excited 
state with A=1 it is possible to evaluate B’— B’’ 
and qg for this state. These results are included in 
table 5. 

Sufficient lines in the P- and /?-branches of 6 
subbands were observed to enable a determination 
of subband centers from the transitions using the 


relation 


ER .(J— l ) 5 | £P AJ) hia 


(B’— BY’) F?+ (D’—D") J (9) 


since AK is a constant for each subband. The highest 
value of AK for any subband is 3. Similarly the bend 
centers of the same subbands can be determined from 
(6). The vo’s obtained by these two methods «gree 
within the experimental error end ere compared in 
table 4. 

These subband centers may then be used to deter- 
mine the quantities in the relation 


2¢) — B42 


[(CA’— B’) —(A”’ 


w+[A’(1 
BAK 


Veub 


[.’(1 ¢) (10) 


BY’ )| kK 


This relation neglects the effect of centrifugel dis- 
tortion but this should not cause serious error since 
only values of A<3 are used. The main error will 
erise because nowhere is 1), considered. However 
since no high quantum numbers could be used in 
the analysis there was no good way to determine 
this constant with eny precision for the excited state. 

In view of the perturbation in the excited state, 
the usual combination relations among the transi- 
tions of a perpendicular band could not be used in 


the anelvsis of this band [12]. However, if the 
perturbation discussed above is the only perturba- 
tion present, the combination relation 


TABLE 4, Subband centers of v3+-v5 band of CH,D 
i AK Q-branch R- and P- 
branches 
0 | », 780. 55 2. 780. 55 
l 1 2,786. 15 2, 786. 17 
] -1 2, 774. 90 2, 774. 88 
2 1 2, 791. 90 2, 791. 90 
4 ] 2, 769. 33 2, 769. 31 
3 l 2, 797. 73 
3 l 2, 763. 76 2, 763. 75 
' —1 2, 758. 30 
Des 1 2, 752. 70 Siew 


T "Os e) 24 Vy 4 [A’(1—2¢)— B’] 


[PO(T x) 
B’') I(J+1)+ (D’—D") F(T +-1)?, 


(Bb (11) 
AV constant, averages out this perturbation and if for 
each Av one plots 


D”’) J?(J+1)? 
versus J(J+-1) 


(POT) + KOS) | 2 (D)’ 


a straight line should result. Such a plot for three 
values of A gives good straight lines and hence seems 
to rule out any further perturbations. 

The value of g obtained from the A=0 subband 
does not appear to be entirely satisfactory for the 
other subbands. However, not enough transitions 
in the other subbands are observed to enable one 
to make meaningful higher order corrections. The 
constants derived from the analysis of this band are 
collected in table 5. Unfortunately the inadequacy 
of the data makes these constants subject to uncer- 
tainties which are larger than would be normal for 
measurements of this precision. The strong over- 
lapping in the high frequency portion of the absorp- 
tion sharply limits the number of “PR and ’R lines 
which can be assigned and the general weakness of 
the “? and “7? lines makes it impossible to assign 
many of these transitions. However it is felt that 
the analysis is correct in all essential features and 
that only the evaluation of the small higher order 
corrections needs improvement. 
vs band of CH,D 


Papier 5. Derived constants from the vs 4 


em! 
vot+[A’(1—2¢) — B’] 2,780. 55 
vi1—¢ By’ 2. &2 
A’— B’ A” — B"”) 0. 02 
RB” 3.88 
B’-—B" —(0. O41 
D}—D; QOOTI 
q O11 


If A’ is calculated from the methane geometry, it 
is then possible to evaluate ¢ and vp. Using A=5.24 
cm! one finds ¢ 0.26. Simple theory predicts [13} 
that the ¢ of this band should be the same as for the 
fundamental v;. Herzberg [9] has estimated ¢; 

0.27, in remarkable agreement. However, it is 
felt that this may be in large measure fortuitous. 
The above figures lead to a value of v =2,776.33 


cm, 
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4. CD.H 


Hybrid band. Only the absorption of CD3H near 
2.600 em~' was studied in detail. This region has 
more than usual interest. The main absorption 
here is due to the first overtone of v5. v5 is a degenerate 
fundamental and its first overtone has the compo- 
nents A, and /7 [9,14]. 

The electric dipole moment selection rules allow 
transitions from the ground state to both the A; and 
FE components of this overtone level. Thus one 
would expect the absorption in this region to show 
the features of both a parallel and a perpendicular 
case. This is, indeed, the case as can be seen in 
figures 8 and 4. There is a collected Q-branch near 
2.565 em~! which is characteristic of a parallel type 
band. It is possible to identify the ?- and ?-branch 
transitions of this band and they are labelled in 
figures 83 and 4. Unfortunately the absorption is very 
weak so the A-components of the ?- and /?-branch 
transitions could not be identified but only the 
general absorption peak representing the aggregate 
absorption of all the A-components. Thus precise 
molecular constants could not be determined for this 
component of the overtone. The band center of this 
component is 2,564.6 em~!'. The bo value is com- 
patible with the previously determined values, 3.278 
em~! and the B’—B” value is about 0.01 em~'" 

The series of (-branches of the subbands belonging 
to the /-component are clearly evident in figure 3, 
the strongest being near 2,590 cem~!. This strongest 
Q-branch belongs to the A=0 subband, and the 
other A-numbering follows automatically, since again 
for this molecule the subbands with A’ a multiple of 
3 are strong due to nuclear spin statistics. In this case 


TABLE 6. 
k=! ke=2 


Ah=+1 AK=-1 AK=1 AW 


AJ=+1 

0) 2,599.19 

] 2 605. 76 

2 2, 612. 36 

3 2, 618. OS 2.613. 06 2, 624. 40 

i 2:'625. 58 2,619. 70 2, 630. 39 

5 2, 632. 20 2, 631. 57 2, 626. 33 2, 637. 59 

( 2, 638, 84 2,638. 16 2, 644.19 

7 2, 645. 52 2,644.79 2,650.81 

s 2, 652. 17 2, 651. 43 2, 657. 41 

9 2, 658. O07 2, 664.03 

10 2, 664. 72 2 670. 66 

1] 2, 671. 38 2, 677. 31 

12 2. 678. OF 2, O83. 9S 

13 

14 

15 

AJ 1 

l 

3 2, 572. 99 2. 567. ON 2 561. 21 

i 2, 572. 45 2, 560, 58 2 554. 71 
2, 559. OS 2, 565. 96 2, 554. 09 2, 548. 24 

b 2, 559. 50 2, 547. 62 2,541.80 

7 2, 541. 21 2,535. 39 

Ss 2, 534. 83 


Lines not listed are obscured by other absorption in this band. 


the ratio is 11:8. Unfortunately the structure of 
these Q-branches cannot be resolved well because 
the band is very weak, and if the pressure in the 
absorption cell is increased the pressure broadening 
eliminates the fine structure. Also the quantity 
(B’ — B’’) is very small in this transition ~ 0.007 em=!. 

The assignments of the P- and R-branch transi- 
tions were made by calculating the approximate line 
positions using the relation 


Yo=Vsew t (A 1—$) — B’]K+[(A’+ B’)—(A’’— B”)] 
Fk? + (B’+ B’’)m+ (B’— B”’)m? (12) 


Here the negative sign of the AK dependent term 
applies for AK=—1 and m=—Jd for the P-branch 
and m=J+1 for the R-branch. The coefficients of 
K and KA? were estimated from the Q-branch posi- 
tions and the coefficients of m and m? were estimated 
from previous work on this molecule [5, 6]. The 
resulting assignments are given in table 6 and shown 
in figures 3 and 4. Assignments could be made in 
11 subbands. Unfortunately, due to weakness of 
the absorption only lines in the branches AK=AJ 
could be assigned with any degree of assurance and 
are the only ones given in the table and the only 
ones used in subsequent calculations. In the A=0 
subband most of the ?-branch lines were masked by 
other absorption and could not be uniquely assigned 
to absorption peaks. The measurements could not 
be extended to longer wavelengths because of the 
weakness of the absorption and the fact that the 
sensitivity of the PbS detector falls off rapidly 
bevond 3.5 uw. Since the Q-branch transitions were 
not resolved and since, in general, only the stronger 


Subband assignments for the E-component of the hybrid band of CD3H 


K=3 K=4 K=5 K=6 


Ak=1 AK=-1 Ah=1 AK=-1 AK=1 AK=-1 AK=1 


2, 637. 12 
2, 643. 70 2, 649, 89 
2, 650. 29 2, 656. 48 2, 662. 74 
2, 656. 89 2, 663. 06 2, 669. 31 2,675. 64 
2, 663. 48 2, 669. 64 2, 675. 86 2, 682. 18 
2, 670. OS 2, 676. 21 2, 682. 42 2, 688. 71 
2. 676. 69 2, 682. 80 2, 683. 99 2. 695. 27 
2, 683. 32 2, 689. 39 2, 695. 55 2, 701. 80 
2, O89. 93 2, 695. 97 2, 702. 12 2, 708. 36 
2, 696. 55 2, 702. 58 2, 708. 67 2, 714. 90 
2, 703. 19 2, 709. 18 2, 721. 42 
2, 709. 83 2, 727. 99 
2, 555, 42 
2, 548. 94 2, 543. 21 
2, 542. 48 2, 536. 74 2, 531. 09 
2, 536. 06 2, 530. 28 2, 524. 63 
2, 523. 86 2, 518.19 
2,511. 78 


2, 505. 39 
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branch of the subbands could be observed it was not 
possible to use the usual combination relations for 
2 perpendicular band in the analysis. Thus an al- 
ternate and less desirable method of deriving the 
molecular constants had to be used. This method 
involved the use of an expression for the line position 
deduced from the energy expression for the two 
energy levels in the transition. In order to include 
the effects of centrifugal distortion an expression of 
thirteen terms had to be used. One would expect to 
fit almost any set of observations with this many 
disposable constants and the use of this expression 


is justified only by the results obtained. The 

expression used. is 

v=wt+ A’ —2¢)— B’]+bk4+ch?*+dkh*4+ eh! 

+ fm +- gm? + hm + im*+jmb+ km ke inet ni? he 
(1:3) 


in which the coefficients are readily related to the 
molecular constants from the expressions for the 
vibrational rotational terms, (1) and (5). Kighty- 
one equations of type (13) were fit by the method o} 
least squares. The resulting values of the molecular 
constants, where they can be compared, are in re- 
markable agreement with those determined pre- 
viously [5, 6]. The agreement between observed and 
calculated frequencies about 0.01 em ~!, and the un- 
certainties in the constants are very small. The 
derived constants are given in table 7. The size of 
the higher order terms indicates that there is no 
l-tvpe doubling interaction in the excited state of 
this band. Using the methane geometry it is pos- 
sible to calculate a value for the ground state 1 ot 
this molecule, A=2.628 em~'. With this value of 
A’ and the constants in table 7 one finds a value 
¢=0.89. Thompson [4] has estimated ¢; from the 
fundamental to be 0.67. Simple theory predicts the 
¢ here should be twice the € of the fundamental: 
however, it is thought that the ¢ value from the 
fundamental band is not sufficiently precise to 
justify comparisons at this time. 

from the 


TABLE 7. Molecular constants of CD3,H derived 


I-com pone nt of the hybrid hand near 2,600 em! 


(1—s)-F 2.977+.002 cm 
1'— BB’ 17-B” 0.0379+0.0018 
6 Sed 3.278+0.001 
BB’ 3.285 +0.001 
D dD; 5OLORK1O 
Dyx 1x 10-5 
Dre 1x10 
pi,.~D"~ ~~ 
(1-2 / 2592.637 +0005 
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The band center, vo, becomes 2,592.64 if one uses 
the above values of A’ and ¢. The separation of 
these two components is: given by g;;[(1+2)’—F, 
with 1=0 in this case. Thus from the component 
separation g55— 8.05 em—!, 


The authors express their appreciation to Carroll 
Dannemiller and Joseph Cameron for carrying out 
the least squares estimates on the 704 computer and 
to &. D. Tidwell and Jessie Kirkland for considerable 
heip in obtaining and reducing the data. 
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Apparatus and techniques are described for the titration of acids and bases by constant- 
current coulometry. The precision and accuracy of the method are indicated by titrations 
of benzoie acid, potassium acid phthalate, adipic acid, hydrochloric acid, and sodium car- 
bonate. Standard deviations of + 0.008 percent have been obtained which indicate that the 
reliability of the method is equal to or exceeds that of classical analytical procedures. 


1. Introduction 


In coulometric analysis, the laws of electrolysis are 
applied to analytical chemistry. The material to be 
analyzed is placed in an electrolytic cell and is neu- 
tralized, oxidized, reduced, or precipitated, as the 
case Way be, by a measured quantity of electricity. 
From this measurement, Faraday’s laws are used to 
determine the amount of chemical substance orig- 
inally present. 

While this technique has been used only for the 
past 20 vears, and particularly in the past decade, 
numerous procedures have been developed, especially 
in microchemical analvsis, where coulometry has 
definite experimental advantages over conventional 
volumetric methods [1, 2, 3]. 

With the evaluation of the faraday constant 
with ever-increasing precision, by purely physical 
methods, coulometry offers the attractive possibility 
of permitting chemical standardizations to be made 
directly with reference to fundamental electrical 
constants. In fact, a proposal has been made to the 
International Union of Pure and Applied Chemistry 
to accept the electron as a fundamental chemical 
standard [4]. 

The National Bureau of Standards has had a con- 
tinuing interest in developing analytical methods of 
the highest accuracy for the determination of the 
purity of materials suitable for chemical standards. 
Because they differ so much from conventional pro- 
cedures, the precision and accuracy of coulometric 
methods are being extensively investigated. A pre- 
liminary report on this program was presented to the 
National Academy of Sciences at its meeting in 
Washington in) November 1956 [5]. The present 
paper is a detailed account of the apparatus and tech- 
niques developed for coulometric titrations of high 
precision. It is shown that the reliability of the 
results obtained is equal to or exceeds that of the 
most careful determinations made by classical analyt- 
ical methods. 


Present address: Delco-Remy Division, General Motors Corporation, An- 


derson, Ind 


2 Figures in brackets indicate the literature references at the end of this paper. 
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2. General Considerations 


2.1. Classification of Methods 


Coulometric titrations may be classified according 
to whether the electrode process directly or indirectly 
influences the composition of the solution. Acidi- 
metric titrations in which hydrogen ion is reduced 
at the cathode are examples of the former type, 
while the precipitation of halide ion by coulometri- 
cally generated silver ion is an example of the latter 
tvpe. 

A further classification may be made on the basis 
of the type of current control. Controlled-potential 
coulometry consists in maintaining the electrode 
potential at a selected value that produces only a 
desired clectrode reaction. In this case, the current 
decreases continuously and exponentially during the 
electrolysis, and asymptotically approaches zero or 
some small constant value. The current is integrated 
during the course of the electrolysis, and is discon- 
tinued when it does not differ significantly from this 
terminal value. 

Constant-current coulometry is) based on the 
maintenance of an accurately measured current and 
the determination of the time required to complete 
the desired reaction. Integration is not required, 
but the electrolysis must be discontinued at the end 
point of the desired reaction. Knd-point determi- 
nation may be made by conventional methods such 
as potentiometric, colorimetric, amperometric, or 
thermometric techniques, depending on the partic- 
ular reaction, 


2.2. Advantages of Coulometric Titrations 


In addition to the possibility of making determi- 
nations based on physical standards, coulometric 
measurements have several obvious advantages. 
Electrochemical reactions may be considered as 
titrations with electrons, a universal “reagent.” 
This reagent may be added uniformly at any desired 
rate by methods more convenient than those re- 
quired to dispense conventional chemical substances. 
Uncertainties in composition of standard solutions 
are eliminated as are volumetric errors. The high 
precision with which electrical measurements can be 








made is a distinct advantage in coulometric methods. 
Methods for the control of electrical circuits are well- 
developed, so that coulometric titrations are espe- 
cially adaptable for automatic operation. 


2.3. Sources of Error 


Tt isof prime importance that the stoichiometry of 
the electrode process and any of its indirect reactions 
be completely defined. Any uncertainty in’ the 
eurrent-chemical reaction relationship will result in 
a direct error in the determination. In = general, 
reactions which proceed with less than 100-percent 
current efficiency are not useful. 

In many cases, the reverse of the desired reaction 
occurs at the second electrode, or the product of 
the reaction at this electrode may interfere with the 
reaction of primary interest. Accordingly, it is 
frequently necessary to isolate the second electrode 
in a separate compartment that is in electrolytic 
contact with the electrolysis chamber. In such cases, 
it is necessary to minimize losses due to diffusion or 
to bulk transfer of the solution. A further source 
of error may arise from ionic migration from the 
electrolysis chamber. The familiar technique of 
addition of a supporting electrolyte will minimize 
this error. 

Errors, both accidental and systematic, can arise 
in the determination of the end point as in conven- 
tional titrations. Errors can arise from faulty cali- 
bration of the apparatus used for the measurement 
of the current and the time. 


3. Apparatus and Procedure 
3.1. Apparatus 


The acidimetric titrations to be described were 
made by the method of constant-current coulometry. 
The cell consists of a platinum cathode and a silver 
anode, with an alkali chloride solution as the support- 
ing electrolyte. The cell reaction results in the 
reduction of one equivalent of hydrogen ion (libera- 
tion of one equivalent of hydrogen gas) at the 
cathode, and the formation of one equivalent. of 
silver chloride at the anode for each faraday of 
electricity. 

Because of the small but finite solubility of silver 
chloride in the electrolyte, a small amount of silver 
is reduced and deposited on the cathode when a 
single compartment cell is used. To overcome this 
source of error, the titration cell shown in figure 1 
was designed. The cathode chamber, C, contains 
the platinum-gauze cathode and also the glass and 
calomel electrodes to indicate the end point pH 
of the solution. An encapsulated magnet bar is 
used to stir the solution. A tube is provided to 
admit nitrogen to the cell to remove carbon dioxide. 
The compartment is sealed with a rubber stopper and 
a bubble trap to exclude air. 

The anode compartment, A, contains a_ silver 
anode of large area so that the anode reaction re- 
sults only in the precipitation of silver chloride on 
its surface, no acid being formed. It is fabricated 
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Figure 1. Coulometric titration cell. 


from a sheet of silver 7.6 by 15 by 0.065 em in thick- 
ness. This is corrugated and placed as far as possible 
from the connecting port, to assure large current 
capacity and good current distribution. Such an 
electrode has a capacity of more than 1,500 coulombs 
and can be regenerated by electrolytic reduction 
of the silver chloride. More conveniently, however, 
the electrode may be renewed by chemical solution 
of the chloride in ammonia. 

The actual compartments are 400-ml beakers of 
the electrolytic type, without lip. They are con- 
nected by a 35-mm tube, approximately 125 mm in 
length, in which are sealed sintered glass disks as 
shown in the figure. Sintered disks 1 and 2 are of 
medium porosity while 3 is of fine porosity. For 
convenience, this bridge is made from commercially 
available “sealing tubes.”” Each of the center com- 
partments has a tube attached to permit emptying 
and filling under nitrogen pressure as will be de- 
scribed. A plug of agar gel, P? (3 percent agar-agar 
in 1 N potassium chloride) is placed at the anode 
side of disk 3 to prevent liquid transfer into or 
from this compartment. 

The circuit used for the maintenance of the con- 
stant current is shown in figure 2. <A storage 
battery, B,, of 48 v serves as the current source. 
It permits a high-resistance circuit so that small 
changes in cell resistance (voltage drop across the 
cell) during the course of the electrolysis will have a 
proportionally smaller effect on the current. A 
voltage dropping resistor, #;, is immersed in an oil 
bath to aid in heat dissipation and stabilize its re- 
sistance during a measurement. The combination, 
R;R,, consists of a fixed resistance and a carbon- 
compression resistor in parallel. It has an overall 
variation of about 20 ohms and is adjusted continu- 
ously during the measurement to maintain the 
current constant. Resistance Ry is used for coarse 
adjustment of the current. 
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Figure 2. Electrical circuit diagram for coulometric titrations. 


A variable resistor, Rs, is adjusted to have, within 
narrow limits, the same resistance as the cell. The 
switch S, when in position (1), permits the current 
to flow through this “dummy” resistor for pre- 
liminary adjustment. Ammeter A is used to facili- 
tate the initial adjustment of the current. Switch 
S., when thrown to position (2), causes the current 
to flow through the cell and simultaneously triggers 
the timing device. 

The current is measured by the voltage drop 
across the standard resistor, R,, a 5- or a 10-ohm 
precision resistor, depending on the value of the 
current. These resistors are placed in a constant- 
temperature oil bath. Saturated standard cells, 
suitably thermostated, serve as the primary reference 
of voltage. The voltage drop across R, is compared 
with the voltage of the standard cell, any difference 
being measured by potentiometer P. 

During a titration, the current is maintained con- 


stant by manually adjusting it in the following 
manner. By preliminary trials, it is possible to 


adjust the value of the dummy resistor to approxi- 
mately that of the cell. Accordingly, when the 
switch is thrown to the “electrolysis” position, 
potentiometer P? can be adjusted quickly to corre- 
spond to the cell current. The potentiometer key 
is then locked in the closed position. Any change 
in the current causes the galvanometer to deflect 
from the null position, and this is immediately 
corrected by adjustment of resistor Rs. 

Initially, the time interval was determined by 
means of an electric clock, driven by a tuning-fork 
regulated constant-frequency device, capable of a 
precision of one part in a hundred thousand. To 
secure higher precision in timing, and to eliminate 
errors due to clutch slippage in the clock, a quartz 
erystal-controlled time-intervel meter (TIM) was 
used in later experiments. The meter is compared 


occasionally with NBS standard time signals to 
check its performance. It is capable of an accuracy 
of one part in a million, and has never differed from 
the time signal by this amount. 

The pH measurements, on which the end point 
determination is based, are made with conventional 
glass electrodes and a line-operated indicator. A 
suitable resistance is connected in series with the 
meter of the instrument and the voltage drop across 
this is measured (14.5 mv per pH unit). This effec- 
tively expands the scale of the instrument since this 
voltage can be readily measured to a few hundredths 
of a millivolt with an ordinary potentiometer. The 
stability of the pH meter during the period of a 
determination of the end point is sufficient to war- 
rant this procedure. 

All apparatus for making the electrical measure- 
ments, including standard cells, standard resistors, 
and potentiometers, was calibrated by the electrical 
division of this Bureau. The value for the faraday 
constant used was 96,495.6 coulombs [7]. 


3.2. General Procedure 


A 1 N solution of potassium chleride is added to 
the anode compartment of the cell. A weighed por- 
tion of solid potassum chloride is placed in the 
cathode compartment and, after flushing with nitro- 
gen, freshly boiled distilled water is addedin sufficient 
quantity to make the resulting solution 1 N. Pure 
nitrogen is then introduced through the connecting 
tube and into the solution to remove dissolved carbon 
dioxide. To insure that the electrolyte is sufficiently 
acid to facilitate this operation, a few drops of acid 
are intentionally added to bring the electrolyte to a 
pH of 5. 

After purging is complete, the added acid is 
neutralized as follows. Stopcocks S; and S: are 
opened to the atomsphere to permit the solution to 
flow into the connecting tube to a depth of about 2 
mm. <A current of 10 ma is passed until the glass 
electrode indicates a pH value of 7.0.5 Nitrogen is 
then admitted into the side tube to “blow” the con- 
tents back into the cathede compartment. This 
procedure is repeated until no further change in pH 
results on mixing and the pH remains at 7.00. 

By appropriate manipulation of the stopcocks, the 
connecting tube is then filled with neutralized 
catholyte. The stopper is momentarily removed to 
introduce the acid to be titrated. 

The current is then adjusted to a value of 200 ma 
and the solution is electrolyzed until the glass elec- 
trode indicates a value very close to the end point. 
At this rate of electrolysis, the acidity of the elec- 
trolyte in the vicinity of the end point changes more 
rapidly than the glass electrode can respond. <Ac- 
cordingly, some experience may be required to insure 
that the end point ts not overrun. Nitrogen pressure 
is then used to push the contents of the connecting 
tube into the cathode compartment for mixing. 


> No appreciable error is incurred by assuming the neutral point of the solution 
to be that of pure water. 
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The connecting tube is then filled with electrolyte 
to a depth of 2 mm, as previously described, and the 
solution is electrolyzed at a current of 10 ma to 
within about 0.5 pH unit of the end point. At this 
current, the glass electrode follows the change of pH 
very closely so that there is little danger of over- 
running the desired value. 

Electrolysis is then carried out, stepwise, for con- 
venient time intervals and the pH is recorded for 
each interval. From these data, & conventional 
differential potentiometric plot is constructed to 
determine the end point. The connecting tube is 
then flushed several times with catholyte until no 
further change in pH results. From the change in 
pH resulting from this flushing, end the titration 
curve, a small correction is made to the end-point 
time already determined. 

A typical end-point determination is shown in 
figure 3. Curve P represents the plot of the poten- 
tial of the glass electrode indicator as a function of 
the elapsed time. Curve DP is the differential plot 
used to locate the end point, identified as point (. 
The tolerance limits for the location of the end point 
that would correspond to an error of +0.002 per- 
cent in the titration are indicated. Point 4 repre- 
sents the time at which electrolysis was discontinued, 








Figure 38. Typical end-point determination. 


TABLE 1. 


Potassium Potassium Benzoic 
hydrogen hydrogen wid # 
Determination number phthalate * phthalate t 
assay assay iSsay 
( 

] 49. 976 49. 979 49. OSS 
2 YY. YXO) 99. 976 99. O91 
3 YG. GTS 99. GSI 4Y, YSY 
4 49. 970 44, 976 100. 001 
5 9Y. 974 YY. O74 YY. YS5 
ty 
% 
Mean 99. 976 99. 977 99. 990 
Standard deviation 0. 004 0. 003 0. 006 


* Timed by electric clock, 
Timed by time interval meter 


while point / indicates the final pH value for the 
solution, after rinsing was completed. The corrected 
end point is then time ('+-(time .4—time Bb). 


4. Titration of Potassium Hydrogen Phthal- 
ate 


4.1. Procedure 


The potassium hydrogen phthlate used in these 
experiments was NBS Standard Sample 84d. Small 
amounts (15 to 20 g¢) were crushed in a mortar and 
dried at 120° © for 2 to 3 hr. The samples were 
stored in a desiccator until used. 

After removal of carbon dioxide from, and neu- 
tralization of the supporting electrolyte, the rubber 
stopper was removed momentarily, and the acid 
(2 g), weighed to 0.01 mg, was dumped into the cell 
from a weighed platinum boat. Reweighing of the 
boat gave the amount of acid introduced into the cell. 

The acid was then electrolyzed at a current rate 
of 200 ma and finally at a rate of 10 ma, as described 
above. 


4.2. Results 


The results obtained for a number of titrations of 
standard sample 84d are given in columns 2 and 3 of 
table 1. The values given represent the purity 
found for the sample on the basis of an equivalent 
weight of 204.216. The electric clock was used for 
measurement of time in the case of the first: group 
while the time-interval meter was used in the second 
series. 

Although separated by a period of time of about 
1 year, and determined by slightly different) tech- 
niques, the agreement of the averages for the two 
series is excellent and shows the high precision of 
coulometric titrations. 

On the basis of very careful titrations by conven- 
tional methods, sample 84d was indicated originally 
as having an assay of 100.04 percent. Subsequently, 
Bates and Wichers [6] compared this material with 
specially purified benzoic acid by differential titra- 
tion and reported an assay of 99.987 percent (¢ 
0.002). The agreement of the values of table 1 with 
this latter value is confirmation that the coulometric 
method is one of high accuracy. 


Summary of results 


Benzoic Hydrochloric Adipic A dipic sodium 
iid | acid 4 icid ad acid 6 « carbonate ® 
iSSay issay iSSay assay iSSay 

GY. 991 49, GR2 49, 935 YY, G85 49. YOU 
OY. 993 100. OO4 YY, 937 YY. OST 49. YO] 
gY, 946 100, O14 99. G37 49. GSO ag, 915 
49. 99] 100. 003 YY, 947 99. OST g9. 910 
99, OSH 100. 002 YY. G42 99. YSU 99. 920 
99. YYS 44. GYD 99. 955 
49. 4G0 49. G51 
49. 94s 
YY, 994 49. 499 YY. G45 OY. YS4 4, 911 


0. 004 0.010 0. OO7 0). OO4 0. O07 


1 Once recrystallized from melt 
« Twice recrystallized from melt 


¢ Assay based on calculated acidity by the pressure-distillation method 
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5. Titration of Benzoic Acid 
5.1. Procedure 


The benzoic acid used was prepared at NBS for 
use as a calorimetric standard. The acid was a 
coarse, crystalline material and has been shown by 
other investigators [5] not to take up water under 
conditions prevailing in this locality. Accordingly, 
it was used without further purification or drying. 

After preparation of the cell as described pre- 
viously, a 1-g sample was weighed in a platinum. boat 
and transferred to the cathode compartment. A 
period of about 30 min was allowed before starting 
the electrolysis, to permit a substantial amount of 
the acid to dissolve. Additional acid dissolves during 
the initial reduction period. When the reaction is 
approximately 98 percent complete, the current is 
discontinued, the connecting tube is blown dry, and 
sufficient time is allowed for the solution of the 
remainder of the sample (about 20 to 30 min). The 
titration is then completed with a current of 10 ma 
as already described. A final rinsing of the cell walls 
and stopper is accomplished with 15 percent, carbon 
dioxide free ethanol, to insure solution of residual 
traces of the acid. 


5.2. Results 


The results obtained in two independent series of 
measurements with different timing techniques are 
shown in columns 4 and 5 of table 1. The values 
reported represent the assay of the acid on the basis 
of an equivalent weight of 122.112. 

By the differential method already referred to, 
Bates and Wichers [6] found the assay of this acid 
to be 99.996 percent (¢=0.001). 


6. Titration of Hydrochloric Acid 


6.1. Procedure 


A solution of hydrochloric acid was prepared by 
weight dilution of constant-boiling acid. Based on 
the pressure during the distillation of the starting 
material and the dilution data, the acid was calceu- 
lated to have a concentration of 0.11485, milli- 
equivalents per gram of solution [8]. 

After the usual preparation of the cell, 50 ¢ of the 
acid was added from a weight buret. Electrolysis 
Was started immediately and continued to about 
99.8 percent of completion. The connecting tube 
was then rinsed and the neutralization was completed 
at the 10-ma rate as described previously. In this 
case, It was not necessary to prepare an end point 
curve, due to the rapid change in pH at the end 
point. Thus, a difference in pH of 0.1 unit at the 
end point corresponded to 0.005 coulombs (0.5 see 
at the 10-ma rate). It was only necessary to attain 
the same end point obtained in the preliminary 
conditioning of the cell, a value easily determined 
within 0.02 pH units. 


6.2. Results 


The results of seven coulometric titrations of this 
acid are given in column 6 of table 1. The reported 
value is the “‘purity”’ of the acid on the basis of its 
calculated concentration from pressure-distillation 
data. Expressed in other terms, the solution was 
found to have a concentration of 0.114850 milli- 
equivalents per gram, as compared with the calcu- 
lated value of 0.11485. indicated in section 6.1. 

While the agreement with the calculated value is 
satisfactory, the precision of the determination is 
disappointing. On the basis of uncertainty of the 
end point, a reproducibility of at least 1 part in 
100,000 would be expected. Undoubtedly, some of 
the error is that inherent in the dispensing of liquid 
samples. In order to estimate the magnitude of 
such errors, a series of measurements were made in 
which 50-¢ samples of water were dispensed from. a 
weight buret into tared weighing bottles. The 
difference between the weight of water dispensed 
and the contents of the weighing bottles showed a 
variation of 3.6 parts in 100,000. If all of this were 
reflected in the uncertainty of the weight of the 
sample, the larger average deviation found in these 
measurements would be partially explained. 


7. Titration of Adipic Acid 
7.1. Procedure 


The adipic acid used was reagent grade material 
which had been further purified by filtration and 
crystallization from the melt by Delmo Enagonio of 
the Bureau. The material was estimated to have 
an assay of 99.9+- percent. The titration procedure 
used was the same as that followed in the case of 
potassium acid phthalate, and the results are reported 
merely to show the adaptability of the method to a 
variety of acids 


7.2. Results 


The results of a series of titrations to determine the 
purity of this sample of adipic acid are given in 
column 7 of table 1. The material was further re- 
crystallized and the purity was found to be higher, 
as recorded in column 8. The smaller value for the 
standard deviation indicated that the homogeneity 
of the sample was also improved by the second 
purification. 


8. Titration of Sodium Carbonate 


8.1. Procedure 


These experiments were performed to extend the 
method to the titration of bases. The sodium car- 
bonate was reagent grade material which had been 
heated in a platinum dish for several hours at a 
temperature of 350° to 400° C. The assay was not 
determined by an independent method. 

The electrolyte used was 1—M sodium sulfate, 
since chlorides would be oxidized at the anode if 
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present in the solution. The agar plug was also 
made with 1—J1/ sodium sulfate. 

The solution was preconditioned as described 
previously. After addition of the base, the platinum 
electrode was made the anode and electrolysis was 
continued until approximately 99 percent of the base 
had been neutralized. The titration was completed 
at the 10-ma rate as before. Sufficient time must be 
allowed near the end point for the removal of carbon 
dioxide. 


8.2. Results 


The results of a series of five titrations of this ma- 
terial are given in column 9 of table 1.) The precision 
obtained is comparable to that found with the acids 
determined by this method. 


9. Discussion 


This investigation shows that the coulometric 
method is capable of a high degree of precision. — In 
view of the results obtained for benzoic acid and 
potassium acid phthalate, for which the assay is 
known, there is substantial evidence that the method 
is also capable of a high degree of accuracy. 

The attainment of precision and accuracy of the 
order indicated by the data requires careful control 
over several sources of error. These may be dis- 
cussed best by a division into errors of procedure and 
of measurement. 

Errors may arise as a result of the pretreatment of 
the acid to be determined, in its introduction into the 
cell, and as a result of losses during the titration. 

Errors associated with the preparation of the 
sample are frequently overlooked. In the case of 
the phthalate, for example, the procedure for drying 
the sample is critical. It is important that only 
small quantities be crushed preliminary to drying, in 
order to release microamounts of entrapped mother 
liquor from the crystals. Furthermore, the crushing 
must be done gently, to avoid generation of appre- 
ciable heat in localized areas which might result in 
decomposition of a part of the sample. Also, the 
sample should be spread out to facilitate the drying 
process. 

To minimize errors associated with the transfer of 
the sample into the cell, it has been found advan- 
tageous to stop the flow of nitrogen during this 
process. This assures that fine spray will not deposit 
on the platinum dish to cause an error in weight and 
also prevents spray material from the 
cathode compartment. 

After the sample has been introduced into the cell, 
it can be lost only through transfer of solution into 
the anode compartment, or escape as spray through 
the nitrogen exit tube. The former is effectively 
prevented by the agar plug in the connecting tube. 
The latter is minimized by maintaining only sufficient 
flow of nitrogen during the titration to guarantee a 
positive pressure of nitrogen in the cell. 

It is very important to remove residual amounts 
of strong acids from the cell, such as would be 
introduced during a cleaning operation. This is 


losses. of 


accomplished by drawing large quantities of hot 
distilled water through the frits. The effectiveness 
of the cleaning procedure can be checked when the 
solution is brought to neutrality at the start of a 
titration. Prolonged downward drift of pH on stand- 
ing and continuing change of pH on rinsing the 
connecting tube indicate the presence of acid in the 
frits. Should this situation arise, it is necessary to 
diseard the solution and to rinse the cell further with 
hot distilled water. 

Errors associated with the electrical measurements 
are easily made much smaller than other errors, 
Standard resistors, certified to a few parts in a 
million, are readily obtained, and normal tempera- 
ture control is all that is required to assure main- 
tenance of their calibration to this accuracy.  Cali- 
brated potentiometers permit measurement of the 
voltage drop, and hence the current, to better than 
one partin ten thousand. For current measurements 
of the highest precision, the arrangement used here 
is preferred, in which the potentiometer measures 
only the difference between the voltage drop across 
the resistor and that of a saturated standard cell, 
This differential method permits meesurements with 
a precision of a few parts in a million. 

No difficulties have been experienced in maintain- 
ing the current constant to better than one part in 
a hundred thousand. The galvanometer used to 
indicate potentiometric balance had a sensitivity of 
0.3 pv/mm. Accordingly, a deflection of 8 mm from 
the null position corresponded to only luv, or 1 ppm 
with reference to variation in the current. 

It was possible to maintain belance within a few 
millimeters throughout a determination, provided 
the variable resistor used for current control had a 
smooth response. In this connection, a carbon com- 
pression resistor, spring-loaded and mounted verti- 
cally, is much superior to multiturn potentiometers 
which have a tendency to become ‘noisy’ with use. 

For accurate mexsurements, it is imperative to 
insure that the current through the standard resistor 
is actually that through the cell. While this may 
seem an obvious condition, the possibilities of extrane- 
ous connection to ground (‘‘ground loops’’), even 
though of high resistance, must be considered. 

Timing errors should not seriously affect the pre- 
cision or accuracy of coulometric titrations. For 
samples of the order of about 10 milliequivalents, 
approximately 1,000 coulombs ere needed, requiring 
about 5,000 see at the 200-ma rate of electrolysis. 
Electric clocks, driven by frequeney-controlled power 
supplies, are inexpensive timing devices which pro- 
vide measurements reliable to a few parts ina 
hundred thousand. Perhaps the largest source of 
error associated with their use is the so-called clutch 
error involved in starting and stopping, which may 
amount to a few hundredths of a second per opera- 
tion. The erystal-controlled electronic time interval 
meter used here is reliable to a part in a million. 
Both techniques require that the timing switch oper- 
ate in close synchronization with the current switch. 

Errors arising from the determination of the end 
point depend on the slope of the pH curve in this 
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region, as Well as the sensitivity of the pH indicating 
device. In the case of the weak acids used in this 
investigation (pK values, 4.5 to 5.5), the uncer- 
tainty of the end point is about 8 millicoulombs, 
corresponding to errors of less than 10 ppm. For 
strong acids the end point error is much smaller 
than this, due to the greater slope of the end point 
curve, and depends largely on the reliability of the 
pH indication. 

Errors in accuracy arise from the deviation of the 
inflection point in a potentiometric titration from 
the stoichiometric equivalence point. According to 
Roller’s equations [9], errors from this source should 
be less than 0.001 percent. There is some evidence 
to indicate that the experimental deviations are 
somewhat greater than the calculated values [6]. 
However, insufficient work has been done to justify 
the assignment of corrections to the values reported 
in this paper. 

10. Summary 


This investigation has shown that coulometric 
acidimetric titrations can be made with a precision 
and accuracy equal to or exceeding that obtainable 
by the most careful conventional chemical methods 
of analysis. Furthermore, reliable results are less 


dependent upon the manipulative skill of the analyst 
than are classical methods. The method has the 
further advantage of being absolute in the sense 
that the results are based on electrical standards 
and no assumptions as to the purity of chemical 
primary standards are involved. Titrations reliable 
to a few hundredths of a percent for both strong 
and weak acids can be readily obtained. The method 
is equally well-adapted to the titration of bases. 
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The heat of formation of titanium tetraiodide was determined relative to that of titanium 


tetrabromide by comparison of their heats of hydrolysis in dilute sulfuric acid. 


The differ- 


ence in the heats of formation may be expressed by the equation: 


Til,(¢) 


AIT(25° C) - 230.91 


By taking the heat of formation of TiBr,(e) as - 


of Tils(e) is calculated to be 385.81 


+ 2 Br, (liq) 


- (0.75 kj/mole ( 


-+4.64 kj/mole ( 


TiBr,(c) + 21,(e), 


-55.19 


0.18 keal/mole). 


616.72 + 4.60 kj/mole, the heat of formation 


91.21 +1.11 keal/mole). The 


heats of hydrolysis of TiBry and TiCl, were similarly measured; the value obtained for the 


differenee (186.77 


1. Introduction 


There are no experimental values available in the 
literature for the heat of formation of titanium tetra- 
iodide. Brewer, Bromley, Gilles, and Lofgren [1] ! 
calculated the heat of formation from estimated heats 
of solution and lattice-energy calculations. They 
obtained —101 +10 keal/mole for the heat of forma- 
tion of Tif, at 25° C. 

The direct combination of elemental titanium 
with gaseous iodine is not feasible because of the low 
vapor pressure of iodine at ordinary temperatures; 
at higher temperatures, the tetraiodide tends to de- 
compose into Til, and Til, [2]. 

The heats of formation of titanium tetrachloride 
and titanium tetrabromide have recently been deter- 
mined by the direct combination of the elements in 
a calorimeter [3, 4}. 

It is possible to hydrolyze Tily, as well as TiCl, 
and TiBr,, in dilute sulfuric acid. We can assume 
that the final state of the Ti‘? ion is the same in 
each case and that the differences in the heats of 
hydrolysis are measures of the differences in the heats 
of formation. We may, therefore, compare the heats 
of hydrolysis of Til, and TiBr, and obtain the heat of 
formation of Til, relative to that of TiBr,. We may 
also compare the heats of hydrolysis of TiBr,y and 
TiC], and obtain the difference between their heats 
of formation. This value should agree with the 
difference between the directly determined heats of 
formation. 

The present investigation is comprised of two 
parts. The first part consists of the measurement 
of the heats of hydrolysis of Til, and TiBry and the 
caleulation of the heat of formation of Til, relative 
to that of TiBr,;. The second part includes the 
measurement of the heats of hydrolysis of TiBr, 
and TiCl, and the comparison of the difference 
between these heats of hydrolysis with the difference 
between the heats of formation. 


- 1.34 kj/mole) is in good agreement with the difference between the 
directly determined heats of formation (187.11 


5.35 kj/mole). 


2. Source and Purity of Samples 


The TiCh, TiBry, and Til, were prepared by the 
Inorganic Chemistry Section of the Chemistry 
Division. The purity of the TiCl, and TiBr, was 
determined to be 99.999 and 99.998 mole percent, 
respectively, from calorimetric freezing-point meas- 
urements conducted by George T. Furukawa of the 
Thermodynamics Section of the Heat Division. The 
purity of the Til, was determined by analysis to be 
99.91 mole percent; this analysis was performed by 
the Inorganic Chemistry Section of the Chemistry 
Division. 

The 1—N sulfuric acid solution used as the ealori- 
metric fluid for the hydrolysis experiments was 
prepared from concentrated, reagent-grade acid. 

The hydrobromic and hydriodic acids were redis- 
tilled from reagent-grade solutions; the constant- 
boiling fractions were collected under nitrogen and 
stored in darkness in completely filled, glass-stoppered 
bottles. The hydrochloric acid was taken directly 
from a fresh bottle of reagent-grade acid. 

The TiCl, and TiBr, were distilled in vacuum into 
Pyrex glass bulbs, and the bulbs sealed. Soft glass 
bulbs were filled with the solid Til, in a helium 
atmosphere. The bulbs were then fitted with poly- 
ethylene caps, removed from the helium atmosphere, 
cooled in ice, and sealed. The hydrochloric, hydro- 
bromic, and hydriodic acid solutions were placed in 
soft glass bulbs by means of a funnel drawn down to 
acapillary. They were then cooled in ice and sealed. 

The quantity of each of the titanium tetrahalide 
samples was determined by weighing the bulbs before 
and after filling. The quantities of the hydrobromic 
and hydriodic acid samples were determined from the 
mass of the sample and the concentration of the stock 
solution. The quantities of the hydrochloric acid 
samples were determined by analysis of the resulting 
calorimetric solution. 


1 Figures in brackets indicate the literature references at the end of this paper. 
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3. Apparatus and Procedure 


The calorimeter used for measurement of the heats 
of hydrolysis and dilution is shown in figure 1. It 
consists of a glass vessel having a silvered, evacuated 
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Glass calorimeter. 


FIGURE 1. 


jacket, a glass-enclosed manganin heating coil, a 
glass stirrer, a glass bulb-crusher, and = a_ glass- 
enclosed platinum resistance thermometer. The 
upper part of the vessel consists of a 60/35 standard- 
taper joint; the volume of the enlarged portion is 
approximately 500 ml. 

The heating coil, stirrer, bulb-crusher, and plati- 
num thermometer are all supported by the brass 
head. The resistance of the heating coil is 124.38 
ohms. The heater leads, which pass through the 
supporting tubes, have branch points, similar to a 
4-lead thermohm, located approximately midway 
between the calorimeter and jacket) boundaries, 
The stirrer shaft is constructed of stainless steel and 
passes through a stainless steel tube fitted at both 
ends with Teflon gaskets which serve as bearings. 
The glass stirrer is sealed into a stainless steel cup 
welded to the stirrer shaft. The calorimetric 
assembly was immersed in a thermostatically-con- 
trolled water bath during each experiment. 

The apparatus for measurement of calorimeter 
temperatures, for calibration with electrical energy, 
and the details of calorimetric procedures have been 
described in previous papers [8, 4]. 

In order to compare the heats of hydrolysis. of 
TiBr, and Til,, it was necessary that the composition 
of the end solution be essentially the same in each 
case. For this reason 0.04 mole of hydriodic acid 
was added to the sulfuric acid solution for the TiBry 
hydrolysis experiments, and 0.04 mole of hyvdro- 
bromic acid was added to the solution for the Til, 
hydrolysis experiments. A’ similar procedure was 
used in the determination of the heats of hydrolysis 
of TiBr, and TiC,. 

The solutions resulting from the Til, hydrolysis 
experiments were carefully removed from the calorim- 
eter and neutralized with ammonium hydroxide 
solution. The solution was boiled to coagulate the 
precipitate which was then filtered and ignited to 
constant weight and weighed as TiQ,. The average 
ratio of the TiO, found to the theoretical quantity 
based upon the weight of sample was 1.00224 
0.0028. 

The solutions resulting from the Til, hydrolysis 
experiments and from the hvdriodic acid dilution 
experiments varied in color from faintly amber to 
brown. After standing for a few hours they became 
noticeably darker in color, presumably because of 
oxidation of the hyvdriodic acid. 

The solutions resulting from the TiBr, hydrolysis 
experiments and the hydrobromic acid dulition 
experiments were clear, but usually turned) very 
faintly amber after standing, probably because of 
oxidation of the hydriodic and hydrobromic acids. 

The solutions resulting from the TiC], hydrolysis 
experiments were very. slightly turbid, but) they 
became clear after standing for a few hours. This 
was probably caused by the rapidity of the reaction 
which, for some unknown reason, varied considerably 
among experiments. 
the reaction caused liquid to be splashed on the 
calorimeter walls above the level of the solutions, 
which made it necessary to reject the experiment. 
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In some cases the violence of 





4. Experimental Results 
4.1. Hydrolysis of Til, 


The results of the electrical calibration experiments 
are given in table 1, where / is the electrical energy, 
in joules, added to the system; A/te is the correctec 
temperature rise in ohms; and /, the energy equiva- 
lent of the “standard” calorimetric system in joules 
per ohm. 


Tasue 1. Electrical calibration of the Vily system 
Experiment E ARe EB. 
j Ohm j/johm 
| 4559. 62 0. 208827 21834. 4 
y a 4554. 26 . 208568 21833. 3 
3. 4541. 07 . 2080388 21828. 4 
Din n setde 4501. 31 . 210309 21831.3 
5 : 4541. 72 . 208048 21830. 3 
Mean 21831. 5 
Standard deviation of the mean +1. 1 


The results of the Til, hydrolysis experiments are 
given in table 2. The quantity Ae is the deviation 
in joules per ohm of the electrical energy equivalent 
of the actual calorimetric system from that of the 
calibrated system. This deviation includes the heat 
capacity of the sample and of the glass bulb, less 
that of the empty bulb used in the calibration. The 
energy in joules, g, evolved by the process is obtained 
as the product of Affe and the actual energy equiva- 
lent of the system. The heat of hydrolysis, —A// 
(27° C), is therefore the ratio of g to the number of 
moles of sample. For these calculations the heat 
capacities of Til,(¢) and Pyrex glass were taken as 
(0.235 and 0.795 j/g°C, respectively. For conversion 
to the conventional thermochemical calorie, one 
calorie has been taken as equivalent to 4.1840 J. 

The heat of hydrolysis obtained in table 2 cor- 
responds to the process: 


Til,(c) + [24H,SO,+ 4HBr-+ 2600 H,O](soln) 
[Tit*+41- +4H* + 4Br- 4+ 24H,SO,+ 2600H,O}] 
(soln), (1) 
AH (27°C) = —217.81 +0.42 kj/mole. 


TABLE 2. Results of the experiments on the hydrolysis of Tils 


4.2. Hydrolysis of TiBr, in Sulfuric Acid Solution 
Containing Hydriodic Acid 


The results of the electrical calibration and hy- 
drolysis experiments are given in tables 3 and 4, 
respectively. The heat capacity of crystalline TiBr, 

, 9F Pe fad ry. ° . 
was taken as 0.358 j/g°C [5]. The resulting heat of 
hydrolysis corresponds to the process: 


TiBr,(c) + [24H.SO,+4HI +2600 H,O](soln) 


[Tit +41-+4H+ +4Br-+24 H,SO, 
+2600 H,O](soln), (2) 
AH (27°C) = — 244.63 + 0.20 kj/mole. 
TaBLE 3. Electrical calibration of TiBrs system 
Experiment E ARc Es 
j Ohm | j/ohm 
1 2064. 70 0. 094642 21815. 9 
2 2460. 74 . 112853 21804. 8 
5 2458. 30 . 112683 21816. 1 
4 2459. OS . 112730 21813.9 
Bees. 2458. 50 - 112733 21808. 2 
Mean 21811.8 
Standard deviation of the mean +2.2 
TABLE 4. Results of the experiments on the hydrolysis of TiBrs 
Experiment Ae ARe q TiBrs —AH(27°C) 
jlohm Ohm j Mole | kj/mole 
1 11.9 | 0.091141 1989. 32 | 0. 00812468 | 244. 85 
2 14.8 . 110283 2407.12 | .00984278 | 244. 56 
3 14.0 . 111086 2424. 65 . 00992151 244. 38 
4 13.6 . 100716 2198. 31 . 00898268 244. 73 
Mean 244. 63 
Standard deviation of the mean a2 ss +0. 10 


4.3. Heats of Dilution of Hydriodic and Hydrobromic 
Acids 


The heats of dilution of the constant-boiling HI 
and HBr solutions in 1—N sulfuric acid are given in 
tables 5 and 6. Because of the small amount of 
energy evolved in these experiments, the actual 
calorimetric system was calibrated just prior to each 
experiment. The quantity /, denotes the energy 
equivalent of the actual calorimetric system. The 
heats of dilution correspond to the following re- 
actions: 


[4H1I+ 29.1 H,O](soln) +[24H,SO,+ 2600 H,O](soln) 
[4HI + 24H.SO,-+ 2629.1 H,O](soln), (3) 


Experiment Ae ARc q Til; ~ AH (27°C) AH (27°C) —11.256-+0.088 kj, 
h Oh Mol kj le = a . ‘ N 
1 24.9 0 073937 1616 00 0), 00742699 17 BS [4H Br+ 20.1 H,¢ )}(soln Pt (24 HSC ds 
2 ry fae . 082475 1802. 82 _ OOS26689 218. OS 
3 20.4 082103 1794.84. 00822223 218, 29 +9600 H.¢ )| (soln) 
4 26.5 071992 1573.60 | . 00724680 217, 14 ‘ 
5 29.9 . 100405 2194. 99 01009720 217. 39 bis ‘ pene 
6 25.0 065496 1431.51 00655530 | 218.37 =[4H Br+ 24H,SO,-+ 2620.1 H,O](soln), (4) 
Mean pee 217.81 AH (27°C) = —31.97 + 0.24 kj, 
Standard deviation of the mean +0). 21 
: oa ———— AH (25°C) = —31.54+0.25 kj. 
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Heat of dilution of aqueous hydriodic acid in 1—N 
sulfuric acid 


TABLE 5. 


Experiment Ee ARc q HI —AH (27°C 

jlohm Ohm j Mole kj/mole 

1 22100. 4 0. 001374 30. 37 0. 0109470 2. 77 

2 21975. 1 . OOOS24 18. 11 . 0064166 2. 822 

3 22082. 4 . GOO776 17.14 . 0060461 2, 835 

4 22075. 3 . 000572 12. 63 . 0044313 2 850 

5 20183. 8 . CO3105 62. 67 . 022347 2 804 

6 20097. 1 . 003661 73. 58 . 026262 ? S02 

Mean 2.814 

Standard deviation of the mean +0). O11 


TaBLeE 6. Heat of dilution of aqueous hydrobromic acid in 


1—N sulfuric acid 





Experiment FE ARe q HBr —AH (27°C 
j/ohm Ohm j Mole kj/mole 
] 0. OO5S812 12s. 44 0. 015962 8. 047 
2 22073. 3 . 003218 71.03 _ OO88925 7. USS 
3 20150. 6 O14119 284. 51 _ 036083 «(7 SSS 
4 19801. 4 012594 249. 38 . 031894 7.944 
Mean 7.903 


Standard deviation of the mean__ +0 030 


* Experiment number 3 was performed at 25° C rather than at 27°C, and the 
value was not included in calculating the mean. The value at 25° C is required 
for calculations in sec. 5.2. of this paper. The value of —A/H/ at 25° C was taken as 
7.885 kj/mole (from experiment number 3). This value is in good agreement with 
that calculated from the mean value at 27° C, using an estimated ACp correction. 


4.4. Heat Capacity of Hydrobromic and Hydriodic 
Acid Solutions 


Because of the lack of data on the heat capacities 
of aqueous solutions of hydrobromic and hydriodic 
acid, these values were determined in separate experi- 
ments. The effective heat capacity of the empty 
calorimeter was determined by measuring the heat 
capacity when filled with pure water, and subtracting 
the known heat capacity of the water from the ob- 
served energy equivalent. The calorimeter was then 
filled with the same volume of the hyvdrobromic acid 
solution, weighed, and the energy equivalent deter- 
mined. By subtracting the heat capacity of the calo- 
rimeter from the observed energy equivalent, the 
heat capacity of the hydrobromic acid solution was 
determined. The same procedure was used for deter- 
mination of the heat capacity of the hydriodic acid 
solution. 

The mean heat capacities obtained for aqueous 
hvdrobromic and hydriodic acid in the range from 
25° to 27° C were found to be: 


[HBr+5.0 H,O], Cp=2.008 j/g°C 


[H1+7.3 HO], Cp=1.845 j/g°C. 


4.5. Heats of Dilution of Sulfuric Acid 


The heats of dilution of the sulfurie acid solution 
by the water included in the aqueous halogen acids 


have been calculated from existing data [5]. The 

values obtained correspond to the following processes: 

29.1 HO (liq) +[24H,SO,+ 2600 H,O+4HI] (soln) 
(24H.SO,+-2629.1 H,O+4H I] (soln), (5) 


AH (27°C) = —0.402 +0.084 kj, 


20.1 H,O (lig) +[24H.SO,+ 2600 H,0+4H Br] (soln) 
[24H.SO,-+2620.1 H,O+4HBr] (soln) (6) 


AI (25°C) = — 0.301 +0.084 kj, 


13.5 H,O (liq) + [24H.SO,-+- 2600 H,O+4HCI] (soln) 
[24H.SO,+ 2613.5 H,O0+4H CI] (soln), (7) 


AIT (25°C) = —0.230 40.084 kj. 


4.6. Heat of Hydrolysis of TiCl, 

The results of the electrical calibration and 
TiCl, hydrolysis experiments are given in tables 7 and 
8, respectively. The value obtained for the heat of 
hydrolysis corresponds to the process: 

2600 H,O+ 4H Br] (soln) 


TiCl, (lig) +[24H,SO,-4 


(Ti#t +4Cl + 4H*++4Br- +-24H,SO, 
+ 2600 H,O](soln), (8) 
AF (25°C) = — 241.81 +0.72 kj/mole. 


TasBie 7. Electrical calibration of the TiClh, system 


Experiment E ARc | OF 

Ohm jiohm 
] 2907. 27 0. 133216 21823. 4 
2 2909. 5S 133345 21819. 5 
3 2903. 97 . 13838034 21828. 8 
4 2910. 32 . 1383377 21820. 3 
5 2903. 01 133112 21808. 5 
6 2900. 78 133003 21809. 9 
Mean cs 2ISIS. 4 
Standard deviation of the mean +3. 2 


TABLE 8. Results of the experiments on the hydrolysis of Til ‘I, 


Experiment A¢ ARc q TiCls AIT(27 
C) 

ohm Ohm j Mole kj/ mole 

] 12.5 0, OS7154 1002. 65 0. OO78529 242. 29 
2 15.7 . 118947 2487. 93 . 0102626 242. 43 
3 20.7 . 146828 3206, 59 . 0132998 241.10 
4 23.0 162914 3558. 27 .O147774 240. 79 
5 12.8 OSU692 1958. OS . OOSO773 242. 42 
Mean 241 81 
Standard deviation of the mean £0). 36 


4.7. Heat of Hydrolysis of TiBr, in Sulfuric Acid 
Solution Containing Hydrochloric Acid 


The results of the electrical calibration and TiBr, 
hydrolysis experiments are given in tables 9 and 10, 


164 





The value obtained for the heat of 


respectively. 
hvdrolvsis corresponds to the process: 


TiBr,(e) + [24H.SO,+-2600 H,0+4HCI] (soln) 
Tit* +4Br- + 4H*+4Cl- + 24H,SO, 
+2600 HO] (soln), (9) 





AH (25°C) 240.80+0.62 kj/mole. 


4.8. Heat of Dilution of Hydrochloric Acid 


The results of the hydrochloric acid dilution ex- 
periments are given in table 11. As in the case of 
the hydriodic and hydrobromic acid dilution experi- 
ments, the actual calorimetric system was calibrated 
electrically for each experiment. The quantity of 
hydrochloric acid was determined by analysis of 
the resulting solution after each experiment. The 
value obtained for the heat of dilution corresponds 
to the process: 


HCL» 13.5 H.O](soln) + [24H.SO,-+ 2600 H,O] (soln) 
24 HLSO,+ 2613.5 HO + 4HCI] (soln), 

(10) 

04.241 


AIT (25°C) + 0.320 kj. 


Panne 9. Electrical calibration of the TiBry system 
Experiment E ARc E 
j Ohm ohm 
l 2910. 17 0. 132796 21914. 6 
4 29007. 41 132761 21899. 6 
} 2902. 85 132508 21907. 0 
i 2002. 5b 132474 21910. 4 
f 2920). 41 138284 21911. 2 
\iean 219US. 6 
Standard deviation of the mean $2.5 





Parte LO. Results of the experiments on the hydrolysis of TiBr, 
Experiment Ae AR q TiBry AH (25°C 
ohm Ohm } Vole kj) mole 
13.7 0. 122050 2675. H2 0. OLLI647 239. 65 
19.4 146890 $221.00 01383545 241.19 
3 15.0 120186, 2832, 22 OLLT395 241.26 
} 11.0 OU36S1 2058. 45 OOSS12 241. 24 
10.8 O99166 2173. 66 0090327 240. 4 
Mean 240.80 
Standard deviation of the mean +0. $1 
Taste Il. Heat of dilution of hydrochloric acid in 1-N 
sulfuric acid 
Experiment OF ARc (] HCl \77(25°C 
ohm Ohm j Vole lj) mole 
22129. 1 0. 025049 554. 31 0. 041259 13. 43! 
22125. 3 024841 549. 61 040435 13. 502 
a 22129. 8 024867 550. 80 (40749 13. 505 
' 22110. 0 . 022551 198. 60 036502 13. 650 
22111.9 . 028629 522. 48 038246 13.610 
Mean 13. 560 
Standard deviation of the mean +0. 040 


812903—59 D 


5. Heats of Formation 


5.1. Titanium Tetraiodide Compared to Titanium 
Tetra bromide 


Appropriate combination of eq (1), (2), (3), (4), 
5), and (6) leads to the process: 


Til,(e¢) + [4H Br 9H.O (liq) 


TiBry(c) 4 


+ 20.1 EH.O}(soln) 4 


[4H 129.1 H.O](soln), (11) 


AIT(27°C) = + 6.00 £0.54 kj. 

The heat capacities of the aqueous hydrobromic 
and hyvdriodic acid solutions, for the range 25° to 
27°C, have been taken as 2.008 and 1.845 j/g¢°C, 
respectively, from section 4.4 of this paper. The 
heat capacities of TiBr, and Til, have been taken as 
0.358 and 0.235 j/g°C, respectively. By means of 
these data we obtain for eq (11): 


AH (25°C) = +-6.29 £0.55 kj. 

The heat of formation of HL in [H1+7.3 H,O}(soln) 
was taken as —50.66+ 0.04 kj/mole [5], and of HBr 
in [HBr+5.0 H,O](soln) as —109.96+0.18 kj/mole 
[5]. 
Thus, we obtain for the difference in heats of for- 
maticn of Til, and TiBr;: 

21.(e), (12) 


Til,(¢) +2Br, (liq) = TiBr,(e) 4 


AH (25°C) 230.91 +.0.75 kj/mole 


55.19+0.18 keal/mole. 


Since the heat of formation of TiBr,(c) has been 
previously measured [4] by the direct reaction of 
titanium with bromine [A/7/f{(25°C) = —616.72 44.60 
kj/mole], we obtain the following value for the heat 
of formation of Til,: 

Til,(e) (13) 


Ti(solid) +21 (e) 


AEFFP (25°C) 385.81 +4.64 kj/mole, 


92.21 +1.11 keal/mole. 


5.2. Titanium Tetrachloride Compared to Titanium 
Tetra bromide 


The appropriate combination of eq (6), (7), (8), 
(9), (10), and (4) leads to the process: 


TiBr,(e) + [4HCl+ 13.5 H,O](soln) + 6.6 H,O (liq) 


TiC], (liq) + 4HBr+20.1 H,O]}(soln), (14) 


AH(25°C) 21.764 1.04 kj/mole. 
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The heat of formation of HCL in [HCl+3.4 H,O} 
(soln) is taken as —151.21+0.17 kj/mole [5] and of 
HBr in [HBr+-5.0 H,O](soln) is taken as —109.96 
+0.13 kj/mole [5]. Thus, we obtain for the differ- 
ence in heats of formatien of TiBr, and TiCl,: 


TiBr,(c) +2CL(g) = TiC, (liq) + 2Br, (liq), (15 
AH (25°C) = — 186.76 + 1.34 kj/mole, 
= —44.64+0.32 keal/mole. 
This value is in very good agreement with —44.72 


+1.28 keal for the difference in the directly deter- 
mined heats of formation of TiBr, (A//f= — 147.40 
+1.10) and TiCl, (AHf=—192.12+0.65) [8, 4]. 
This gives added confirmation to the values obtained 
for the heats of formation of TiCl, and of TiBry. 

Gross, Hayman, and Levi [6, 7] have also reported 
direct determinations ef the heats of formation of 
TiBr,(AHf=—148.10+0.25 keal/mole) and TiCl, 
(AHf=—191.45+0.30 keal/mole). The difference 
between these values is —43.35+0.35 keal/mole. 
which is in reasonably good agreement with the 
value obtained in the present investigation. 
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Reactions of Pentafluorohalobenzenes *’ 


Walter J. Pummer and Leo A. Wall 


(June 15, 1959) 


Both pentafluorobromobenzene and pentafluoroiodobenzene are reactive intermediates 
and can be used to introduce the perfluorophenyl (CsF;—) group into a variety of new 


compounds. 


The chief methods used are the Grignard coupling or addition reaction as 


well as the Ullmann condensation reaction. Some of the new compounds prepared were pen- 
tafluorobenzonitrile, perfluorodiphenyl], pentafluorophenyl-a-ethanol, and pentafluorostyrene. 


1. Introduction 


Pentafluorobromobenzene is obtained as one of the 
products from the synthesis of hexafluorobenzene by 
the pyrolysis of tribromofluoromethane [1]... The 
isolation of this material is complicated by the 
presence of various fluorobromo olefinic and satu- 
rated aliphatic compounds boiling in the same range. 
These undesirable components can be effectively 
eliminated by subjecting the residue (bp >110° C) 
from the initial pyrolysis to another pyrolysis. 


! This paper is based on work sponsored by the Bureau of Aeronautics and the 
Oflice of Naval Research, Department of the Navy, Washington, D.C. 

? Presented as part of the Fluorine Symposium, Division of Industrial and 
Engineering Chemistry, 134th meeting of the American Chemical Society. 
Chicago, Il., Sept. 1958 

Figures in the brackets indicate the literature references at the end of this 


ice the preparation of the article, work similar in part, has been published 
. Nield, R. Stephens, and J. C. Tatlow, J. Chem. Soc. 1959, 166. 











This procedure converts the contaminants to more 
hexafluorobenzene and to higher boiling analogs 
without any apparent destruction of the penta- 
fluorobromobenzene. A relatively pure sample can 
then be obtained by simple fractionation. Penta- 
fluorobromobenzene and pentafluoroiodobenzene can 
also be prepared by the bromination or iodination 
of pentafluorobenzene in 65 percent oleum.* 
Pentafluorobromobenzene (1) and pentafluoroiodo- 
benzene (IIT) are useful reactive intermediates 
Since they contain labile bromine or iodine atoms, 
these compounds can be used to synthesize numerous 
new aromatic fluorocarbons containing the per- 
fluorophenyl grouping, C;F; The methods chiefly 
used are the Grignard coupling or addition reaction 
and the Ullmann-type condensation. The following 
scheme will serve to illustrate some of these reactions: 
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Perfluorodiphenyl (V) was obtained in good yields 
by the Ullmann condensation of pentafluorobromo- 
benzene (I) in a sealed Pyrex tube at 180° C. The 
Grignard coupling method with cobalt chloride gave 
a low yield. 

The preparation of perfluorobenzonitrile (IV) 
directly from pentafluorobromobenzene (1) following 
the Braun method [2], the use of pyridine and cuprous 
cyanide, gave unsatisfactory results. The reaction 
is extremely sensitive to heat, and after a short time 
at 138° C black insoluble solids appeared and no 
identifiable products could be obtained. The nitrile 
(IV) is also obtained from the reaction of (III) and 
cuprous cyanide. In this case the iodine is even 
more readily removed than the bromine, and once 
the reaction was initiated, the products were re- 
moved immediately by vacuum distillation. This 
procedure led to fair yields of pentafluorobenzo- 
nitrile (IV), based on its hydrolysis product, penta- 
fluorobenzoic acid (VITT). 

Pentafluorobromobenzene (1) readily forms a 
Grignard reagent, pentafluorophenylmagnesium bro- 
mide (II). Acidification of this reagent vields 
pentafluorobenzene (VI). The addition of acetalde- 
hyde to (II) forms pentafluorophenyl-a-cthanol 
(VII), which when dehydrated over alumina pellets 
at 350° C gave pentafluorostvrene (IX). Tarrant 
[3] emploved the coupling of aromatic Grignard 
reagents and fluorochloro olefins to produce various 
stvrenes halogenated in the side chain. Several 
attempts were made to utilize this concept for the 
synthesis of perfluorostyrene via the coupling of (11) 
and trifluoroiodoethene [4] or its precursor, 
1,1,2-trifluoro-1-chloro-2-iodoethane. Instead of the 
desired coupling effect, an exchange reaction occurred, 
and pentafluoroiodobenzene (IIT) was obtained. 

A similar attempt at preparing perfluorostyrene 
involved a mixed Ullmann reaction between penta- 
fluorobromobenzene (1) and trifluorovinyl bromide. 
Only small amounts of perfluorodiphenyl (V) were 
isolated. No mixed coupling products were ob- 
served. Reactions involving the Grignard reagent 
(II) and tetrafluoroethyvlene are at present being 
explored. 

Initial attempts to carbonate the reagent (IT) 
by pouring it into an ether-solidified carbon dioxide 
slurry failed to produce any pentafluorobenzoic acid 
(VIII) after acidification. Only pentafluorobenzene 
(VI) was obtained. In this respect the Grignard 
reagent resembles — pentachlorophenylmagnesium 
chloride in that it also failed to carbonate [5]. 


2. Experimental Procedures * 


2.1. Pentafluorobenzene (VI) 


To 1 g (0.041 g-atom) of magnesium turnings in 
10 ml of anhydrous ether was added 1 ml of penta- 
fluorobromobenzene (I) and a small crystal of 
iodine. Local heating was necessary to initiate the 
reaction. When the reaction was sufficiently in 


The boiling points given here are not corrected. Chemical analyses were 


performed by E. Deardorff of the Analytical Section, and the mass spectrometer 


analyses were obtained by Miss E, Quinn of the Mass Spectrometer Section 
of NBs. 


progress the flask was cooled in an ice-water bath 
and the remainder of compound (I) (total 10 ¢ 


. =) 
(0.04 mole) was added over a half-hour period. The 
Grignard solution is dark brown in color. When 


all the magnesium appeared to have been consumed, 
the solution was allowed to come to room temper- 
ature and was stirred for an additional one-half 
hour. The Grignard solution was then poured on to 
50 g of solidified carbon dioxide. After 1 hr, 100 ml 
of 10 percent hydrochloric acid was added. No 
pentafluorobenzoic acid (VIII) was observed, but 
there was obtained 3.5 g (51%) of pentafluoro- 
benzene, bp 85° to 86° Cynp'**=1.3925; reported 
bp 88° to 89°; np*=1.3931 [6]. 


2.2. Pentafluoroiodobenzene (III) 


To 200 ¢ (1.5 moles) of iodine and 1 kg of 65 
percent oleum, after stirring for 1 hr at room tem- 
perature, was added 255 g of reduced hexafluoro- 
benzene (a mixture containing 45% of hexafluoro 
benzene, 40% of pentafluorobenzene, and 10% of 
tetrafluorobenzene) [7|.. The mixture was. stirred 
for 4 hr at 55° to 60° C. It was allowed to come 
to room temperature overnight. The flask was 
cooled in an ice bath, and gradually 1 liter of ice 
water was added. It was diluted further with 
another liter of ice water, then decolorized with 
aqueous sodium bisulfite, and the products (177 ¢) 
separated. After drying (Na SO,) and distillation, 
there was obtained 82 ¢ of unreacted hexafluoro- 
benzene, 36 g of pentafluoroiodobenzene (bp 73° to 
mm, Np’ =1.4990), and 48 ¢@ of higher 
iodinated fluorobenzenes. Analysis—calculated for 
C.F;1: C, 24.5; 1, 43.1. Found: C, 24.7; 1, 42.1. 


2.3. Perfluorodiphenyl (V)—-Ullmann Method 


Five grams (0.02) mole) of pentafluorobromo- 
benzene (1) and 2.5 ¢ (0.039 g-atom) of activated 
copper powder [8] were sealed in’ an evacuated 
5-mm Pyrex tube 20-cm long. The tube was placed 
in a furnace at 180° to 190° C for 48 hr. At the 
end of this period the temperature was slowly 
raised to 290° and held at this temperature for 6 hr. 
It was allowed to cool. Upon cooling, the contents 
solidified in long white needles. The tube was 
cooled in liquid nitrogen and opened. There was 
obtained 2.95 g (87%) of perfluorodiphenyl, mp 
63° to 66° CC. A reerystallization from methanol, 
followed by sublimation at 50°/1 mm afforded 
white plates, mp 67.5° to 68° C. Analysis 
calculated for CyoFyy: C, 43.05. Found: C, 42.9. 


2.4. Perfluorodiphenyl (V)-Grignard Coupling 
Method 


Sena 
19° /30 


The Grignard reagent, prepared from 10 ¢ (0.04 
mole) of compound (I), 1 g (0.04 g-atom) of magne- 
sium turnings, and 25 ml of anhydrous ether, was 
filtered into a dropping funnel and added slowly to a 
refluxing mixture of 0.4 ¢ (0.003 mole) of anhydrous 
cobalt chloride and 3 ml of ethyl bromide in 3 ml of 
anhydrous ether. The addition required 45 min. 
The mixture was then refluxed for 2 hr more, cooled, 
and acidified with 50 ml of 10 percent hydrochloric 
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acid. The ether layer was separated, dried (Na SQO,), 
and distilled. There was obtained 2.5 g (38%) of 
pentafluorobenzene, bp 84° to 86° C, and 0.5 ¢ of 


on] 
& O07 


perfluorodiphenyl (V) (15%), mp 67° to 68° C. 


2.5. Pentafluorophenyl-a-Ethanol Alcohol (VII) 


To the Grignard reagent, as prepared from 10 ¢ of 
compound (I), was added rapidly at 0° C, 6.6 g (0.15 
mole) of acetaldehyde. The solution was stirred at 
0° C for 1 hr more, and finally decomposed by the 
addition of 50 ml of 6N hydrochloric acid. The 
ether laver was separated, dried (Na :SO,), and 
distilled. There was obtained 7 g (81%) of penta- 
fluorophenvl-a-ethanol (VIT), bp 80° to 82°/37 mm; 
n?— 1.4426. Analvsis-calculated for CsH;O;F: C, 
45.3; H, 2.35. Found: C, 44.0; H, 2.30. 


2.6. Pentafluorostyrene (IX) 


Two grams of pentafluorophenyl-a-ethanol (VIT) 
was passed through a 10-mm diameter by 20-em 
glass tube containing 5 ¢ of alumina pellets [9] at 
345° to 350° C with the aid of nitrogen gas. The 
product (1.4 g) was caught in a solidified carbon 
dioxide trap and later distilled after drving (Na SOx). 
There was obtained 0.6 ¢ (83%) of pentafluoro- 
stvrene (IX), bp 140° to 141° C, n'$=1.4414. Mass 
spectrometer analysis of the product indicated a 
slight contamination with the starting material (VIT). 


2.7. Attempted Preparation of Perfluorostyrene 

(a) The Grignard reagent (IT) was prepared from 
10 ¢ (0.04 mole) of pentafluorobromobenzene (1), as 
described previously. To the refluxing Grignard 
solution was added 9.76 ¢ (0.04 mole) of 1,1,2-tri- 
fluoro-1-chloro-2-iodoethane in 5 mil of anlwdrous 
ether. No noticeable heat effect was observed. 
The solution was refluxed for 3 hr. At the end of 
this time a black solid appeared on the sides of the 
flask. The mixture was allowed to stand at room 
temperature overnight, and finally it was decom- 
posed by addition of 50 ml of 6N hydrochloric acid. 
The ether laver was separated, dried (Na SO,), and 
distilled. There was obtained 5.78 g (50%) of 
pentafluoroiodobenzene (IIT), bp 77° to 78°/35 mm, 
n2°— 1.4990, 

Mass spectrometer analysis of this product showed 
it to be 98 percent pure. Also, a carbonlike solid, 
which has not been characterized, was obtained. 

(b) The Grignard reagent (IT) was prepered from 
10 ¢ (0.04 mole) of pentafluorobromobenzene (1) in 
the usual fashion. To the Grignard solution, main- 
tained at 0° C by means of an ice bath, was added 
12.4 ¢ (0.06 mole) of trifluoroiodoethene. The solu- 
tion was stirred for 1 hr more and then refluxed for 
2 hr. Decomposition was effected in the usual 
fashion. Distillation of the residue, after removal of 
the ether, afforded 2 ¢ of pentafluorobenzene (VT), 
1 ¢ of pentafluoroiodobenzene (IIT), and 4.5 ¢ of 
unreacted pentafluorobromobenzene (1). 

(c) Two and one-half grams (0.01 mole) of com- 
pound (1), 3.71 g (0.06 mole) of activated copper, 
and 3.22 ¢ (0.02 mole) of trifluorovinyl bromide were 
sealed in an evacuated 5-mm tube 25-em long. The 


tube was placed in a furnace at 100° C, and the 
temperature was gradually raised to 200° C. This 
temperature was maintained for 5 hr. No apparent 
reaction was observed. The temperature was raised 
slowly to 260° to 270° C, and the tube was kept at 
this temperature for 15 hr. The copper became 
straw colored. The tube was cooled in liquid 
nitrogen and opened. The mixed Ullmann reaction 
produced only 0.1 ¢ (2%) of perfluorodiphenyl (V). 
No mixed coupling products were observed. 


2.8. Perfluorobenzonitrile (IV) 


(a) From pentafluorobromobenzene (I), 5 g (0.02 
mole) of compound (1), 3.58 g (0.02 mole) of cuprous 
evanide, and 1.6 g (0.02 mole) of pyridine were 
refluxed together for 1 hr. The blackish mixture 
was cooled, and 50 ml of 6N hydrochloric acid was 
added. There was obtained a carbonlike insoluble 
solid, along with 30 percent (1.5 g) of recovered 
pentafluorobromobenzene (1). No pentafluoroben- 
zonitrile (IV) was isolated. 

(b) From pentafluoroiodobenzene (IIT), 5 g (0.017 
mole) of compound (II1), 1.6 g¢ (0.009 mole) of 
cuprous cyanide, and 1.34 g (0.017 mole) of pyridine 
were placed in a small flask and gradually heated. 
The contents began to darken with increasing 
temperature. At 150° C the contents were a black 
viscous liquid. This temperature was maintained 
for 5 min more. The flask was allowed to cool to 
100° C, and the apparatus was rearranged for distilla- 
tion. The products were removed under reduced 
pressure obtained by a water aspirator. There was 
obtained 3.07 g of vellow liquid, bp 185° to 190° 
C, np® 1.4764. The liquid, on standing, gradu- 
ally turned brown and probably contained some 
unreacted starting material (III). To show that 
the pentafluorobenzonitrile (IV) had indeed formed, 
half of the above liquid, 1.5 g, was hydrolyzed with 
75 percent sulfuric acid at 180° C. Upon pouring 
on ice, followed by an ether extraction and removal 
of the solvent, there was obtained 0.2 g (16%, 
based on reacted material) of pentafluorobenzoic acid 
(VIII), mp 101° to 108° C (reported mp 104° to 
105°) [10]. There was obtained 0.4 ¢ of unreacted 
starting material. 
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Tritium-Labeled Compounds I. 


Radioassay of Tritium- 


Labeled Compounds in “Infinitely Thick” Films With’a 
Windowless, Gas-Flow, Proportional Counter ' 


Horace S. Isbell, Harriet L. Frush, and Ruth A. Peterson 


(May 5, 1959) 


A simple, sensitive, and reliable technique has been devised for the radioassay of non- 
volatile, water-soluble tritium compounds. The substance to be analyzed is dissolved in 
an aqueous solution of a thickening agent, preferably sodium O-(carboxymethy]) cellulose or 
sodium alginate. The solution is placed in a shallow planchet, and after evaporation of the 
water, the resulting film, which is “infinitely thick” to tritium beta particles, is counted with 
a 27, windowless, gas-flow, proportional counter. By means of an empirical factor, deter- 
mined with a substance of known radioactivity, the counts are converted to microcuries. 
In a film having a thickness of 0.7 mg/cm2, the counting efficiency is about 4 percent; the 
standard deviation from the mean, obtained in a series of routine Measurements, was less 


than 2 percent. 
0.01 microcurie of radioactivity. 


An assay can readily be made with tritium-containing material having 
The method, which is applicable to nonvolatile, water- 


soluble solids, solutions, or liquids, is suitable for routine analyses. 


1. Discussion 


Tritium is one of the cheapest and most versatile 
radioisotopes for use as a tracer in chemical reactions. 
However, its widespread use has been hampered by 
a lack of convenient methods of analysis. In con- 
nection with the development of methods for the 
preparation of position-labeled, radioactive carbo- 
hydrates containing tritium, a simple, rapid tech- 
nique was needed for the assay of radioactivity of 
nonvolatile, trittum-labeled materials. Liquid sein- 
tillation counters [1],2 although very satisfactory for 
assaving dissolved samples, are both expensive and 
unavailable to many workers. The conversion of 
materials to gases, for counting in ionization cham- 
bers [1] is laborious for routine analyses. Several 
workers have counted tritium-labeled materials from 
“infinitely thick” films with proportional counters 
(2 to 7]. A tritium-containing solid, if available in 
relatively large quantity, may be packed into a 
cupped planchet and counted [2, 3, 4]; this low-effi- 
ciency method has chiefly been used for comparisons, 
rather than for the determination of absolute activ- 
ities. Alternatively, an infinitely thick film of a 
nonvolatile liquid has been employed [5], or a solid 
film has been deposited on a planchet by evaporation 
of the solvent from a solution [6, 7]. The last- 
mentioned method is probably the most versatile, 
but it suffers from the fact that crystallization on the 
planchet can render the surface area of the film un- 
certain and variable [5, 7], and may even cause con- 
tamination of the counter during flushing [6]. 

For routine radioassay of nonvolatile, water- 
soluble, tritiated materials in a windowless, gas-flow, 
proportional counter, a technique has been developed 
in this laboratory that avoids difficulty caused by 


! Part of a project on the development of methods for the production of radio- 
active carbohydrates, sponsored by the Division of Research of the Atomic 
Energy Commission. 

* Figures in brackets indicate the literature references at the end of this paper. 


crystallization and permits the counting, from uni- 
form, infinitely thick films, of materials available in 
small quantity. The radioactive substance is dis- 
soived in an aqueous solution of a thickening agent, 
preferably sodium O-(carboxymethyl)cellulose or 
sodium alginate.’ The amount of sample and thicken- 
ing solution are so adjusted that the thickness of the 
film formed on the planchet after evaporation of the 
water is more than the maximum range of the tritium 
beta particles, (infinitely thick, approximately 0.7 
mg/cm? for the material used [9]). The films adhere 
well to the planchets, and the thickener prevents 
crystallization. The counts are converted to micro- 
curies (uc) by means of an empirical factor determined 
under the same conditions with a sample of known 
‘adioactivity.!. Although the counting efficiency of a 
film having a thickness of 0.7 mg/cm? is only about 4 
percent, the reproducibility of the method cor- 
responds to a standard deviation of less than 2 
percent. A satisfactory determination can be made 
with as little as 0.01 we of a tritium-labeled material 
on a 2-in. planchet. Because of its simplicity the 
method is suitable for routine analyses. The ac- 
curacy of the method depends on the absolute activity 
of the material used as standard. For the work 
reported here, p-glucose-/-t was standardized by 
comparing it with the NBS tritium oxide standard. 

In the course of the work, a study was made with 
films of less than infinite thickness. Although high 
counting efficiencies can be obtained, uncertainty 
arising from uneven distribution of the sample makes 
counting from thin films less satisfactory than count- 
ing from infinitely thick films by the method 
described. 

3 A critical study of the energy loss and penetration range of beta particles is 
given in [8]. In general, the range increases with the atomic number of the mate- 
rial penetrated. In the film-counting method described here, the films are of 
organie materials, and have approximately the same elemental composition. 
In the presence of elements of higher atomic number, the greater depth of pene- 
tration would require a thicker film for infinite thickness, and use of a different 
factor for converting counts to microcuries. 


4 The counting of carbon-14 from thin films of agar has been reported by Mc- 
Cready [10]. See also a modification of the method [11]. 
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2. Materials and Apparatus 
2.1. Tritium-Labeled Materials 


p-Glucose-1-t was prepared by sodium amalgam 
reduction of p-glucono-6-lactone in tritiated water 
(12| in the presence of sodium acid oxalate [13]. 
The sample used as a tritium standard was assayed 
by comparing its activity in a phosphoric acid- 
phosphoric anhydride solution with that of the NBS 
standard tritium oxide sample No. 4926; the measure- 
ments were made with a proportional counter [14]. 
p-Mannitol-1-f was prepared by reducing p-mannose 
with tritiated lithium borohydride [14]. 


2.2. Solution of Sodium O-(Carboxymethy]l)cellulose 


The thickening agent was commercial, medium- 
viscosity sodium O-(carboxymethyl)cellulose, CMC- 
12MP, a product of Hercules Powder Company, 
Wilmington, Del. The CMC stock solution ordi- 
narily used for the radioassay of solid samples was 
prepared by dissolving 1 g of CMC, 0.5 g of anhy- 
drous p-glucose, and 10 mg of eosin® in sufficient 
water to give 100 ml of solution. The p-glucose 
serves as an organic “ballast”? material and plasti- 
cizer. It was omitted entirely, or in part, in thicken- 
ing solutions for use with large samples of weakly 
radioactive material. p-Glucitol has been used with 
equal success as ballast in place of p-glucose. 


2.3. Solution of Sodium Alginate 


The sodium alginate used was_ high-viscosity 
Algin, obtained from Kelco Company, San Diego, 
Calif. The stock solution was prepared from 0.5 ¢ 
of Algin, 1 g of anhydrous p-glucose, 10 mg of eosin, 
and sufficient water to give 100 ml of solution. 


2.4. Counting Equipment 


Radioactivity measurements were made with a 
27, windowless, gas-flow, proportional counter, 
Model PC—3, obtained from Nuclear Measurements 
Corp., Indianapolis, Ind.; the slide holder of the 
instrument had been modified by inserting tubes for 
water circulation [15]. Commercial gas, consisting 
of 90 percent of argon and 10 percent of methane, 
was used in the counter. The gas was dried by suc- 
cessive passage through soda-lime and anhydrous 
calcium sulfate. The 2-in., stainless-steel, cupped, 
flat-bottomed, counting planchets used in this work 
were also obtained from Nuclear Measurements 
Corp.; they have an effective area of 21.4 em?. Some 
of these planchets were modified by enclosing a 
10-cm? area with a circular groove beyond which the 
solution did not spread. 


5 The dye serves to show the uniformity of the film. 
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3. Recommended Procedure for Radioassay 
of Tritium-Labeled Compounds in Infinitely 
Thick Films 


The solution of the radioactive sample in the 
aqueous thickening agent can be prepared by any 
convenient method. Ordinarily, radioactive solids 
are weighed, and then dissolved in a known volume 
of the stock thickening solution. <A radioactive 
liquid, usually less than 0.5 ml, is mixed with a 
measured amount of the thickening solution. How- 
ever, if the liquid is sufficiently radioactive, it is 
measured in a dilution pipet and diluted to a definite 
volume with the thickening solution. Approxi- 
mately 1 ml of the solution to be assayed is trans- 
ferred to a clean,® 2-in. cupped planchet (or 500 
microliters is transferred to the 10-em? area of the 
planchet), which is then placed under an infrared 
lamp. After the solvent has evaporated, the planchet 
is kept for about 1 hr in a desiccator over a saturated 
solution of potassium acetate in contact with the 
solid phase. This rather unusual procedure gives 
films of reproducible moisture content, and does not 
lower the count substantially. The sample is then 
counted in the proportional counter for 3 periods, 
with a total of at least 10,000 counts. Before each 
period the counter is flushed with the counting gas 
for 30 see. It occasionally happens that successive 
counts show a marked trend, presumably caused by 
the accumulation of a static charge on the film. If 
this occurs, the sample is removed to the desiccator 
for a short time, and then recounted. The back- 
ground count is determined with a_ film-covered 
planchet placed in the counting chamber. 

The voltage for optimum precision must be deter- 
mined under the conditions used. The optimum 
voltage for assay of solid samples with the equip- 
ment used in the present study was approximately 
2,050 volts. 


4. Discussion of the Procedure and Results 


The activity of tritium in the sample is given by 
the relationship: 


dps 
| ar 


‘Ds x? (1 
7." 


ue 


where cps is the net count, m is the total solids in a 
volume of solution containing the radioactive sample, 
and & is a constant determined by the same method 
for a compound of known tritium content. 

The value of m must be determined accurately, 
although the exact amount of solution placed on the 
planchet is not critical if the film formed from it is 
infinitely thick. Because of hydration, the value 
of m calculated from the ingredients in the formula- 
tion may differ slightly from that determined by 
weighing the film produced, under the conditions of 


6 Care must be used to obtain a grease-free surface on the planchet. In this 
laboratory the planchets are cleaned in a solution of sodium hydroxide or 
“Versene’’ immediately before use. 
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the method, from a definite volume of the solution. 
Determination of the weight of the film in each 
analysis is laborious. This step can be avoided by 
determining the weight of the film formed on a 
planchet from exactly 1 ml of the stock solution of 
CMC or Algin under the conditions of the analysis. 
Because of the viscosity of the solution, a pipet, 
calibrated to contain 1 ml, is emploved, and the 
solution adhering to the pipet is carefully rinsed out 
and added to the planchet. The film, obtained after 
evaporation of the water, is conditioned in the usual 
manner and weighed. Several determinations are 
made, and the average value is used for calculating 
m, the weight of the film containing the radioactive 
substance, on the assumption that m equals the 
weight of the film from the amount of stock solution 
used plus the weight of the radioactive sample. All 
nonvolatile matter in the formulation must be 
included in the calculation of m. 

The value of & depends on the surface area, back- 
scattering, efficiency of the counter, geometry, and 
other factors; it has the dimension of weight7!. To 
obtain accurate results in the analysis of a series 
of compounds, & must be determined by assay of a 
standard tritium compound under the conditions 
used for the unknowns. With materials of similar 
character, the amount of radiation seattering is 
constant. Because scattering varies with the et- 
fective atomic number of the reflecting material, 
heavy clements such as barium increase internal 
scattering and may lead to small errors. These 
can be avoided by adding, to the standard, com- 
pounds of any heavy elements that may be present 
in the material to be assayed. Ions of certain 
metals, es, for example, lead, form inscluble com- 
pounds with both CMC and Algin. This undesirable 
effect can be overcome by adding Versene to the 
film-forming mixture. 

Table 1 presents results obtained with two samples 
each of p-glucose-/-¢ and pb-mannitol-/-f. The stand- 
ard deviations ere based on the deviations in counting 


rates for the separate planchets from the mean 
counting rate for the group of planchets. Sample 
G,; was used as a standard, and its specific activity 
(0.0865 ye/mg) was carefully determined by the 
comparative method described under materials and 
apparatus. Use of this material as a control gave 
values of k under various conditions, i.e., for CMC 
and <Algin films and for the 10-cm? and 21.4-em? 
(2-in.) planchets. By use of these constants, the 
specific activities of G., Mi, and M, were determined. 
The counting rates for films having thicknesses 
between 0.69 and 1.48 mg/cm? showed no significant 
differences. 


5. Assay of Tritium Compounds in Films of 
Less Than Infinite Thickness 


Table 2 and figure 1 show the relative counting 
rate of a preparation of p-glucose-t as a function of 
source thickness. Each sample was spread on a 
freshly cleaned, 2-in., stainless-steel planchet, and 
the solvent was evaporated under an infrared lamp. 
The planchet was cooled in a desiccator over anhy- 
drous calcium chloride, and then counted with the 
equipment described before. The first seven deter- 
minations were made on films from solutions of 
p-glucose-t in water, and the last three on films from 
solutions of p-glucose-t containing sufficient Algin to 
give the weight recorded. The radioactivity in thin 
films can be estimated by use of the data of table 2 
and figure 1. 


6. Discussion of Absolute Radioactivity 


The counting rate of a trittum source in an efficient, 


windowless, gas-flow, proportional counter can be 
represented by an equation of the type [6]: 


TABLE 1. Counting rates of tritium compounds in “infinitely thick” films 


Sample Film Counting Data 
| Specific 
activity 4 kX 105 
of sample 
Substance 4 Weight > Thicken- Area Total solids} Thickness No. of Average Standard | 
(m’) ing agent (m) planchets count deviation ¢ 
mg cm mg mg/cm? cps % pe/m’ 
(iy 0. 435 CMC 10.0 y He 0.77 5 50.1 1. 04 (0. 0856) 9. 65 
Gy 87 CMC 21.4 15.4 42 5 108.7 0. 86 (. 0856) 4.45 
Gy ALS Algin 10.0 8.1 Sl 1 15.1 . 58 | (. 0856) 9. 72 
Gy 83 Algin 21.4 16.2 76 ) 95.3 . 87 (. 0856) 4. 60 
Gy . 82 CMC 21.4 14.8 69 5 101.3 . 44 . 0814 | (4. 45) 
Gy 1. 23 CMC 21.4 22.1 1.03 5 100. 5 63 | . 0804 | (4. 45) 
(iy 0), 82 CMC 10.0 14.8 1.48 11 17.6 1.15 . 0829 | (9. 65) 
| 
(ig . 82 Algin 21.4 16.7 0.73 5 95.0 1. 57 . 0837 (4. 60) 
(ig 1.23 Algin 21.4 23.5 1.10 3 94.5 a 8 0831 | (4. 60) 
M}-_.- 0. 036 CMC 21.4 14.0 0. 65 4 289. 4 0.14 5.01 | (4. 45) 
My 2. 50 CMC 21.4 16.7 ys. 5 34.1 1. 40 0. 010 (4. 45) 


» G, and G2 are different preparations of p-glucose-/-t; My) and Mo are bD-mannitol-/-t. 
» m’ is the weight of the radioactive sample contained in the film of weight m. 


2 


> 
< 100 
V n—1 


1 ye=cpsXmXk, Values in parentheses were used as standards. 
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dps= 


where dps is the activity of a source, and /,, fo, and 
fs are correction factors for the geometry of the 
detector, for backscattering, and for self-absorption, 
respectively. The geometry factor, f,, with a flat 
planchet is approximately 0.5. The backscattering 
factor, f>, is a function of the atomic weight of the 
backing material and depends somewhat on the 
character of the radiation and on the geometry of 
the counter [16]. But with an infinitely thick 
source, scattering from the backing material may 


Counting rate of D-glucose-t as a function of film 
thickness 


TABLE 2. 








Volume } | Counting 

D-Glu- of film- Total | Film | Average | Standard | efficiency 

cose-t* | forming | solids thickness | count > | deviation | eps/dps 

solution | 
| 
mg ml mg mg/cm? cps / | // 

0. 02 0.02 | 0.02 | 0.000934 30.9 2.3 48.8 
05 05 | .05 . 00234 74.3 2.9 46.9 
. 10 .10 | .10 . 00467 142.2 3.4 44.9 
. 20 . 20 | . 20 . 00934 232.7 5.7 36.7 
. 50 . 50 | . 0 . 0234 499 4.4 31.5 
1.00 1.00 | 1.00 0467 | 750 4.0 23. 7 
2.00 | 2.00 2. 00 0934 | 1103 3.4 17.4 
1. 00 1. 00 ¢ 9.00 0. 420 4 1971 2. 1 6.9 
1. 00 1. 00 17. 00 . 794 2108 2.0 3.9 
1.00 1.00 | 25. 00 1,170 2108 0.6 ref 


a The activity of the D-glucose-t was 0.0856 wc/mg. The “effective area’ of 
the planchet was 21.4 cm?, 

b Average of three planchets, each counted te 10,000 counts. 

e In the last three experiments, the amount of D-glucose-t was constant, and 
the weight of film was increased by use of an Algin solution. 

4 The observed cps and the total dps in the last three experiments were multi- 
plied by the dilution ratio, i.e., by 9, 17, and 25, respectively, to reduce all results 
to a common specific activity. 
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Ficure 1. Counting efficiencu as a function of film thickness. 


be neglected because the only reflected radiation 
that reaches the counter originates in the source 
itself. With sources composed of carbon, hydrogen, 
and oxygen, the scattering factor is of the order of 
5 percent, and can be neglected in obtaining an 
approximate estimate of the true activity. Prior 
workers have shown that the self-absorption factor, 
fs, for the beta radiation is given by the relationship: 


Ay, 1—e7#T ; 
)={= Sr (3) 
Ao pl : 


where yp is the self-absorption coefficient and A? is 
the radiation from a source containing a unit weight 
of the radioactive material with a thickness 7; Ao 
is the limiting value for A7 as 7 approaches zero. 
If the weight of the source is m mg and the area ts 


. m i : 
acm’, T= "4 and eq (3) can be rewritten: 
( 
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When 7 is infinite with respect to beta particles, as 
in the analytical method, 1—e747=1, and eq (4) 
becomes: 

mAr a 


Ay BR 


~t 


When 7 aproaches zero, f, approaches 1, and from 
eq (2) and the definition of Ao: 


2 dps bs 


Cps 
A(2) 
mn fe | MW 


Since mA,y is the cps of a source weighing m milli- 


a Ed To (6) 


grams, eq (5), by substitution for mAyp and Ay 
becomes: 
ad CPSs = 
== - / 6 (74) 
M dps f ade 
( )> te ‘ to 
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Assuming 27 geometry (f,=0.5) and no back- 
scattering (f,—1), values of 1/u can be calculated 
from eq (7) for all of the measurements recorded in 
table 1. However, there is considerable uncertainty 
in the value of a for determinations made with the 
2-in., commercial, flat-bottomed planchets, because 
of the curved edge. Hence, the best value of 1/u is 
that obtained from the measurements for the stand- 
ardized sample G, in the 10-cm? areas. In the first 
experiment of table 1, a@=10.0, cps=50.1, and 
(dps/m) * fg= (0.435 & 0.0856 & 37,000/7.7) * 0.5. 
From these values there is obtained: 


0.0560 mg/cm’. 


The value of 1/u can be used to obtain 7),. (the 
half-value thickness of the radiation) and /P (the 
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range) from the relationships [9, 17] 7),.=0.693/u 
and R=5/u. With the above value of 1/u, 7y.= 
0.039 mg/em? and R=0.28 mg/em*. This value of 
R can be compared with experimental values, 
reported in the literature, of 0.23 mg/cm? [18], 0.5 
mg/cm? [19],70.6 mg/cm? [6], and with the value of 
0.7 mg/cm?’ based on an empirical energy-range 
relationship [9]. The variation in the results arises 
from experimental difficulties and from uncertainty 
in the absolute disintegration rate of the source. 

It is of interest to note that extrapolation of the 
observed counting rates of figure 1 to zero thickness 
gave a Value of approximately 50 percent of the abso- 
lute disintegration rate, based on the NBS tritium 
standard. With a stainless-steel planchet, under 
conditions somewhat comparable to those used here, 
Rydberg’s data [6] for thin films indicate a back- 
scattering factor of about 1.25. Thus the count 
would be expected to be approximately 62 percent of 
the absolute rate, with a stainless-steel planchet, 
2r geometry, and no self-absorption. A low value 
might arise from unequal distribution of the material 
on the planchet. This factor is of no importance 
with thick samples, but could give rise to large errors 
with thin samples. Such errors make accurate assay 
of tritium materials in thin samples impracticable. 

Assay of infinitely thick films is, however, highly 
satisfactory, and is recommended for routine, labora- 
tory analysis. In eq (1), the empirical constant & 
includes corrections for the area of the source, 
geometry of the counter, self-absorption, backscat- 
tering, and conversion of total dps to microcuries. 
By use of eq (1) and (7), # may be related to a/u thus: 


F 7 Value derived from the average of two half-thickness values experimentally 
determined by Dorfman. 
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k =— B 
aX fo fo 37,000 


However, because of uncertainties in the values for 
a, fe, and f,, the constant k, for the calculation of 
absolute radioactivity from observed counts, should 
. experimentally determined under the conditions 
of use. 
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Tritium-Labeled Compounds II. General-Purpose Appar- 
atus, and Procedures for the Preparation, Analysis, 
and Use of Tritium Oxide and Tritium-Labeled 
Lithium Borohydride' 


Horace S. Isbell and Joseph D. Moyer 
(June 11, 1959) 


A general-purpose manifold is deseribed, which is useful for numerous procedures 
involving tritium gas and tritium-labeled materials. 

Methods are giver in detail for (a) converting tritium gas to tritium oxide, (b) preparing 
tritium-labeled lithium borohydride, and (c) conducting a variety of reactions in a closed 
system. Auxiliary equipment is shown, including water traps, reaction flasks, and apparatus 
for preparing solutions and making filtrations in closed systems. Methods are presented for 
assaying tritium-labeled water by (a) dissolving it in a phosphoric anhydride-sulfurie acid 
solution and counting in a windowless, gas-flow, proportional counter, or (b) converting it 
to hydrogen-/ and assaying it in an ionization chamber. The use of tritium-labeled lithium 
borohydride is illustrated by the preparation of p-galactitol-1-t from p-galactose. Subse- 
quent papers will deseribe the use of this apparatus in the synthesis of tritium-labeled 
carbohydrates. The detailed description of techniques and apparatus should be helpful 
to others who are interested in using tritium. 


1. Introduction 


In a program sponsored by the Division of Re- 
search of the Atomic Energy Commission at the 
National Bureau of Standards, efficient methods 
have been developed for the production of C!- 
labeled carbohydrates [1 to 16].2~) Attention has now 
been directed to the development of methods for the 
synthesis, analysis, and study of  trittum-labeled 
materials. As the first step in the program, a con- 
venient method was developed fer the assay of non- 
volatile, water-soluble tritium compounds [17]. Un- 
doubtedly, this and other improved methods of 
analysis [18 to 21] will lead to rapid expansion in the 
use of tritium compounds for research.* 

The production of tritium-labeled carbohydrates 
required the development of apparatus and#tech- 
niques for conducting a variety of operations in closed 
systems. These include: (a) Dilution of tritium 
with hydrogen, (b) removal of samples of gas for 
radioassay, (¢) conversion of tritium to trittum oxide, 
(d) purification of materials by filtration and erystal- 
lization in closed systems, (e) exchange of tritium 
with hydrogen in solids, (f) recovery of hydrogen-t 
released in chemical reactions, and (g) circulation ot 
hydrogen-t, in catalytic reductions. The following 
sections describe convenient equipment for conduct- 
ing these and other procedures. 

Detailed directions are given for the preparation 
of tritium-labeled lithium borohydride. This versa- 
tile reagent makes possible a wide variety of reactions 
for introducing tritium into specific positions in 

1 Project sponsored by the Division of Research of the Atomic Energy Com- 
oo in brackets indicate the literature references at the end of this paper. 
3 Published research on deuterium and tritium compounds prior to 1956 is listed 


in a series of annual bibliographies prepared by the Thermodynamics Section of 
NBS [22] under the auspices of the Atomic Energy Commission. 


organic compounds. Use of tritium-labeled lithium 
borohydride is ilustrated by the preparation of 
p-galactitol-1/-¢* from p-galactose. Other applica- 
tions of the reagent will be described in future publi- 
cations. 

A convenient procedure is described for the r&dio- 
assay of water-t in a proportional counter. The 
water-t is dissolved in a phosphoric anhydride-sul- 
furie acid solution and counted in a windowless, 
gas-flow, proportional counter in a manner similar 
to that used for assaying C-labeled compounds in 
solution [23]. 


2. Apparatus and Procedures 


2.1. General-Purpose Manifold ° 


The general-purpose manifold shown in figure 1 is 
constructed from  high-vacuum-grade _ ball-joints, 
stopeocks, and other parts. In figure 1, the stop- 
cocks in the primary line are numbered; those con- 
necting to auxiliary equipment are lettered. The 
commercial, automatic Toepler pump, 7’, is operated 
by an auxiliary pump activated by a lock-up or latch 
relay. <A fritted-glass filter, just above the check 
valve on the pressure side of the Toepler pump, 
prevents entrainment of particles of mercury. To 
conserve space, the copper oxide furnace (for con- 
verting hydrogen-t to water-t) is placed above and 
behind the main line; it is kept in place, but is by- 
passed when not in use. Interchangeable, low- 
temperature traps (fig. 2) are used for collecting and 
storing water and other volatile liquids. Part (a) 

4 The product is not a meso compound, but a true derivative of p-galactitol. 
5 The manifold and auxiliary equipment are more or less conventional, and ideas 


for their construction were drawn from many sources. Textbooks by Glascock 
[24] and by Sanderson [25] were particularly helpful. 
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General-purpose manifold, 
Sintered-glass filter. 
Automatic Toepler pump. 
Pirani vacuum gauge. 
Mercury manometer. 
Hydrogen-t storage flask. 
Gas buret. 

Tritium in break-seal tube. 
Gas-sampling tube. (See 
Water-storage flask. 
Copper oxide furnace. 
Cold traps. (See fig. 2). 
Break-seal receiving tube. 
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Interchangeable traps. 
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of this trap is ordinarily used alone, but for more 
efficient removal of condensable gases, a multiple 
trap (b) is included [26]. The two components are 
immersed in the same freezing bath, and can be used 
in combination (on the same outlets as a single trap). 
The traps are shown symbolically at (/—/J) and 
(L-M) in figure 1. The gas buret, 2B, is used not 
only for measuring and dispensing hydrogen, but 
also for transferring gas to the reaction flasks. At 
the close of each preparation, hydrogen-t in the sys- 
tem is either converted to water-t in the copper 
oxide furnace (/-G@), or is pumped into flask A 
attached to the manifold; any residual gas is stored 
in the Toepler pump. The arrangement given in 
figure 1 is used for the production of tritiated water. 
Other arrangements are used for special purposes; 
these will be described in subsequent reports _on 
specific products. 


2.2. Dilution of Tritium with Hydrogen 


The apparatus of figure 1 is assembled, made 
vacuum-tight,® and evacuated. The system is then 
filled with dry hydrogen, introduced at /7, and the 
quantity of hydrogen required for the desired dilution 
is drawn into gas buret B. The stopcock to B is 
closed and the system is evacuated to a pressure of 
less than 20 microns of mercury. With all stopcocks 
closed, the ‘‘break-seal’”’ joint connecting the tritium 
sample, (, is ruptured by means of a steel ball oper- 
ated by a magnet. Then, with stopcocks 1, 2, and a 
closed, the tritium and the hydrogen are mixed in 6 
by opening stopcocks 6 and ¢ and raising and lowering 
the leveling bulb containing mercury. Finally, the 
hydrogen-t is transferred to the evacuated reservoir 
A for storage, leaving any desired amount in BP. 
Residual hydrogen-t in the line is drawn into the 
Toepler pump and held for future use. The empty 
tube, (, is removed and can be replaced by other 
apparatus, 





® Stopcocks must be of high quality and carefully lubricated with vacuum 
grease. Ball-joints are sealed with vacuum-grade sealing wax. 


2.3. Radioactivity Assay of Hydrogen-t 


The gas in the manifold can be sampled at any 
desired time by means of the gas-sampling tube 
shown in figure 3(a). The tube is attached to the 
manifold at any available outlet. With the stopcock 
to the manifold closed and stopcock 1 of figure 3(a) 
open, the gas-sampling tube and connection are 
evacuated through stopcock 2 by an auxiliary pump. 
After stopcock 2 has been closed, the sample tube is 
filled with hydrogen-t from the manifold. The 
temperature and pressure are then measured. Next, 
1 is closed (to trap a sample of the gas between 1 
and 2). Residual gas in the line connecting the 
manifold is drawn into the Toepler pump, and the 
stopcock of the manifold is closed. Finally, the 
sample tube is removed from the manifold. 
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Fiaure 3. (a) Gas sampling tube and (b) Tube connected for 
transfer of sample to tonization chamber. 


For radioassay, this gas is transferred to the 
ionization chamber with the apparatus of figure 3(b) 
in the following manner. The ionization chamber, 
A, is evacuated through stopcock 2. Then 2 is 
turned to allow gas from the sample tube to flow 
into the chamber. The line from a cylinder of 
hydrogen to 1 is first evacuated and then filled with 
hydrogen. Hydrogen is allowed to pass slowly, 
through 1 and the sample tube, into the ionization 
chamber, until atmospheric pressure is reached; 
then the stopcock on the ionization chamber is 
closed. The ion current arising from the sample in 





the ionization chamber is measured with a vibrating- 
reed electrometer in the customary manner [21]. 
The amount of tritium in the sample is calculated 
with a factor that is obtained by assaying hydrogen-t 
prepared from the NBS tritium oxide standard or 
from any source of known tritium content. 


2.4. Conversion of Hydrogen-t to Water-t 


To convert hydrogen-t to water-t, the copper 
oxide furnace (/-G) of figure 1 is heated to 500° C. 
The stopcocks are set so that hydrogen-t can pass 
from the gas buret B through the furnace and the 
water trap (J—J) (cooled by a dry-ice bath). The 
Toepler pump is started and the gas is circulated 
until the pressure in the system becomes constant. 
If desired, the water-¢ that collects in the trap can 
be diluted with water distilled from flask FE. 
Finally, the trap containing the water-¢ is isolated 
by closing stopcocks 7 and 7 and opening stopcock 4. 
The water-t is stored in the trap until it is used. 
Any residual hydrogen-t in the system is diluted 
with hydrogen from 6. The mixture is then oxidized 
by passing it through the copper oxide furnace and 
a second water trap (Z-—M) (cooled by a dry-ice 
bath). The gas is again circulated by the Toepler 
pump until the pressure in the system becomes 
constant. Finally, stopcocks Z and m are closed and 
the water-t of low activity is stored in the trap until 
needed. 

Water-t can be distilled from the trap into a re- 
ceiver attached to the manifold at k. Water of high 
activity is collected and sealed in a_ break-seal 
receiving tube; it must be used in a closed system. 
Water of low activity is collected in a receiving 
flask such as A of figure 4. The receiving flask has 
a side arm closed with a rubber cap suitable for the 
addition or withdrawal of liquid by means of a 
hypodermic syringe. Under a suitable hood, sam- 
ples of low-activity water can be withdrawn from 
the flask by a micropipet of the syringe type, and 
used as desired. 
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Figure 4. Adapter and flask for use with injection technique. 
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2.5. Assay of Water-t With a Proportional Counter 


A definite quantity of water-t, less than 0.5 ml 
containing at least 10 we of tritium, is placed in a 
5-ml volumetric flask. The water-¢ is frozen, and 
the flask is partially filled with a solution prepared 
from 175 g of sulfuric acid and 25 g of phosphoric 
anhydride. The mixture is warmed to room tem- 
perature; the volume is adjusted to 5 ml with the 
sulfuric acid solvent, and the solution is mixed. 
Then approximately 1 ml of the solution is trans- 
ferred to a stainless-steel cell and counted in a 
a 27, windowless, gas-flow, proportional counter. 
With the cell and equipment used for assay of carbon- 
14 [23], the counting rate was found to be 5.75 eps 
per we of tritium per ml of solution. This corre- 
sponds to a counting efficiency of only about 0.0155 
percent. The low efficiency limits the method to 
assay of relatively high-activity materials. 

The phosphoric anhydride-sulfuric acid solvent 
‘an be used for assay of a variety of tritium com- 
pounds. To avoid absorption of moisture, the 
solution must be transferred to the counting cell 
quickly, and the counting gas must be dried with 
phosphoric anhydride. 


2.6. Assay of Water-t as Hydrogen-t in an Ionization 
Chamber 


Methods used for gas counting of tritium prior 
to 1958 are reviewed in [19]. In the present study, 
satisfactory results have been obtained with the 
apparatus of figure 5 and the procedure described 
below. This consists of diluting water-t with aqueous 
acid, and generating hydrogen-t by allowing the 
acid to react with metallic magnesium. 

A definite volume, or weight, of water-t, containing 
about 3 ye of tritium, is diluted to 10 ml with 5- 
percent aqueous sulfuric acid and the solution is 
placed in flask A of figure 5. About 1 ¢ of mag- 
nesium turnings is placed in flask B. Then the ioni- 
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Apparatus for preparation and assay of hydrogen-t 


FiIGurReE 5. 


A Dropping funnel containing sample of water-¢ in 5-percent 
aqueous hydrochloric acid. 

B Flask containing 0.1 g of magnesium turnings. 

( Water trap cooled by dry-ice—ethyl alcohol bath. 

D Manometer. 

FE Jonization chamber. 


evacuated through stopcock 2 to a pressure of less 
than 0.1 mm of mercury. After the vacuum line 
has been closed, the tritiated aqueous acid is slowly 
added to B, until sufficient hydrogen-+t has been 
generated to bring the system to atmospheric pres- 
sure. The pressure and temperature are read, and 
stopcock 3 is closed. Then the ionization chamber 
is connected to a vibrating-reed electrometer and 
the drift rate is measured in the customary manner, 
The radioactivity of the sample is evaluated by 
comparison of the drift rate with that in a control 
experiment run with a water sample of known 
tritium content (e.g., NBS Standard Sample 4926), 
The results are adjusted to correct for variation of 
temperature and pressure from standard conditions. 


2.7. Apparatus for Dissolving, Filtering, and Freeze- 
drying Materials in a Closed System ‘ 


The apparatus of figure 6 can be used for a variety 
of operations. Flasks A and FF are connected by a 
tube having a sealed-in, fritted-glass filter, C, and 
are joined to flask / through a 24/40 § joint and a 
large-bore stopcock, Y). In most applications, the 
material to be purified is placed in A, and the solvent 
in #&. Flask F is cooled in liquid nitrogen and the 
system is evacuated through stopcock 1. After stop- 
cock 1 is closed, sufficient solvent is distilled from # 
into A to dissolve the material. Stopeock Dis 
closed, and the unit A—F? is turned so that A is above 
RP. Then /} is cooled, and the solution passes from A 
through the filter to PB. The residue on the filter is 
washed by placing a cork ring and a small quantity of 
powdered dry ice on the top of A; solvent condenses 
in A and ultimately draims back into /. After 
extraction of the material on the filter by this means, 
the solvent in / is transferred to / by distillmg or 
freeze-drving. (For details, see the next section.) 
Such volatile solvents as ether, dioxane, liquid 
ammonia, liquid sulfur dioxide, and water may be 


used. Nonvolatile reagents may be placed in either 
A or J}, and solvents and reagents may be distilled 


into either A or / from F or from sources attached 
at stopcock 2. Solutions may be stirred with a mag- 
netic stirring bar, placed in either A or /} before 
evacuation. The device may be attached to the 
general-purpose manifold at 2, and volatile sub- 
stances may be introduced or removed at will. 


2.8. Preparation of Lithium Borohydride-t 


As shown by Brown, Kaplan, and Wilzbach [28], 
when lithium borohydride is heated with tritium gas 
at about 200° C, tritium-hyvdrogen exchange takes 
place and tritium-labeled lithium borohydride is 
formed. For efficient reaction, the lithium  boro- 
hydride must be in a finely divided condition. Com- 
mercial lithium borohydride can be converted to a 
suitable powder with the apparatus of figure 6. To 
obtain dry ether about 0.2 ¢ of lithium aluminum 
hydride and 100 ml of reagent-grade ether are placed 
in flask # of figure 6. When the evolution of hydro- 

7 An excellent review by R. 8S. Tipson [27] describes a variety of schemes for 


effecting crystallization and filtration in closed systems. Some of these are 
similar to, but not the same as, that described here, 
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Fiaurk 6. Apparatus for preparing solutions of materials in 
a closed system. 
A and B_ § Flasks. 


C Filter construeted from funnel and a 20-mm, sin- 
tered-glass filter disk. 
D 10-mm Right-angle stopeock. 
FE 200-m1 > Flask. 
F Cork ring to hold dry ice in place. 
land 2 Precision-grade stopeocks. 


gen has subsided, 2 g¢ of commercial lithium boro- 
hydride is transferred to A (in a dry box filled with 
nitrogen). Then the apparatus of figure 6- is 
assembled and is connected to a vacuum pump 
through stopcock 1 and a trap cooled by liquid nitro- 
gen. The ether in # is frozen with liquid nitrogen, 
and the system is evacuated. ‘Phe vacuum line is 
closed off, and the liquid-nitrogen bath is moved 
from EF to A. The ether from F distills into A. 
When most of the ether has been transferred, stopcock 
D is closed and A is brought to room temperature. 
The lithium borohydride dissolves; the resulting 
solution is filtered to remove insoluble matter by 
turning the unit, carrving A and /, through 180° 
about the horizontal axis and then cooling /. This 
causes the liquid to drain from A into /. When the 
solution has collected in #2, a cork ring and a few 
grams of dry ice are placed on A. Then ether dis- 
tills from 7? into A and collects above the fritted-glass 
filter. When solvent. sufficient for washing the 
residue on the filter has condensed, the cork ring and 
dry ice are removed and flask #2 is cooled. The 
solvent drains into &, washing the residue on the 
filter. Finally, the solvent is removed from /} by 
opening stopeock D, cooling / with liquid nitrogen, 


and warming /?. Care must be used to avoid 
bumping. When all of the solvent has been trans- 


ferred to E, the last trace of solvent is removed 
through the vacuum line by heating /# (in a silicone 
oil bath) to 200° CC. Finally, stopeock D is closed 
and dry nitrogen is introduced through stopcock 2. 


The resulting product is a friable powder suitable for 
“labeling.”’ 

To make tritium-labeled lithium borohydride by 
tritium exchange [28], flask B, containing finely 
divided lithium borohydride, is connected to a source 
of tritium as shown in figure 7. The apparatus is 
attached to the general-purpose manifold (fig. 1 atc), 
made vacuum tight, and evacuated. Stopcock 1 and 
the stopcock to the manifold are closed, the break- 
seal to the tritium is broken, and flask B containing 
the lithium borohydride is heated in a silicone oil 
bath at 200° C. After 96 hr, the flask is opened to 
the manifold, and the hydrogen-tritium mixture in 
B is transferred by the Toepler pump to a receiver 
attached to the manifold. All traces of tritium are 
removed by introducing hydrogen from the manifold 
and re-evacuating. Finally, dry nitrogen is intro- 
duced into flask B from a source connected to the 
manifold. Then flask B is removed, stoppered, and 
kept in a dry box until the labeled lithium borohy- 
dride is to be used. The specific activity of the 
product depends on the relative proportions of 
tritium and lithium borohydride used. A prepara- 
tion derived from 2 ¢g of lithium borohydride and 4 
curies of tritium gave a product having an activity 
of 2.94 curies. Thus, 73 percent of the tritium 
used was fixed in the hydride. 





FicurE 7. Apparatus for treating solids with tritium gas. 
2.9. Preparation and Use of Solutions of Lithium 
Borohydride-t 


Under nitrogen in a dry box, the desired amount 
of lithium borohydride-t (from flask B, described in 
the preceding section) is placed in flask A of figure 4. 
The flask is closed (in the dry box) with the adapter, 
B. The side arm of A is closed with a rubber 
stopper, J), suitable for perforation with a hypo- 
dermic needle. The flask and adapter are weighed 
before and after introduction of the lithium boro- 
hydride. The flask is then connected, as shown in 
figure 4, to flask C, containing tetrahydrofuran and 
sufficient lithium aluminum hydride to remove all 
moisture and other reactive contaminants. About 
20 ml of tetrahydrofuran is used for each 100 mg 
of lithium borohydride in flask A. The tetrahydro- 
furan in Cis frozen in liquid nitrogen, and the system 
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(including A) is evacuated. Then the tetrahydro- 
furan from C is distilled into A, by warming C and 
cooling A. After the solvent has been transferred 
to A, stopcock 2 is closed and A (including adapter 
B) is disconnected from C. 

Reductions can be performed directly in flask A 
A solution of the material to be reduced is added 
through the stopper D by means of a hypodermic 
syringe fitted with a stainless-steel needle. When 
reduction is complete, the reaction flask is connected 
to the manifold of ‘figure 1 and any hydrogen-t 
present is transferred by the Toepler pump to a 
reservoir attached to the manifold. The nonvolatile 
product is isolated in a manner appropriate to the 
particular substance (see sec. 2.11). In the event 
that several reductions are to be performed, portions 
of the lithium borohydride solution in A are with- 
drawn with a graduated syringe and injected into 
several evacuated reaction flasks (like that shown 
by A of fig. 4) containing the material to be reduced. 


2.10. Analysis of a Solution of Lithium Borohydride-t 


For analysis, a known volume of the solution of 
lithium borohydride-t in tetrahydrofuran is injected 
into 5-percent aqueous hydrochloric acid contained 
in the gas buret of figure 8. The injection is made 
through a rubber cap, D, near the bottom of the 
buret. A small amount of mercury is placed in D, 
to prevent leakage. The hydrogen-t is collected in 
the gas buret, and the radioactivity is measured by 
transferring the gas from the gas buret to an evac- 
uated ionization chamber, with the arrangements 
shown in figure 8. After introduction of the sample, 
the buret is flushed with hydrogen generated by 
injection of a solution of nonradioactive sodium 
borohydride at Finally, the ionization chamber, 
(’, is filled to atmospheric pressure with hydrogen, 
and the radioactivity is measured with an electrom- 
eter in the conventional manner. 

The specific activity of lithium borohydride-¢ can 
be estimated from the specific activity of a product 
made from it, as, for example, p-galactitol-1-t 
described in the next section. This method of 
analysis, although useful, is subject to error from 
isotope effects. 


2.11. Preparation of D-Galactitol-/-t 


A magnetic stirring bar, two millimoles of p- 
galactose, and one millimole of sodium carbonate are 
placed in a 50-ml flask of the type shown in figure 4. 
The flask is connected to the general-purpose mani- 
fold (fig. 1) and evacuated. Then the connection 
to the manifold is closed; the flask is cooled by an 
ice bath, placed on the driving mechanism of the 
stirrer. The stirrer is started, and 2 ml of ice water 
is added by injection with a hypodermic needle, 
through the rubber stopper of the flask. This is 
followed by addition, in the same manner, of 2 ml 


of a solution containing 12 mg of lithium boro- 
hydride-t (1.47 me/mg) in dry _ tetrahydrofuran. 


Stirring is continued for about 1 hr at ice tempera- 
ture, and then the solution is allowed to warm to 
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FiGuRE 8. Apparatus for assay of metal hydrides. 


room temperature. After 18 hr, 1 ml of 5-percent 
aqueous acetic acid is added to decompose any un- 
reacted lithium borohydride; the solution is frozen 
in a liquid-nitrogen bath, and any hydrogen-t present 
is collected in the manifold. Finally, the connection 
to the manifold is closed and the flask is removed 
from the manifold. The solvent is removed by the 
freeze-drying technique, and, ultimately, transferred 
to a closed container used for disposal of radioactive 
waste. The residue containing the product is dis- 
solved in 10 ml of water, and the solution is evape- 
rated to dryness in a rotary still. The residue is 
dissolved in water and the solution is passed through 
a column containing 10 ml of cation-exchange 
resin.’ The effluent is evaporated to dryness on 
a vacuum rotary still; the residue is dissolved in 15 ml 
of methanol and the solvent is evaporated. Dis- 
solution in methanol and evaporation are repeated 
4 times to remove all of the boric acid as methyl 
borate. The residue consists of p-galactitol-/-t and 
a small quantity of unreacted p-galactose. It is re- 
crystallized from 15 ml of hot e ‘thanol. Ina prepara- 
tion conducted as described, the crude product 
contained 16 me of tritium, and the yield of re- 
crystallized p-galactitol-/- was about 90 percent 
of the theroretical. The p-galactitol-/-t was assayed 
for radioactivity by the method of Isbell, Frush and 
Peterson [17]. A 0.00025-mg sample in a sodium 
O-(carboxymethyl) cellulose film weighing 19.8 mg 
gave 126 counts per second. This corresponds to an 
activity of 44 we/meg. 


§ Amberlite IR_120-H, Rohm and Haas Co,, Philadelphia, Pa. 
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Preparation of Trichloride and Tetrachloride 
of Molybdenum 


Dwight E. Couch and Abner Brenner 
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Molybdenum trichloride was successfully prepared in quantity by the reduction of 


molybdenum pentachloride with hydrogen, 


The most satisfactory yields were obtained 


with a 4 to 5 mole excess of hydrogen at pressures of 100 psi or higher and at a temperature 


of 125°C. 


Molybdenum tetrachloride was prepared by direct reaction of molybdenum trichloride 


with molybdenum pentachloride in a sealed tube or steel bomb maintained at 250° C. 


X-ray 


patterns of the various chlorides were obtained. 


1. Introduction 


Blomstrand [1]! prepared low-valence molybdenum 
halides by passing hydrogen over the pentachloride 
ina hot tube. The process has been modified many 
times. Liechti and Kempe [2] used a similar pro- 
cedure and claimed to have prepared molybdenum 
tetrachloride by condensing the vapors resulting 
from the thermal decomposition of the trichloride. 
Michael and Murphy [3] prepared molybdenum tet- 
rachloride and pentachloride by reacting the di- and 
trioxide with carbon tetrachloride at 250° C. Linder, 
Haller, and Helwig [4] prepared molybdenum di- 
chloride by the reaction of molybdenum with phos- 
gene at about 600° C. However, Biltz and Fendius 
[5] claimed that the tetrachloride had never been 
prepared. They prepared the dichloride by dispro- 
portionation of the trichloride. Wardlaw and Webb 
[6] studied molybdenum pentachloride in’ organic 
solvents and found that it reacted with dry prvidine 
to form pyridine-molybdenum tetrachloride, but 
they were not able to remove the pyridine from this 
salt completely. Forland [7] obtained a patent for 
the preparation of molybdenum chlorides by the 
direct chlorination of molybdenum disulfide. 

Hellriegel [8] prepared molybdenum di- and. tri- 
chloride from molybdenum and molybdenum penta- 
chloride. Senderoff and Brenner [9] tried several 
reported methods in an effort to prepare the lower 
molybdenum chlorides. They found that the direct 
reduction of pentachloride by hydrogen [1] was very 
slow and gave only a small vield of molybdenum 
trichloride, and also reported the reaction between 
molybdenum and phosgene to be unsatisfactory. 
They prepared the trichloride by a procedure similar 
to that of Hellriegel [8] and also prepared the di- 
chloride from the trichloride. Senderoff and Labrie 
[10] attempted to prepare molybdenum trichloride 


—— 


Figures in brackets indicate the literature references at the end ofthis paper. 


by reduction of the pentachloride with hydrogen, 
aluminum, and sodium, but failed to develop a suit- 
able procedure. They prepared molybdenum tetra- 
chloride containing organic material by heating 
molybdenum pentachloride with cetane or paraffin 
and obtained the trichloride by heating this product. 

Since the tetrachloride and trichloride were re- 
quired in moderate quantities for electrolytic studies, 
attention was directed toward the development of 


more convenient methods of making these com- 
pounds. As molybdenum pentachloride is now com- 


mercially available, it was used as the starting 
material. The trichloride was readily prepared by 
hydrogenation of the pentachloride under pressure 
at about 125°C. The reaction was unusual in that 
it proceeded to completion only with solid penta- 
chloride. The tetrachloride was prepared by heat- 
ing the trichloride and the pentachloride together in 
a sealed tube. 


2. Preparation of Molybdenum Trichloride 


Ordinary high pressure equipment was used for 
all hydrogen reductions discussed in this paper. It 
consisted of a stainless steel bomb equipped with 
glass liner. The bomb head was equipped with 
a pressure gauge and a three-way valve for evacu- 
ating the bomb and for subsequently adding the 
hydrogen. The bomb was heated with an external 
resistance furnace. All chemicals were handled in 
an inert atmosphere to prevent contamination of 
the molybdenum compounds with moisture and 
oxvgen. After assembly, the bomb was evacuated, 
flushed once with hydrogen, refilled with hydrogen, 
and heated. 

The determination of the yield of molybdenum 
trichloride was accomplished by mixing the reaction 
product with 1:1 hydrochloric acid-water solution 
in which the trichloride is insoluble, then filtering 
and drying the residue in a vecuum desiccator. This 
procedure gave the yield to within 2 percent, an 
accuracy sufficient for this work. 
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A few preliminary reductions of the pentachloride 
were made with an excess of hydrogen at a temper- 
ature of 200° C and a pressure of 5 to 10 atm (see 
table 1). An inspection of the reaction product 
indicated that only the upper surface of the penta- 
chloride had reacted and a solid crust of molybdenum 
trichloride was formed on top of the liquid penta- 
chloride. The reaction was slow because of the 
low rate of diffusion of hydrogen through this crust 
of molybdenum trichloride. Experiments 1 and 
2, table 1, show that the procedures reported in the 
literature are not satisfactory for preparing mo- 
lybdenum trichloride. 


TaBLe 1. Preparation of molybdenum trichloride by reduction 
with h ydrogen 
| H> Wt Mole ratio |, MoC] 

Exp. Temp Time | pres- | MoCls |Wt loss | Ho: MoCls in 
sure used product 

of hr psi g g % 
Ls 125 to 160 4 14.7 50 None 
2. _| 250 to 275 12 14.7 200 5 
3. 200 65 1375 100 25:1 79 
4 145 1 90 211 13.5 0.4:1 4 
5, 135 0.5 1100 50 9.6 36:1 SO) 
6. 125 0.5 100 dO 10.7 4:1 SO 
7 125 16 55 100 12 0.5:1 S 
8 125 16 150 50 12 5:1 98 
9 125 18 450 50 12 16:1 OS 
10) 125 12 $50 50 11 18:1 97 
11 - 85 16 1600 50 0.5 58:1 None 
12 25 16 1750 100 None 31:1 None 


The reduction of pentechloride to trichloride was 
98 percent complete after 16 hr at 125° C with a 5 to 
1 mole ratio of hydrogen to pentachloride. The 
temperatures and pressures were not critical as 
long as there was an appreciable excess of hydrogen 
and the temperature was between 120°C and the 
melting pomt of molybdenum pentachloride, 190° C. 
Results of successive treatments with fresh hydrogen 
indicated that once the molybdenum trichloride had 
formed it was stable at the conditions used, 125° C, 
and was not reduced to the dichloride. 

The last traces of molybdenum pentachloride and 
hydrogen chloride were removed from the trichloride 
by heating it in an inert atmosphere or in a vacuum 
after the reduction. Generally, this was accom- 
plished by heating to 250° C under vacuum for about 


15 min. Chemical analysis of a typical preparation 
gave 46 percent of molybdenum and 52 percent of 
chlorine. (Theoretical: 47.5 percent of molybdenum 


and 52.5 percent of chlorine.) The trichloride pre- 
pared in this way was practically insoluble in aqueous 
hydrochloric acid. However, if it was allowed to 
stand in contact with water for several hours, some 
oxidation took place as shown by the formation of 
a blue molybdenum solution. It also formed 
oxvchlorides if allowed to remain in contact with 
moist air for several weeks. The trichloride was 
insoluble in dimethyl formamide, acetone, pyridine, 
acetic anhydride, ethyl alcohol, ethy] ether, and 
dimethyl cyanamide, but was soluble in ethyl 
pyridinium bromide at 110° C. 


The molybdenum trichloride prepared by hydro- 
genation involved a solid phase reaction, and non- 
volatile impurities remained in the product. It was 
desirable to purify molybdenum trichloride by sub- 
limation to compare the properties of the solid-phase 
reaction product with that deposited from the 
gaseous phase. However, molybdenum trichloride 
decomposes on heating according to the reaction: 


heat 
12 MoCl,; ———> Mo,Cl,2.(s) +6 MoCl,(q) 


and on cooling, 


cool 


2 MoCl,(g) ———> MoCl,(s) + MoCl;(s)[11] 


Three sections of glass tubing connected by 24/40 
standard taper jomts were used for preparing the 
molybdenum trichloride. One end, fitted with an 
inert gas inlet tube, was charged with melybdenum 
trichloride. The tube was then heated to about 
500° C while a slow stream of dry helium passed 
through it. The molybdenum trichloride decom- 
posed to yield a bright yellow residue of dichloride 
and a brown vapor. When allowed to cool, this 
vapor disproportionated to form a dark solid which 
was composed of trichloride and pentachloride. 
After the molybdenum trichloride was decomposed, 
the center section of the tube was heated to 300° C 
and the molybdenum pentachloride distilled into the 
third section. The apparatus was then placed in 
an inert atmosphere chamber and disassembled. — In 
this way it was possible to isolate all three of the 
reaction products, molybdenum di-, tri-, and penta- 
chiorides. 

Excessive heating of the molybdenum dichloride 
did not contaminate the trichloride in the experiment 
since the dichloride also decomposes to the tetra- 
chloride when heated above 700° C. In faet, both 
molybdenum trichloride and pentachloride can be 
prepared by heating the dichloride according to the 
following equations [11]: 


2 MoCl,(s) ——> Mo(s)+ MoCl,(q) 


cool 
2 MoCl,(y) ———> MoCl,(s) + MoCl,;(s). 


3. Preparation of Molybdenum 
Tetrachloride 


The preparation of tetrachloride by the dispro- 
portionation of trichloride, as reported by Liechti 
and Kempe [2], was attempted. Chemical analysis 
of the sublimate from this reaction indicated that 
the product was not molybdenum tetrachloride; 
approximately one-half of it was insoluble in water 
and was found to be molybdenum trichloride. Thus, 
the material they prepared was apparently an equi- 
molar mixture of trichloride and pentachloride. 
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Molybdenum tetrachloride prepared by the authors 
according to the procedure described by Senderoff 
and Labrie [10] was assayed for molybdenum tri- 
chloride. The product prepared with cetane at 
165° C contained 40 percent of material insoluble in 
1:1 HCl, assumed to be trichloride; that prepared 
with paraffin at 145° C contained only 6 percent of 
insoluble material. Both of these preparations con- 
tained organic contaminants. The tetrachloride 
ywrepared according to the procedure described by 
Michael and Murphy [3] contained no detectable 
trichloride, was completely soluble in 1:1 hydro- 
chloric acid, and was slightly soluble in carbon 
tetrachloride. 

An attempt was made to prepare tetrachloride by 
the controlled reduction of the pentachloride with 
hydrogen. With smaller amounts of hydrogen than 
are theoretically necessary to produce molybdenum 
tetrachloride (experiment 4, table 1), the product 
contained 4 percent trichloride, 36 percent of tetra- 
chloride, and 60 percent of pentachloride. A second 
reduction (experiment 7, table 1) using the stoichio- 
metric amounts of hydrogen needed to produce 
tetrachloride gave a mixture of halides containing 
8 percent of trichloride, 74 percent of tetrachloride, 
and 18 percent of pentachloride. It was evident 
from these reactions that the molybdenum tetra- 
chloride prepared in this way would always be con- 
taminated with both trichloride and pentachloride. 

Molybdenum tetrachloride was successfully pre- 
pared by direct reaction of pentachloride and 
trichloride in a sealed tube. This reaction was 
carried out by sealing about 5 g of the finely ground 
and mixed halides in a 14-mm ID gage-glass tube 
filled with helium or argon. All operations were 
performed in an inert atmosphere. The complete- 
ness of the reaction was determined by an assay for 
molybdenum trichloride based on its insolubility in 
1:1 hydrochloric acid. Molybdenum tetrachloride 
and pentachloride are soluble. ‘Table 2 shows the 
results of several different preparations. Experi- 
ments number 3 and 4, table 2, indicate that the 
reaction was more nearly complete at 250° C than 
at 400° C, 


TaBLeE 2, Preparation of molybdenum tetrachloride by direc 
reaction of the trichloride with pentachloride 


Exp. Mole ratio Temp Time MoCl, in 
MoCls/MoCl product 
a hr % 

i. pure MoC]l 250 150 99 
2 : 2.5: 250 70 trace 
3: 2:1 400 | 20 9 
Re a 250 | 20 3 
5 1.5:1 | 150 | 150 14 
i) Rost 250 | 150 2 
7 | ap 9 HO 250 | 70 | ) 
Sa Lk 250 150 | 3 
9 1.2:1 250 | 500 | { 
10) 1:1 250 | 40 10 
11 IE 250 | 70 | 7 
ie i] 250 | 150 | 4 


+ 100-2 batch in a steel bomb, 


The high yield of molybdenum tetrachloride is not 
entirely in agreement with the calculated data 
presented by Quill [11], according to which the tetra- 
chloride is not stable when cooled from temperatures 
above 125° C. It is probable that the pressure in 
the tube and the excess of pentachloride shift the 
equilibrium in a favorable direction. 

Since an excess of molybdenum pentachloride was 
required to convert all the trichloride to tetrachloride 
in a reasonable time, it was necessary to find a 
method of separating the pentachloride from the 
tetrachloride. A large number of anhydrous organic 
solvents were tested, but none was satisfactory. In 
fact, no solvent was found which would dissolve 
either halide in appreciable quantities without re- 
acting with it.? 

The products from the reaction of molybdenum 
pentachloride with molybdenum trichloride were 
subjected to vacuum sublimation at 120°C. The 
results showed that the pentachloride could be 
distilled out without causing any thermal decompo- 
sition of the tetrachloride. If any disproportiona- 
tion was taking place during these sublimations, it 
would have been evident from the formation of the 
trichloride. However, an assay of both the subli- 
mate, which was the pentachloride, and the residue, 
which was the tetrachloride, showed only traces of 
trichloriode. A chemical analysis of the residue 
from the sublimation showed that it contained 38 
percent of molybdenum and 59 percent of chlorine. 
(Theoretical: 40.3 percent of molybdenum and 59.7 
percent of chlorine.) 


4. X-Ray Studies 


Table 3 shows the d spacings and relative inten- 
sities of the molybdenum penta-, tetra-, tri-, and 
dichlorides. These data have not been reported 
previously. The molybdenum tetrachloride pre- 
pared by reacting the trichloride with the penta- 
chloride and that prepared in accordance with 
Michael and Murphy [3] and Senderoff and Labrie 
[10] had similar patterns. Diagrams of the molyb- 
denum trichloride prepared by hydrogen reduction 
and by thermal decomposition were identical. 
;Furthermore, after the trichloride was washed with 
dilute hydrochloric acid and then dried in a vacuum 
desiccator, it still gave the same pattern showing 
that it was not affected by short-time treatments 
with dilute hydrehloric acid, but if allowed to- 
stand in contact with moist air the pattern of tri- 
chloride was changed. The X-ray pattern of molybde- 
num dichloride which had been washed with dilute 
nitric acid or allowed to stand for a long period of 
time and dried, differed from that of the pure di- 
chloride. This was rather surprising since the com- 
pound apparently did not take up moisture and 
showed no visible signs of reacting with the nitric 
acid. 


ee 


2 Roger J. Labrie of the NBS has recently found that titanium tetrachloride 
will dissolve molybdenum pentachloride but not molybdenum tetrachloride. 
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TaBLe 3. X-ray diffraction data on various molybdenum The authors acknowledge the financial support 
halides showing d spacing and relative intensities taken with of the Springfield Armory and the assistance of 
aaa 4 , aera ee “a A. M. Brown, C. H. Schreyer, and W. H. Metzger, 

Jr., in making some of the preparations. 


MoCl,; MoClh; MoClae MoCl; MoCls 
(Continued) 
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Infrared Spectra of Asphalts 


Burton D. Beitchman 


(March 18, 1959) 


A procedure for preparing thin films of airblown asphalts and a spectroscopic study of 
such films throughout the infrared region of 2.5 to 15 microns are described. The relation- 
ship of transmittance of several absorption bands with the durability of the asphalt is dis- 


cussed. Exposure to ultraviolet radiation 


produced changes in transmittance at four 


wavelengths which, in a series of asphalts studied, showed statistically significant correlations 


with durability. 


1. Introduction 

Roofing asphalts very greatly in durability 
depending upon the source of the particular asphalt. 
One method for measuring durability that is cur- 
rently being used is based upon the time required for 
development of cracks over 50 percent of the surface 
of an asphalt film when such a film is exposed to 
natural weathering conditions or to radiation from a 
carbon are in an accelerated weathering test.! 

Stewart ? recently reported an infrared study of 
components from three esphalts before and after 
natural weathering. These asphalts varied in dura- 
bility, and Stewart pointed out chemical properties 
in the components that might be associated with the 
durabilities of the asphalts. 

An infrered study of the whole asphalt appeared 
desirable since the general characteristics of asphalts 
could be readily indicated without resorting to a 
fraetionstion procedure. The use of whole asphalt 
permits easy study of changes effected in the asphalt 
by exposure to heat, ultraviolet, and other forms of 
radiation. 

Schwever * has reported a study of films of asphalt 
residua. These films were too thick (0.1 mm) to 
permit study below 8 yw and he therefore resorted to 
solution techniques to obtain absorption data in this 
region, 

The procedure described in this paper offers a 
simple means for studying whole asphalt. 


2. Procedures 


An infrared study was made of 28 airblown roofing 
asphalts from widely different sources. The soft- 
ening points 4 of these asphalts ranged from 205° F 
to 241° Fand the penetrations at 77° F from 11 to 23. 

Films of asphalt were prepared by placing a 
small drop of molten asphalt on a sheet of cello- 
phane film. A second sheet of cellophane was 


1S. Hi. Greenfeld, ASTM Bull. 193, 50 (1953). 

2 J. E. Stewart, J. Research NBS 58, 265 (1957), R P2759. 
H. E. Schweyer, Anal. Chem., 30, 205 (1958). 

‘ASTM Method 1D36-26, 

5 ASTM Method D5-82. 


placed over the asphalt, and the two sheets were 
then quickly placed in a hydraulic press in which 
the upper platen was maintained at approximately 
90° C1: A pressure of 16,000 psi was applied to 
the sheets and maintained for about 20 sec. The 
sheets were removed from the press, and the thick- 
ness of the asphalt film and two cellophane sheets 
was measured with a dial-thickness gage having a 
foot about 4% in. in diameter with which thickness 
could be estimated to the nearest 0.01 mil. By 
determining the thickness of the two cellophane 
sheets, the thickness of the asphalt film was obtained 
by difference. To ensure greater reproducibility 
of results, areas of nearly uniform thickness were 
predetermined and indicated by marking the bound- 
aries of these areas on the cellophane film. The 
asphalt films were obtained in various thicknesses 


ranging from about 0.5 to 2.5 mils. By allowing 
the asphalt film and two sheets of cellophane to 
soak in distilled water for several minutes, one of 
the cellophane sheets could be removed easily. The 
asphalt film and the remaining cellophane sheet 
were soaked again in distilled water for several 
minutes in order to remove the second cellophane 
film. The asphalt film was mounted on a film sup- 
port made of cardboard or plastic which had a 
rectangular opening, 1.43.5 em. The prede- 
termined area of nearly uniform thickness was 
placed over the opening. After allowing the asphalt 
films to dry, infrared spectra were taken. 

The films prepared by this method were not 
perfectly uniform in thickness and thickness ranges 
have been recorded in this paper rather than average 
values. This must be considered a limitation where 
the transmittance (7)? of the unexposed films have 
been used for quantitative comparisons. The 
changes produced at four wavelengths by exposure 
to ultraviolet radiation did not appear to be de- 
pendent upon film thickness over the range of thick- 
ness employed. Thus, correlations which were 
made of changes in absorption with durability do 
not have this limitation imposed upon them. 

6S. H. Greenfeld reported the use of a press for making films with a range of 
5 to 50 mils thickness, NBS Building Materials and Structures Report 147 (1956). 


7 The symbol 7 has been used for optical transmittance in accordance with the 
publication ASA Z 10.6 (1948) of the American Society of Mechanical Engineers. 
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3. Results 


The infrared spectra of films of asphalts obtained 
from widely different sources showed a marked 
similarity in their absorption characteristics. All 
the asphalts had absorption bands at the same or 
very nearly the same wavelengths and their spectra 
differed primarily in the transmittance values of 
some of these bands. 

The structural groupings which appeared to be 
present in all the asphalts examined are hydroxyl, 
—CH,—, C—CHs3, carbonyl, —(CH2),— and aro- 
matic rings. There were also absorption bands 





which indicated aromatic substituents or fused 
rings. Other weak but reproducible bands were 


noted at 8.66 mw (1,155 em) and 9.71 yw (1,030 
em~'). The structural assignments given to the 
various bands are listed in table 1. These are for 
the most part, the same groupings identified by 
Stewart ? in fractions of three asphalts. 

Figures 1 and 2 show representative spectra of 
two of the asphalt films. Table 1 gives values 
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of log 1/7, determined by the base-line method, for 
absorption bands from the spectra of seven asphalts. 
The less durable asphalts appear to have slightly 
more hydroxyl absorption than the more durable 
asphalts. The C—CH, band intensity at 7.25 yu 
(1,380 em!) seems also to show a relationship to 
durability. Uncertainties in the average thickness 
produce corresponding uncertainties in the trans- 
mittance. In the case of the hydroxyl band at 
2.91 uw (3,430 em7!), the uncertainties in average 
thickness do not invalidate the observed general 
correlation with durability. However, in the case 
of the band at 7.25 » (1,380 em), uncertainties in 
the value of log 1/7 may be sufficient to obscure 
any correlation with durability. 

Duplicate samples of asphalts A and F are included 
in table 1 to demonstrate the reproducibility of the 
spectra. The relation of film thickness to absorption 
intensities appears to follow the Beer-Lambert Law 
as illustrated in figure 3. 

Figures 4 and 5 show spectra of two asphalt films 
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Infrared spectrum of asphalt 1D (1.8 to 2.1 mils thick). 
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Infrared spectrum of asphalt B (1.9 to 2.3 mils thick). 
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TaBLE 1, Log 1/r values obtained from the infrared spectra of asphalt films 





——$__——_ —_ 


Frequency _.--.cm—!__} 3, 430 | 


Wavelength .....----.---.--------------+-+2--- w.-| 2.91] 3.38] 5.88] 6.25 
| | } | 
| | 
| faite | Aro- 
Assigned structural groups ..:|—08 


—_—— 


| Soften- | Penetra- | Dura- | Thickness 
Asphalt | ing | tionat | bility} of film 
point®| 77°F. | 
°F | Days | Mils | 

EE. 231 18 | 23] 1.8to2.0/ 0.065 |  (¢) 0.127 | 0.206 
C s. 230 oe 3 45 1.7 to 1.9 . 050 (¢) -050 | .158 
G 2% 13 51| 1.7to1.9| 1041]  (e) “074 |. 200 
2. 216 | 12 66 | 1.8to2.1 . O19 (¢) . 305 . 240 
A... 220 13 73 1.6to 18 . 030 (¢) . 071 one 
A 220) 13 73| 1.6to1.8| .032| (¢) .071 | . 166 
B. 229 | 14.5 73 | 1.5to1.9] .025 (¢) .074 | .200 
Beas : 221 | 19 75 | 1.7to2.0| .033 (¢) .036 | .192 
e.. 221 | 19 75} 1.7t020| .081 (¢) .038 | . 208 


8 Data obtained from asphalt suppliers. 
b Durability by accelerated weathering. 
e Intense absorption. 
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2,960 | 1,700 | 1,600 


| C—CH; |} C=O | matic 
| | 





























| 1,470 | 1,380 | 1, 155 1,030 | 961 | 866 811 | 750 720 
6.80] 7.25) 8.66) 971 | 10.41 | 11.55 | 12. 33 | 13.33 13. 89 
| | C—O, | | 
| —CHe—| | S=O | Naph- | Aromatic 
C—CH;|C -CH3)|-_---.-_| or SiO | thenic | substitution —(CH2).—- 

\ | 

Log 1/7 values (determined by base-line method) 

z eee Gs ES | | 
0.948 | 0.518 | 0.034 | 0.021 | 0.020 | 0.061 | 0.066 | 0.110 0. 072 
1.040 | 542} 034] 2024] .016 | .039 | 060 | .114 079 

| 1.040} .478| .042/ .026| .019] .096} .068| .116 . 080 
1.102} .476| .044] .030/ .019 | .085| .062] .024 .131 

| 942} 415) 2085 | :012 |) -017 | 078 | 2070 | -117 100 

| .977| .409| .028} .018| .013| .072| .060| .101 . 089 
1. 182 411} .026} .026| .012] .100] .061 | .108 115 
1.068 | .458/} .038/ .020} .020| .083 | .068| .114 - 108 
1. 064 480 | .036| .014| .020/ [082 | 076 | .123 100 











Data obtained by J. Falzone and S§., Ishihara, Research Associates for Asphalt Roofing Industry Bureau, NBS. 
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FIGURE 3. 
— 1.1 to 1.3 mils thick, - 


before and after exposure to radiation from a carbon 
are for 19.5 hr. These spectra show characteristic 
changes that were also observed in other asphalts. 
Table 2 shows some of the changes produced in the 
values of log 1/7 of four asphalts. Changes in the 
values of log 1/7 for bands at 3.38 w (2,960 ¢m™) and 
6.80 uw (1,470 em!) are not shown in table 2 (the 
films used for illustrative purposes in this table were 
too thick to show the changes in these regions). 
All the asphalts increased in carbonyl and hydroxy] 
absorption and the less durable asphalts generally 
showed a greater increase in carbonyl and hydroxy! 
absorption and absorption at 8.66 (1,155 em™) 
and 9.71 mw (1,030 em™) than the more durable 
asphalts. 

Plots of durability values versus changes in log 
1/7 at 5.88 uw (1,700 em7!), 2.91 pw (3,430 em"), 8.66 u 
(1,155 em™'), and 9.71 w (1,030 em) appeared to 
show inverse relationships. Because of the small 
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Infrared spectra of asphalt 1. 


- - 1.8 to 2.0 mils thick. 


magnitude of these changes at the? latter ¥ three 
wavelengths, statistical correlations based on their 
rankings were used rather than the actual values, 
and these were compared with the rankings of 
durabilities. On the basis of values obtained from 
spectra of 11 asphalts, Spearman’s rank correlation 
coefficients were calculated for all four correlations 
of durability with changes in the value of log 1/7 
(table 3). 

These correlation coefficients are significant at the 
1 percent level for 9 degrees of freedom. As for 
the correlation of durability with changes in trans- 
mittance at 5.88 y, a linear relation could be assumed 
to exist between durability and the reciprocal of 
changes in log 1/7. The product-moment corre- 
lation coefficient was calculated to be 0.779, which 
is significant at the 1 percent level for the number 
of samples involved. 
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Changes produced in log 1/7 values of asphalts by ultraviolet exposure 


TABLE 2. 
Frequency -- em! 3,430 | 1, 700 
Wavelength____._.- p..| 2.91 5.88 
Assigned structural groups. Ol C=0 
Soft- Penetra- | Dura- Thickness 
Asphalt ening tion at | bility > of film 
points 77° Fa | 
ae Days Mils 
Q 227 18 35 | 2.2 to2.6 0. 055 0. 262 
2°... - 2.2 to 2.6 . 130 . 565 
im 224 21 38 2.2 to 2.5 O19 . 135 
Re. a ee 2.2 to 2.5 . 034 325 
Ss 22 19 50 2.3 to2.6 . 031 . 153 
Se 2.3 to 2.6 061 . 343 
- he 223 23 72 2. 30 to 2. 56 . O45 . 205 
4 US a 2. 30 to 2. 56 . O50 . 321 


® Data obtained from asphalt suppliers 
b Durability by accelerated weathering. 
¢ Exposed 19.5 hr. 


Data obtained by J. Falzone 


1, 600 1, 380 1 355 1,030 961 S66 | 811 750 720 
6. 25 7. 25 &. 66 Y. 71 10. 41 11.55 12. 32 13. 33 13. 89 
C—O, + nae se 
Aro- SO, Naph- Aromatie substitution 
matic | C—CH } or SiO thenic CIT2)4 
Log 1/7 values (determined by base-line method 

0. 421 0. 036 0. 028 0. OOS 0. 071 0. O84 0. 096 0. 101 
. 558 . 144 . 150 O19 . O66 O76 . 105 . 125 
. 284 0.60 . 026 . 05 . 008 067 . 037 090 004 
. 301 54 . 040 10 . 020 OS4 . O61 105 . 110 
. 362 4 038 . 027 O19 . 072 . 039 O76 O86 
. b54 59 . 030 098 O10 042 . 043 090 099 
377 5S 036 . 032 017 062 . 079 . 109 132 
. 407 5S . 036 OHS 024 O50 066 .119 114 


and G. MeDonald, Research Associates for Aspha!t Roofing Industry Bureau, N.B.S. 
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Correlations of durability with changes in log 1/r 
at certain wavelengths 


TABLD 3: 


Speariman’s 
rank cor- 
relation 

coefficient 


Wavelength (wave 
number) 


5.88 w (1,700 em!) 0. 782 
8.66 w (1,155 emo!) 945 
9.71 w (1,030 em-!) 836 
2.91 uw (3,480 em!) | L815 


4. Discussion of Results and Conclusions 


A simple technique for studying the structure of 
whole asphalts by infrared spectroscopy has been 
developed. The use of this method on 28 asphalts, 
which come from widely different sources and which 
vary greatly in durability as roof coatings, showed 
the asphalts to be remarkably similar in absorption 
characteristics. In the unexposed asphalt the only 
absorption band that showed some degree of corre- 
lation with durability was that produced by hydroxyl 
absorption. In general, the less durable asphalts 
exhibited stronger hydroxyl absorption than the 
more durable asphalts. The intensities of the band 
at 7.25 u (1,380 em7!) attributed to C-CH, absorption 
appeared to show some correlation with durability, 
but because of variations in film thickness, uncer- 
tainties existed which tended to obscure the signifi- 
cance of such a correlation. 

The exposure of asphalt films to a carbon are pro- 
duced oxidative changes, which were largely the 
formation of carbonyl and hydroxyl groups and also 
changes in absorption at 9.71 » (1,080 ¢m7') and 8.66 


mw (1,155 em~'). The durabilities of the asphalts 
showed inverse relationships with changes in absorp- 
tion at each of these wavelengths. By calculating 
Spearman’s rank correlation coefficients, each of the 
correlations was shown to be statistically significant. 
A product-moment correlation coefficient was calcu- 
lated for the correlation of durabilities with the 
changes of log 1/7 at 5.88 uw (1,700 em~'). This corre- 
lation coefficient also indicated that the correlation 
was highly significant. 

The absorption band near 9.71 » (1,030 em~') was 
observed in fractions ofasphalt by Stewart.? Stewart 
considered C—O, S=O, or SiO as the most likely 
greups to be responsible for this absorption. He 
also observed that weakening of this band in exposed 
asphaltic fractions appeared to be consistent with 
oxidation of sulfoxides to sulfones. The asphalt 
films exposed to the carbon are in this study, how- 
ever, showed increases in the band near 9.71 uy. 
The correlation of durabilities with changes in car- 
bonyl and hydroxyl absorption and also with changes 
in values of log 1/7 at 8.66 w (1,155 em) and 9.71 u 
(1,030 em~') suggest that the changes at these wave- 
lengths may have been caused by changes in con- 
centration of C—O linkages. 


The author is grateful for the assistance given by 
Victoria E. Sanchez and Dwight G. Moore, Jr., in 
obtaining spectra and thanks members of the Floor, 
Roof, and Wall Coverings Section for their coopera- 
tion. Statistical analyses of the correlations indi- 
cated in this paper were made by Hsien Hsiang Ku. 


WasuHinaton, D.C. (Paper 63A2-—14). 
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Publications of the National Bureau of Standards 


Selected Abstracts 


A coulometric-titration coulometer, 5S. W. 
Smith and J. K. Taylor, J. Research NBS 68C, No. 1, 
65 (1959). 

A highly precise coulometer is described which permits inte- 
gration of currents of 160 coulombs or more with a precision 
of about 1 part in 100,000. The current to be integrated 
oxidizes hydroquinone in an electrolysis cell, producing 
quinone and acid. The quinone is then reduced by constant- 
current coulometric-titration, the end point being indicated 
by pH measurements. 


A refined X-band microwave microcalorimeter, 
G. F. Engin, J., Research NBS 68C, No. 1,77 (1959). 


The microcalorimetrie method for evaluation of the efficiency 
and substitution error of a bolometer mount proposed by 
Macpherson and Ixerns has been the object of further study 
and refinement at the Boulder Laboratories of the National 
Bureau of Standards, and an improved instrument based on 
this technique has been recently placed in operation. The 
new microcalorimeter design features: (1) Greatly improved 
ambient temperature control, permitting higher sensitivity 
and resolution, (2) improved d-e instrumentation, (3) im- 
proved mechanical construction giving better repeatability, 
(4) relocation of the thermopile such that it is no longer 
attached directly to the bolometer mount, thus providing 
flexibilitv in the choice of termination, and (5) a more com- 
prehensive treatment of this calorimetric substitution or 
equivalence error. These features permit the determination 
of the effective efficiency of a bolometer mount to an accuracy 
of better than 0.2 pereent. 


A tilting air-lubricated piston gage for pressures 
below one-half inch of mereury, U. O. Huiion, 
J. Research NBS 68C, No. 1, 47 (1959). 


\ deseription is given of a tilting dead weight piston gage 
constructed at the National Bureau of Standards for ranges 
of differential pressure up to about 0.5 ineh of mereury. A 
resolution of better than one part in a hundred thousand of 
full seale has been obtained by use of the toolmaker’s sine 
bar method of angle measurement. The seale is a linear 
funetion of the sine. The instrument can be calibrated 
from basie measurements of length and weight, is rugged, 
and may be constructed in almost any laboratory mechanical 
shop. Sources of possible errors in reading are discussed in 
detail. Comparative tests with certain other gages or ma- 
nometers are cited wherein linearity was found to be within 
one part in 10,000 and agreement within 2.5 parts in 10,000. 
The uses of the gage are briefly discussed. 


Compact multianvil wedge type high-pressure 
apparatus, EK. C. Llovd, U. O. Hutton, and D. P. 
Johnson, J. Research NBS 63C, No. 1, 59 (1959). 


The paper describes apparatus for generating high pressures 
by applieation of anvil forees against each face of a solid 
polvhedron. The equipment constructed utilizes a tetra- 
hedron of pyrophyllite and four tungsten corbide anvils as 
deseribed earlier by H. T. Hall. External foree is applied 
to one of the anvils and wedge reaction forces act on the 
other anvils, permitting the equipment to be used in a 
conventional hydraulic press. Examples of results of use of 
the equipment are given. Success of the design has been 
shown by the repeated generation of pressures in excess of 
100,000 atmospheres (and higher pressures are anticipated), 
indicating that large existing presses might be used to generate 
pressures of this order in volumes of several cubic inches. 
For higher pressures a two-stage multianvil apparatus is 


proposed in which the second-stage anvils are embedded in 
a large pyrophyllite tetrahedron to obtain necessary support- 
ing forces. 


Comparison of absolute intensities of [OI] 
5577 in the auroral and subauroral zones, F. E. 
Roach, J. W. McCaulley, and C. M. Purdy, 
J. Research NBS 68D, No. 1, 19 (1959). 


The distribution of 5577 zenith intensities is compared for 
Fritz Peak, Colorado, in the subauroral zone and Thule, 
Greenland, near the geomagnetic pole. At the Colorado 
station, the absolute intensity is bright enough to permit 
visual detection during 2 percent of the time; at the Green- 
land station 27 percent of the time. The distribution curves 
have a general similarity, suggesting a phenomenonological 
similarity in the excitation mechanisms at the two stations. 


Compressible turbulent boundary layers with 
heat transfer and pressure gradient in flow 
direction, A. Walz, J. Research NBS 68B, No. 1, 53 
(1959). 


Derivation of a method of approximate calculation using 
integral conditions obtained from basic partial differential 
equations for viscous compressible fluid. Generalization of 
semiempirical laws for turbulent wall friction and dissipation 
for application in the compressible case with heat transfer. 


Creep of cold-drawn nickel, W. D. Jenkins and 
C. R. Johnson, J. Research NBS 68C, No.1, 1 (1959). 


Creep tests were made in tension under constant loads at 
temperatures of 300°, 700°, 900°, and 1200°F on specimens 
of nickel initially cold-drawn 40-percent reduction in area. 
The relations between the experimental results and the 
equations proposed by other investigators are evaluated 
and discussed in some detail. None of these equations was 
found to be suitable to express or predict the creep-test 
results of this investigation with high accuracy as changes 
in structure accompanying the creep processes are inade- 
quately deseribed by these equations. However, conform- 
ance to the parabolic strain-time law was obtained over 
limited ranges of stresses and strains. The effects of cold- 
drawing on the creep properties are discussed, and certain 
etching techniques were employed that gave evidence of the 
presence and movement of dislocations in the structures. 


Current and potential relations for the cathodic 
protection of steel in a high resistivity environ- 
ment, W. J. Schwerdtfeger, J. Research NBS 63C, 
No. 1, 37 (1959). 


In order to evaluate criteria for the cathodic protection of 
bare low carbon steel in a high resistivity environment, 
specimens were exposed in the laboratory for a period of 
two months to a soil having a resistivity of about 20,000 
ohm-em. Emphasis was put on the protective potential 
criterion (—0.77 v_ with reference to a saturated calomel 
half-cell) with and without IR drop caused by the protective 
current included in the potential. A criterion based on cur- 
rent density indicated by the cathodic polarization curve 
was also studied. The current required for protection was 
about three times the magnitude of the corrosion eurrent; 
therefore, the corrosion reaction was either under anodic 
control or an equivalent type of control which was caused 
by high resistance at anodic areas. The best degree of pro- 
tection was achieved on the specimen controlled at —0.77 
v without the IR drop. Applied current indicated by the 
break (change-in-slope) in the eathodie polarization curve 
agreed quite well with the actual current required to main- 
tain polarization at —0.77 v. 


195 








Diffraction of electromagnetic waves by smooth 
obstacles for grazing angles, J. R. Wait and A. 
M. Conda, J. Research NBS 63D, No. 2, 181 (1959). 


The diffraction of electromagnetic waves by a convex cylin- 
drical surface is considered. Attention is confined primarily 
to the region near the light-shadow boundary. The complex- 
integra! representation for the field is utilized to obtain a 
correction to the Kirchoff theory. Numerical results are 
presented which illustrate the influence of surface curvature 
and polarization on the diffraction pattern. Good agreement 
with the experimental results of Bachynski and Neugebauer 
is obtained. The effect of finite conductivity is also con- 
sidered. 


Effect of light and water on the degradation of 
asphalt, L. R. Kleinschmidt and H. R. Snoke, -/. 
Research NBS 68C, No. 1, 31 (1959). 


The effects of light only and of light and water on the degra- 
dation of asphait by accelerated and outdoor exposures were 
determined. Specimens exposed to light only gained weight 
and developed a surface film of material insoluble in common 
asphalt solvents. Specimens exposed to light and water lost 
weight beeause of the leaching out of water-soluble degrada- 
tion products. Since the loss in weight was greater for the 
thicker coatings, it was concluded that a migration of light 
degradable components occurred in specimens exposed to 
light and water. 


Electron beam magnetometer, L. Marion, L. 
B. Leder, J. W. Coleman, and D. C. Schubert, 
J. Research NBS 68C, No. 1, 69 (1959). 


A theoretical investigation of the electron opties of an elee- 
tron beam deflection method for detecting small magnetic 
fields is presented. It is shown that remarkably high sensi- 
tivity can be reached. A laboratory model of sueh a magne- 
tometer was constructed, and it was demonstrated that the 
theoretical estimate of sensitivity 310-5 amperes per 
oersted, could be attained in practice. A discussion of the 
possible improvements which could extend the sensitivity 
of the device is also given. 


Lens design: A new approach, Q. \. Siavroudis, 
J. Research NBS 63B, No. 1, 31 (1959). 


This paper is concerned with the problem of defining and 
calculating the aberrations of an optical svstem and _ their 
application to lens design. A method of defining total 
aberrations is presented which involves the derivatives of 
image space ray coordinates with respect to object space 
ray coordinates. Their calculation makes use of the prop- 
erties of Jacobian matrices and Herzberger’s fundamental 
optical invariant. 


Origin of [Ol] 5577 in the airglow and the 
aurora, F. E. Roach, J. W. MeCaulley, and E. 
Marovich, J. Research NBS 68D, No. 1, 15 (1959). 

The distribution of 5577 zenith intensities at stations in the 
subauroral zone is found to be unimodal with no discon- 
tinuity at the visual threshold. This is interpreted as 
evidence that 5577 airglow and 5577 aurora may have a 
common origin. 

Relations between summation methods and 
integral transformations, W. Greub, J. Research 
NBS 63B, No. 1, 1 (1959). 

Relations between the Lototskvy method of summation and 


those of Borel, Euler, and Knopp are obtained by associating 
an integral transformation with the series transformation. 


Tables of transport integrals: A supplement, 
W. M. Rogers, W. J. Hall, and R. L. Powell, -/. 
Research NBS 68B, No. 1, 23 (1959). 


Tables of values of the transport integrals, defined by 


"© @2"de 
e a qh - a > 
J w) 1, nd 2? 


were published in a previous National Bureau of Standards 
Circular for values of n from 2 through 17. In this report 
values are given for the integrals where n is 18 and 20 and 
where rs ranges from 0.2 to 50.0 in steps of 0.2. 


Transmission and reflection by a parallel wire 
grid, \l. T. Decker, J. Research NBS 68D, No. 1, 87 
(1959). 


A comparison is made at X-band frequencies of the theoretical 
and measured transmission and reflection coefficients for a 
parallel wire grid. The methods used are applicable to the 
measurement of these factors for various building materials. 


Non-Bureau Publications 


An improved 8-hydroxyquinoline method for 
the determination of magnesium oxide in 
portland cement. H. A. Berman, ASTA/ Bull. 
No. 237, 51 (Apr. 1959), 


To determine MgO in portland cement, a rapid method 
(optional in Federal and ASTM. Specifications) involves 
precipitation of the Mg salt of 8-hydroxyquinoline, which is 
then filtered, acidified, and titrated with WBrO,. This 
method gives high results near the specification limit (4-5% 
MgQ), poor operator reproducibility, and excessive specifica- 
tion “failures’” which require time-consuming rechecks by 
the referee ammonium phosphate method. The high results 
are due to coprecipitation of excess 8-hydroxyquinoline with 
the Mg salt at high concentrations. Low results are obtained 
at low reagent or Mg concentrations (0-2% MeO). Pre- 
cipitation conditions were studied. Although they can be 
adjusted to minimize the errors, this approach is not practi- 
enble when the MgO content of the sample is not known in 
advance. However, 2 double precipitation procedure with 
8-hydroxyquinoline has been developed which is as accurate 
and reproducible as the ammonium phosphate method and 
considerably faster. A constant quantity of 8-hydroxy- 
quinoline is added to a definite volume of solution containing 
the Mg (400 ml), the precipitate is filtered, redissolved, and 


reprecipitated with 1 ml. of the reagent solution. The 
improved procedure is proposed as a referee method. After 


separation of SiQs, FesOQ3, AlOs;, MnoOg and CaQ, six cement 
samples can be analyzed for MgO in 2.5 hours. 


Degradation of cellulose acetate films, W. K. 
Wilson and B. W. Forshee, SP Journal 15, No. 2. 
146 (Feb. 1959). 


The degradation of several cellulose acetate film formulations 
has been studied in connection with a project on the develop- 
ment of specifications for a laminating film for the protection 
of archival documents. An unplasticized film and a typical 
plasticized film were studied over a temperature range from 
60° C to 177° C in oxygen and in nitrogen. The degradation 
was principally oxidative in nature and was followed by 
measuring the change in intrinsie viscosity, which was corre- 
lated with changes in the stress-strain diagram. A study of 
the degradation of several commercial films at 124° C showed 
that they varied widely in their stabilitv. For highest 
stability, plasticizers containing an aliphatie ether linkage 
should be avoided and acid acceptors and antioxidants 
should be incorporated into the film formulation. 


Formation of ozone from atomic oxygen at low 
temperatures, R. A. Ruehrwein and J. 5S. Hashman, 
J. Chem. Phys. 30, No. 3, 823 (Mar. 1959). 


Oxygen was passed through an electrodeless discharge and 
was condensed, either alone or mixed with diluent gases, at 
liquid helium temperature. On warming through a tempera- 
ture range of 13-30° Kk. the deposits darkened considerably 
in color (blue) and the amount of ozone formed was deter- 
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mined. From these observations it was concluded that about 
11 percent is the maximum concentration of oxygen atoms 
that can be trapped in the solid and that the reaction 
0+ 0O.= Oy, has an activation energy of 1 keal. At high con- 
centrations of oxygen atoms in the discharged gas, ozone 
was formed upon condensation by the predominant reactions 
(1) OF O=O,% and (2) O+0,=0;. The experimentally ob- 
served conversion to ozone of 58 percent agreed well with 
that calculated from a rate law with Ky= Ko. 


Note on measurement of sine-wave response of 
lenses, R. E. Stephens, J. Opt. Soc. Am. 49, No. 4, 
413 (Apr. 1959). 

In the procedure of Lamberts for measuring the sine-wave 
response of a lens the phase angle is measured by determining 
its cosine. For angles less than 45 degrees the precision 
achieved by measuring the cosine is less than is possible by 
measuring the sine. By displacing the scanning slit one- 
quarter wavelength transversely data for computing the sine 
of the phase angle may be obtained. 


Selection of glasses for three-color achromats, 
R. E. Siephens, J. Opt. Soc. Am. 49, No. 4, 898 (Apr. 
1959). 

If two thin lens elements of two different glasses are com- 
bined into a thin two-element system, then the chromatic 
properties of the system may be deseribed by the effective 
values of VI[V=(N—1)/(Nr—Ne)], and P[P=(AN/Nr—No)). 
The effective value, P, is a linear function of V2 On a graph 
of P versus V, the straight line drawn through the points for 
two glasses is the locus of all possible effective values of V 
and P for combinations of these two glasses. By means of 
such a graph the effective value P for a combination of two 


Other NBS 


Journal of Research, Section 63B, No. 1, July— 
September 1959. 75 cents. 


Relations between summation methods and integral 
transformations. Werner Greub. (See above ab- 
stracts). 

On a modification of Watson’s lemma. F. Ober- 
hettinger. 

Principal submatrices of a full-rowed non-negative 
matrix. K. Goldberg. 

Zeros of certain polynominals. A. J. Goldman. 

Tables of transport integrals: A) supplement. 
William M. Rogers, William J. Hall, and Robert 
LL. Powell. (See above abstracts). 

Lens design: A] new approach. Orestes N. Stav- 
roudis. (See above abstracts). 

Analytical integration of the differential equation for 
water storage. Vujica M. Yevdjevich. 

Compressible turbulent boundary layers with heat 
transfer and pressure gradient in flow direction. 
Alfred Walz. (See above abstracts). 


Journal of Research, Section 63C, No. 1, July— 
September 1959. 75 cents. 


Creep of cold-drawn nickel. William D. Jenkins and 
Carl R. Johnson. (See above abstracts). 

Friction and endurance of prelubricated and unlubri- 
eated ball bearings at high speeds and extreme 
temperatures. Hobart S. White. 


glasses may be matched to P of a third glass whose V differs 
from V by a useful amount. An achromatic lens system of 
the three glasses may then be computed, which will have a 
common focal length for three wavelengths. 


Some radiological applications of gamma ray 
transport theory, M. J. Berger and L. V. Spencer, 
Radiation Research 10, No. 5, 552 (May 1959). 


Gamma ray transport theory is applied to the calculation of 
energy dissipation in bounded water media under conditions 
approximating those found in radiological problems. A com- 
bination of analytical, numerical and random sampling 
techniques is used. Two typical situations are considered: 
a medium in which an isotropic source of radiation is dis- 
tributed, and a medium irradiated by an external source 
consisting of a narrow pencil-beam. The spatial pattern of 
energy dissipation is worked out in detail for sources emitting 
1.28-Mev photons. 


Thermal decomposition of poly(vinyl chloride), 
R. R. Stromberg, S. Straus, and B.G. Achhammer, 
J. Polymer Sci. 35, No. 129, 355 (Mar. 1959). 
Poly(vinyl chloride) polymerized using gamma radiation, 
benzoyl peroxide, and 2,2’-azo-bis(isobutyronitrile) as initia- 
tors, was thermally decomposed in vacuo. The gaseous decom- 
position products as analyzed by mass spectrometry consist 
primarily of hydrogen chloride for decomposition tempera- 
tures below approximately 300° C. At temperatures approxi- 
mating 400° C various hydrocarbons are evolved along with 
the hydrogen chloride. Energetic considerations are dis- 
cussed and a temperature dependent activation energy is 
described. A free radical mechanism for the dehydrochlori- 
nation reaction that appears consistent with the reaction 
rate data is postulated. 


Publications 


Effect of light and water on the degradation of 
asphalt. L. R. Kleinschmidt and H. R. Snoke. 
(See above abstracts). 

Current and potential relations for the cathodic 
protection of steel in a high resistivity environ- 
ment. W. J. Schwerdtfeger. (See above ab- 
stracts). 

A tilting air-lubricated piston gage for pressures 
below one-half inch of mereury. U. O. Hutton. 
(See above abstracts). 

Compact multianvil wedge type high-pressure ap- 
paratus. E. C. Lloyd, U. O. Hutton, and D. P. 
Johnson. (See above abstracts). 

A coulometric-titration coulometer. Stanley W. 
Smith and John K. Taylor. (See above abstracts). 
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Leder, J. W. Coleman, and D. C. Schubert. (See 
above abstracts). 

A refined X-band microwave microcalorimeter. 
Glenn F. Engen. (See above abstracts). 


Journal of Research, Section 63D, No. 1, July— 

August 1959. 70 cents. 

Preliminary results of the National Bureau of 
Standards radio and ionospheric observations dur- 
ing the International Geophysical Year. David 
M. Gates. 

Origin of [OI] 5577 in the airglow and the aurora. 
Franklin E. Roach, James W. McCaulley, and 
Edward Marovich. (See above abstracts). 


197 








Comparison of absolute intensities of [OI] 5577 in 
the auroral and subauroral zones. F. E. Roach, 
J. W. McCaulley, and C. M. Purdy. (See above 
abstracts). 


Origin of ‘‘very low-frequency emissions.” R. M. 
Gallet and R. A. Helliwell. 
Climatology of ground-based radio ducts. Bradford 


R. Bean. 

Power requirements and choice of an optimum 
frequency for a worldwide standard-frequency 
broadcasting station. A. D. Watt and R. W. 
Plush. 

Measurments of phase stability over a low-level 
tropospheric path. M. C. Thompson, Jr. and 
H. B. Janes. 

System loss in radio wave propagation. 
Norton. 

Mode expansion in the low-frequency range for 
propagation through a curved startified atmos- 
phere. H. Bremmer. 

Transmission and reflection by a parallel wire grid. 
Martin T. Decker. (See above abstracts). 

Synoptic variation of the radio refractive index. 
B. R. Bean and L. P. Riggs. 

Low-frequency propagation paths in arctic areas. 


A. D. Watt, E. L. Maxwell, and E. H. Whelan. 


Journal of Research, Section D. Radio Propagation, 
63D, No. 2, September-October 1959. 70 cents. 


Stratification in the lower ionosphere. C. Ellvett 
and J. F. Watts. 

Effect of small irregularities on the constitutive re- 
lations for the ionosphere. kK. C. Budden. 

lonospheric investigations using the sweep-frequency 
pulse technique at oblique incidence. V. Agy and 
K. Davis. 

Fields in electrically short ground systems: an experi- 
mental study. A. N. Smith and T. E. Devaney. 

Diffraction of electromagnetic waves by smooth 
obstacles for grazing angles. J. R. Wait and A. 
M. Conda. (See above abstracts). 

Very-low-frequency radiation spectra of lightning 
discharges. W. L. Taylor and A. G. Jean. 

Radio-wave scattering by tropospheric irregularities. 
A. D. Wheelon. 

Study at 1046 megacycles per second of the reflection 
coefficient of irregular terrain at grazing angles. 
R. E. MeGavin and L. J. Maloney. 

Synoptic study of the vertical distribution of the 
radio refractive index. B. R. Bean, L. P. Riggs, 
and J. D. Horn. 


Kenneth A. 


Analytical integration of the differential equa- 
tion for water storage, V. M. Yevdjevich, J. 
Research NBS 63B, 43 (1959). 


Climatology of ground-based radio ducts, B. R. 
Bean, J. Research NBS 68D, 29 (1959). 


Effect of small irregularities on the constitutive 
relations for the ionosphere, Kk. G. Budden, 
J. Research NBS 68D, 135 (1959). 


Fields in electrically short ground systems: an 
experimental study, A. N. Smith and T. E. 
Devaney, J. Research NBS 68D, 175 (1959). 


Fractional factorial experiment designs for 
factors at three levels, W. S. Connor and M. 
Zelen, NBS Applied Math. Series 54 (May 1959). 


Tables of osculatory interpolation coefficients, 
H. E. Salzer, NBS Applied Math. Series (May 1959), 


Friction and endurance of prelubricated and 
unlubricated ball bearings at high speeds and 
extreme temperatures, H. S. White, J. Research 
NBS 68C, 19 (1959). 


lonospheric investigations using the sweep- 
frequency pulse technique at oblique incidence, 
V. Agy and Kk. Davis, J. Research NBS 63D, 
151 (1959), 


Low-frequency propagation paths in arctic 
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